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Abstract 

We describe the main building blocks of a generic automated package for the calculation of 
Feynman diagrams. These blocks include the generation and creation of a model file, the graph 
generation, the symbolic calculation at an intermediate level of the Dirac and tensor algebra, 
implementation of the loop integrals, the generation of the matrix elements or helicity ampli- 
tudes, methods for the phase space integrations and eventually the event generation. The report 
focuses on the fully automated systems for the calculation of physical processes based on the 
experience in developing GRACE-loop which is a general purpose code applicable to one-loop cor- 
rections in the Standard Model. As such, a detailed description of the renormalisation procedure 
in the Standard Model is given emphasizing the central role played by the non-linear gauge fix- 
ing conditions for the construction of such automated codes. These new gauge- fixing conditions 
are used as a very efficient means to check the results of large scale automated computations 
in the Standard Model. Their need is better appreciated when it comes to devising efficient 
and powerful algorithms for the reduction of the tensorial structures of the loop integrals and 
the reduction of the > 4 point-function to lower rank integrals. A new technique for these 
reduction algorithms is described. Explicit formulae for all two-point functions in a generalised 
non-linear gauge are given, together with the complete set of counterterms. We also show how 
infrared divergences are dealt with in the system. We give a comprehensive presentation of 
some systematic test-runs which have been performed at the one-loop level for a wide variety 
of two-to-two processes to show the validity of the gauge check. These cover fermion-fermion 
scattering, gauge boson scattering into fermions, gauge bosons and Higgs bosons scattering 
processes. Comparisons with existing results on some one-loop computation in the Standard 
Model show excellent agreement. These include e^e^ tt,W^W^ , ZH; 77 tt,W^W~ ; 
ej — > eZjvW and W^W~ W^W~. We also briefly recount some recent development con- 
cerning the calculation of one-loop corrections to 3 body final states cross sections in e'^e~ with 
the help of an automated system. 
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1 Introduction 



1.1 The need for automation in the Standard Model 

Much of the success of the Standard Model, SAi , of the electroweak mteraction rests on the 
results of the various precision measurements, notably those of LEP and SLC. These precision 
measurements required the knowledge of higher order quantum corrections. Although the latter 
are rather involved, calculations are still under control since the bulk of the observables pertain 
to two-body final states. In fact due to the present available energy, the most precise predictions 
relate to fermion pair production, a calculation which is far easier to handle than that for W 
pair production even if one leaves out the fact that for the latter one needs a full 4-fermion final 
state calculation. Next generation machines will involve much higher energies and luminosities 
opening up the thresholds for multiparticle production and/or the need to go beyond one and 
two-loop radiative corrections. On the other hand even when one only considers three particles 
in the final state, the complexity increases tremendously especially within the electroweak frame- 
work. So much so that even a process like e~^e~ — > VePeH which would be the main production 
mechanism for the Higgs at the next linear collider and where the tree-level calculation receives 
a contribution from only a single (non-resonant) diagram, a full one-loop calculation has only 
very recently been completed^ [21 E]- For such processes, hand calculations become quickly 
intractable and very much prone to error. Moreover, a complete hand calculation for such pro- 
cesses is not possible, even for the tree-level cross sections, as one has to resort to numerical 
methods for the phase space integration. Especially for QCD processes, to alleviate some of the 
major hurdles in the calculation of matrix elements for physical observables beyond leading and 
next-to-leading order, one has devised some powerful alternatives to the standard diagrammatic 
Feynman approach jlj, with most recently the development of the twistor-spacejSlEllZllHl- How- 
ever most of them, if not all, involve at most one massive particle and a single parameter, the 
QCD coupling constant. Moreover the techniques work because of the exact gauge symmetry 
of QCD and thus, apart from a handful processes, these methods can not be carried over to the 
electroweak theory where the computations involve a variety of masses and scales. Faced with 
these difficulties the need for computers is even more evident for the calculation of electroweak 
processes. 

Ideally one would like to automatise the complete process of calculating radiative correc- 
tions and multi-particle production starting from the Lagrangian or the Feynman rules to the 
cross section. Automation is, in principle, feasible since most of the ingredients of perturbation 
theory are based on well established algorithms. With the increase in computer power and stor- 
age, together with possible parallelization, one could deal, in a relatively short time, with more 
and more complex projects thus bypassing the problem of huge output files being produced, 
at least in the intermediate stages. The idea of automation of the calculations in high-energy 
physics is not new. It dates back to the 1960's when packages such as SCHOONSCHIPl^ and then 
REDUCE [T?H ITT] had b een developed. These are symbolic manipulation codes that automatise the 
algebraic parts of a matrix element evaluation, like traces on Dirac matrices and contraction of 
Lorentz indices. Such codes have evolved quite a lot with applications not only optimised for 
high-energy physics like FORMAT?' but also more general purpose programs like Mathematical^ 
and Maple |14,. Generation of QED Feynman graphs at any order in the coupling constant 
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was automatised in the late 70'Sjl5,. One of the first major application of these early devel- 
opments in this field was the calculation of the anomalous magnetic moments of the electron 
and the muon^Bj. The first automatic system incorporating all the steps for the calculation 
of a cross section, from Feynman graph generation, amplitude generation through a REDUCE 
source code that produces a FORTRAN code, phase space integration and event generation 
with BASES/SPRING[T7| is GRAND[TH1. It was hmited to tree-level processes in QED. In the early 
nineties, a few groups started to develop packages aiming at the automatisation in the SA4 |19j . 

1.2 Different approaches and levels of automation 

A hand calculation of a process at a certain order in perturbation can follow different methods, 
approaches, tricks and sometimes relies on approximations. It is no wonder then that these same 
issues and variations have translated into the automation of these calculations and have led to 
the implementation of a few softwares with varying degrees of automation, different domain of 
application while exploiting different programming languages, see Ref. PHI for a survey of some 
of these systems. Some of the systems are collections of software tools used to help only certain 
aspects of the hand calculation. Example are codes that only generate the Feynman diagrams, 
like QGRAppj or codes for the algebraic and analytic manipulations on loop diagrams but on a 
diagram by diagram basis like XLOOPS[.22'. Others are designed for specific applications |23[ l2l] . 
like QCD corrections to some electroweak processes for example. The report will concentrate 
only on the fully automatised systems that are able to output a source code for the numerical 
calculation of cross sections without any intervention by the user, apart of course from providing 
the input which consists in specifying the process. In reviewing the characteristics of these codes 
and the different steps that go into building these tools, we will see that some of the specialised 
codes we have just mentioned could be considered as a module in the long chain that goes 
from the Lagrangian to the cross section. To go into the details of how the various steps are 
implemented we will have to be more specific, since there is hardly any standardisation of either 
the methods, the algorithms or the computer language. Therefore we will most of the time refer 
to the experience we gained in developing GRACE. We therefore present the case of automation in 
the SM at the tree-level and concentrate more on the one-loop level. A fully automatic system 
beyond one-loop has not been constructed yet. 

1.2.1 Automation at tree- level 

In the usual diagrammatic Feynman approach followed by most of the automated systems, the 
cross sections can be obtained by computing directly the unpolarised squared matrix elements or 
in terms of the (helicity) amplitudes using spinor techniques. Although the computer algorithms 
for these two techniques can be quite different, in both cases one needs the Feynman rules 
and Feynman graphs to be generated. The automatic systems, GRACE ^25 , CompHEP'26", the 
FeynArts-FeynCalc-FormCalc package (771 OHl 0^ IHUll^. Madgraph^ and fdc SJ,, follow 
the diagrammatic approach with applications to both the SM and its super symmetric version. 
A detailed description of the Minimal Supersymmetric Standard Model, MSSA4, of GRACE is 
found in 

It is also possible to make do without Feynman diagrams and arrive even more directly at the 
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cross sections. This can be achieved through an iterative solution of the equations of motion 
or by solving the Dyson-Schwinger equation[35 l l36l l37|. This approach leads to faster codes. 
However, it has not been extended beyond tree-level. Two |25ll29j of the diagrammatic approach 
codes have on the other hand been extended to one-loop. 

Depending on the method in the diagrammatic approach, Dirac matrices and spinors can be 
treated as symbolic objects before being converted to numerical quantities with the help of a 
symbolic manipulation system as will be explained in section IIIIL One can also take a more 
"numerical" approach where these objects are combined into numerical entities. This is treated 
in section Kill 

Even at tree-level, one problem is the size of the output file when one is dealing with multi-leg 
processes. Automatic systems produce, for a complicated process, a huge source code^ which 
in turn requires a large amount of CPU time for the numerical evaluation. This necessitates a 
large disk space, a human control over a large number of batch jobs or could even necessitate 
to split the source code into small pieces in order to make compilers work. In addition to the 
problem of size, it may be necessary to write specific kinematics routines when the amplitude 
has a complicated structure of singularities, see section Of course, the non-diagrammatic 
approach also requires a proper phase space integration routine. 

1.2.2 Automation at one-loop level 

The problem with the size of the output files and the integration over phase space are exacer- 
bated for one-loop processes. These are however not the major hurdles for extending a tree-level 
code to one- loop. One first needs to master all the theoretical background related to the renor- 
malisation of a model or a theory. A consistent renormalisation procedure that gives all possible 
counterterms, which would have to be implemented in the automatic code to tackle any one-loop 
process, needs to be clearly defined. A symbolic treatment of space-time dimension is inevitable 
for the regularisation of ultraviolet divergences. Infrared divergences will have to be regularised 
either by a small fictitious mass or through dimensional regularisation, DR, [381 139j . A major 
investment has to do with a fast and efficient algorithm for the loop integrations, especially the 
algorithm for the reduction of the tensorial structures to the scalar A^-point functions and the 
reduction of the A^ > 4 scalar functions to lower A^-point functions for codes that allow multi-leg 
one-loop integrals. Here also, almost each code reverts to a specific technique and algorithm. 
For one-loop calculation in the electroweak theory, 2 — > 2 processes are now easily and fully 
automatised as will be made evident in the report. Although there is a large number of Feynman 
graphs for 2 — > 3 processes, automatic systems such as the package FeynArts-FeynCalc-FormCalc 
with the extension j4()j of the one-loop library LoopToolsj41jj or GRACE-loop have shown the fea- 
sibility of an automatic one-loop calculation for 2^3 processes in the SA4, where human inter- 
vention is kept to management of the large number of files and batch jobs. Most important pro- 
cesses for Higgs production at the linear collider, e+e~ — > Ve^eH 0121131) e+e" — > e~^e~H [1^ . 
e+e" ^ ZHH El, e+e" ^ tiH gSlilHlilZl, 77 ^ tiH^ as well as e+e" ^ e+e-7 ^ 
have been computed thanks to the automatic systems. This is also the case of the most re- 

^The size of the problem grows rapidly as the number of external particles increases. For example, 
2 ^ N tree process in (p^ model has {2N — 1)!! Feynman graphs. 
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cent calculations of the complete one- loop electroweak corrections to 2 ^ 4 processes, namely 
e^e~ — > uuHH^Ol and some specific 4-fermion channels in e'^e~ 4 f |51p. 

1.2.3 Checking the correctness of the results 

An automatic system produces some numbers as the results of a black-box calculation which 
may not necessarily be correct. A user may feed some input data which the authors of the 
system have not thought of. There may be bugs in the program which have not been detected 
with the tests made by the developers. Compilers may have problems, especially with highly 
optimizing options. Even if the generated program is logically correct, numerical cancellations 
may produce wrong results. The Monte-Carlo integration package may give some false value just 
because of insufficient sampling points. Thus systematic procedures of checking the results are 
indispensable for automatic systems. These procedures will be classified into three categories: 

i) Checks by the computer system 

If the generated code is set-up so that it can be run with different accuracies (double, 
quadruple precision), it would be easy to detect problems related to numerical cancella- 
tions. Good FORTRAN compilers supply options for changing precisions without modifi- 
cation of the source code. It is also a good idea to run the program on other machines with 
different compilers or architecture. For the Monte-Carlo integration, one can increase the 
number of sampling points in order to test the stability of the results. 

ii) Self-consistency checks within the automatic system 

If the theoretical model has a free parameter, related to its symmetry, which does not 
alter the physical results, it will be used to check the results. Physical quantities in gauge 
theories are independent under the change of gauge parameters. When an automatic 
system includes gauge parameters as variables in the generated code of the numerical 
calculation, one can explicitly check the gauge invariance of the obtained results. When 
the system keeps the regularisation parameters of the ultraviolet or infrared divergences, 
one can confirm the cancellation of these divergences explicitly. The incorporation of 
these self-checking procedures is one of the most important feature for the reliability of 
the automatic calculation. 

iii) Comparison with other calculations 

This is a standard procedure provided another independent calculation exists or can be 
performed using a different automatic code. 

Although the ultraviolet and infrared tests are rather straightforward to implement, the 
gauge parameter check requires a very careful and judicious choice of the gauge-fixing function. 
This is especially true in automatic codes for one-loop amplitudes and cross sections. A few 
tree-level codes have the gauge parameter check incorporated through the usual linear 't Hooft- 
Feynman gauges or give the possibility to switch to the unitary gauge. For example, in GRACE and 
for tree-level processes, the gauge-parameter independence check has been applied successfully 
by comparing the results in the unitary gauge to those of the 't Hooft-Feynman gauge. An 

first investigation of these processes using an automatic code was done in|52|. 
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agreement up to ~ 15(~ 30) digits in double (quadruple) precision has been reached for a few 
selected points in phase space therefore confirming, at this stage, that the system works very 
well for tree-level processes. However, as we will see, these types of gauges are not suited at 
all for one-loop calculations in the electroweak sector and explain, in part, why there are at 
the moment only two general purpose codes for one-loop calculations. None of them exploits or 
is defined for a general linear gauge. The latter tends to considerably increase an already 
very large size of the file corresponding to each of the numerous one-loop diagrams. Not only 
the expressions get large compared to the usual ^ = 1 't Hooft-Feynman gauge but also call for 
extending the algorithms for the reduction of some new tensor and A^-point function integrals. 
These are at the heart of a one-loop calculation especially that their evaluation is very much 
time consuming. Generalised non- linear gauges [53l I54j. still defined with the gauge parameter 
^ = 1, are on the other hand very well suited for an automatic code of one-loop amplitudes. 
This will become clear when we will go through the different stages of a one-loop calculation 
and the different modules that are required for the construction, or the extension to, one-loop 
amplitudes in the electroweak theory. The implementation of the generalised non-linear gauge is 
therefore crucial in GRACE at one- loop. This is also the reason it has a quite central place in this 
report and deserves that we summarise, already at this stage, some of its salient features and the 
simplifications it brings when implemented in an automatic code for one-loop processes. Note 
that the package FeynArts-FeynCalc has the SM defined in the background-field gauge 
beside the usual linear ^ = 1 Feynman gauge. 

1.3 Importance of judicious gauge- fixing for automated one- 
loop calculations 

A computation in a general gauge or unitary gauge brings about unnecessary complications 
and sometimes troublesome numerical unstabilities especially when one deals with several gauge 
bosons. Take for instance the propagator of the W gauge boson, of momentum k and mass M^y. 
In a general R^ gauge, with the gauge fixing parameter ^w, it writes 



In the 't Hooft-Feynman gauge, = 1; only the "transverse" part consisting of the metric tensor 
Qfiu contributes and leads to a straightforward contraction of neighbouring vertices. Numerical 
instabilities are due to the contribution of the "longitudinal" k^ki, part of the gauge propagators. 
Moreover, the longitudinal tensor structure considerably inflates the size of each intermediate 
result, for example with n intermediate heavy gauge bosons instead of performing n operations 
one performs 2" operations. Since the longitudinal expressions involve momenta, they can 
contribute terms that increase with energy and which require a subtle cancellation among various 
diagrams. A situation which is most acute in the unitary gauge, obtained by formally taking 
^VK — > oo in Eq. 11.11 These problems are of course exacerbated in loop calculations and, as is 
known, calculations and renormalisability itself are arduous if not problematic in the unitary 
gauge[SH]. Within GRACE one of the problems in these gauges (unitary, or general linear type 
gauges) is that the library containing the various loop integrals is designed assuming that the 
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numerator for the propagator of the vector particles is g'^'^ . For instance, the hbrary for the 
three-point vertex functions is implemented with only up-to the third-rank tensor and therefore 
the library applies equally well with fermion loops, gauge loops or a mixture of these if the 
calculation is performed in the 't Hooft-Feynman gauge. In any other gauge one would have, for 
the vertex functions alone, had to deal with a 9th rank tensor! Again this not only creates very 
large expressions but also introduces terms with large superficial divergences that eventually 
need to be canceled extremely precisely between many separate contributions. Fortunately 
one can also work with a class of gauge-fixing conditions that maintain all the advantages 
of the usual 't Hooft-Feynman gauge with exactly the same simple structure for the gauge 
propagators. The new gauge parameters modify some vertices involving the gauge, scalar and 
ghost sector and at the same time introduce new vertices. In fact by judiciously choosing some 
of these parameters the structure of the vertices can get even simpler than with the usual 
linear gauge- fixing conditions. The class of gauges we are referring to exploit non-linear gauge 
fixing conditions |53l I54j . Apart from the possible simplifications that these gauges bring, we 
have considered a generalised class of non-linear gauges so as to perform in an efficient way 
the gauge-parameter independence checks within the GRACE system. Actually the generalised 
gauge we choose depends on 5 parameters j54j. Therefore not only this allows for a wide range of 
checks but since the different parameters affect different parts of the bosonic sector one can check 
different classes of contributions to a single process and thus more easily track down any bug. 
There are other welcome features of these checks. They serve as powerful tools on every step 
of the automated computation, from the correct implementation of the model file which in fact 
can be checked mostly at tree-level to the correct implementation of the tensor integrals. The 
reduction of the latter to scalar integrals is most prone to error. We will show how the tensor 
reduction is carried out in GRACE. The gauge check allows therefore to test that the reduction of 
these integrals into the scalar integrals is implemented properly PJ. Additional tests like those of 
infrared finiteness further verify the scalar integrals. Another advantage of the non- linear gauge 
checks over those that may be attempted within a linear gauge is that on a diagram by 
diagram basis, the gauge-parameter dependence in our checks are polynomials in the non- linear 
gauge parameters whereas in the linear gauge the dependence is buried within logarithms 
and rational functions. We will show how one can exploit this fact for a very powerful gauge 
check. 

1.4 Plan and outline of the review 

The aim of this paper is to describe in some detail the workings of a code for the automatic cal- 
culation of cross sections at one-loop in the SM based on the Feynman diagrammatic approach. 
Though a few codes will be reviewed, the details and algorithms that enter the construction of 
such codes are based primarily on the experience we gained while developing GRACE. The next 
section, section |2l will first give a general overview of the main building blocks of a generic 
automated package for the calculation of Feynman diagrams starting from a Lagrangian down 
to an integrated cross section. A short review of some specialised purpose software is given 
in section 12.11 We then discuss in some detail how the model file is implemented in 12.21 and 
how Feynman graph generation is achieved automatically in 12.31 Amplitude generation at tree- 
level, both through the spinor technique and directly through the squared matrix, is described 
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in section |31 This section will also very briefly present alternative automatic codes that make 
do without Feynman graphs at tree-level. It is followed by the implementation of phase space 
integration and event generation. The extension to a one-loop automatic system is described 
in section |3I The issue of the large size of the generated code especially at one-loop and file 
management is discussed in this section. Section |1] will also emphasise the importance of internal 
self-checks on the correctness of the results of an automatic code, in particular the importance 
of the non-linear gauge. This naturally leads the way to two chapters serving as the theoretical 
background. First, in section [S] quantisation in a non-linear gauge is briefly outlined while a 
detailed description of the renormalisation procedure in the SM within a non-linear gauge is 
exposed in section El Explicit formulae for all two-point functions in a generalised non-linear 
gauge are given, together with the complete set of counterterms. We will also point out some 
of the issues of renormalisation and gauge dependence when dealing with unstable particles and 
the problems with the implementation of the width of these particles in section l631 All of this is 
needed for the implementation of the model file in the SM. The need for a non- linear gauge is 
better appreciated when it comes to devising efficient and powerful algorithms for the reduction 
of the tensorial structures of the loop integrals and the reduction of the > 4 point-function to 
lower rank integrals. A new technique for these reduction algorithms is described in section [Tj 
Section El is devoted to how the ultraviolet and infrared finiteness checks are dealt with in the 
system. In this same section we also give a comprehensive presentation of some systematic test- 
runs which have been performed at the one-loop level for a wide variety of two-to-two processes 
to show the validity of the gauge check. These cover fermion-fermion scattering, gauge boson 
scattering into fermions, gauge bosons and Higgs bosons scattering processes. SectionElis a con- 
vincing testimony of the power of a fully automatic system at one-loop since comparisons with 
existing results on some one- loop computations in the SAi show excellent agreement. These 
include e+e" ^ ti, W+W, ZH; 77 ^ ti, W+W; eZ, uW and W+W- W+W- . The 

final section, section ^1 contains our conclusions. We also provide some detailed appendices. 
In particular we provide the full set of Feynman rules within the generalised non-linear gauge 
as well as the library for the counterterms. Full results for all the self-energy diagrams of all the 
particles in the model including the Goldstone sector is also relegated to an appendix. 



2 Overview of an automatic system: 
GRACE as an example 

The different components and steps that go into the calculation of a cross section, or even the 
corresponding Monte-Carlo event generator, in a code for the automatic evaluation of Feynman 
diagrams are depicted in Fig. ^ Automatic systems not based on the Feynman diagrammatic 
approach have a different architecture, moreover the technique they are based on have not been 
extended beyond the tree-level. We take GRACE-loop as an example. Although not all the 
modules are present in all such codes, especially as concerns the one- loop part, this should give 
an overall view of how such systems work and what inputs are required to make them function. 
Details of the different parts and components of the packages will be reviewed in the next sec- 
tions together with the theoretical background. In the following, the most important modules 
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Figure 1: GRACE System Flow. 
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will be enclosed in boxes for easy reference to the flow chart of Fig. ^ 



Theory 



Particles and Vertices 



The system first requires the implementation and definition of a model through a model file that 
gives the particle content, the parameters as well as the vertices as derived from the Lagrangian. 
The implementation of the model file is described in more details in section 12.21 



User Input 



The user, on the other hand, sets as input the incoming particles and the final particles and 
specifies the order, in perturbation theory, at which the cross section is to be calculated. 



Diagram Generator Diagram Drawer 



GRACE will first generate, through its own routine grc the full set of diagrams (tree and loop as 
well as counterterms for the latter) with the possibility of a Postscript output of the Feynman 
diagrams with the help of the utility diagram-drawer gracef ig. Most of the postscript files for 
the Feynman graphs in this report have been produced by gracef ig. The FeynArts|28j package 
based on Mathematica performs similar steps. Section 12.31 will review some of the issues and 
part of the algorithms for the automatic generation of the Feynman diagrams while appendix |j] 
will give some more technical details about graph generation. 



Matrix Element Generator 



The matrix element generator encodes all the information on the diagrams. For the tree-level 
process the system generates a FORTRAN code which describes the helicity amplitudes using the 
CHANEL library and routines |2Z1- More details about the algorithm together with a worked-out 
example will be presented in section |3I We will only briefly describe the approach that avoids 
Feynman diagrams altogether in section 

For the computations at one-loop, one first generates a symbolic manipulation source code, 
based on REDUCE. 11 or F0RM|12j that writes, for each set of tree and loop diagrams, the in- 
terference term 7^'°°P7^*'^<''=t _ helicity formalism option is also possible here. Only then 



the FORTRAN source code is generated and the cross section computed with the help of the 
loop library and the counterterm library that performs the integration over the Feynman pa- 
rameters and takes into account the counterterm constants. The symbolic manipulation for 
the loop calculation performs a number of important tasks, such as Dirac algebra (taking 
traces for fermion loops) and tensor manipulations in n = 4 — 2e dimension if DR 38t i39j is 
used^, introducing Feynman parametric integrals, shifting loop momenta appropriately etc.. 
FevnCalc/FormCalc|29l 1301 131j carries out a similar function with the help, for the one-loop 
integrals, of the LoopTools library 

The implementation of all these intermediate steps that are necessary at one-loop is discussed 
in more detail in sectional To fully appreciate the issues at stake, it is important also to refer 
to the theory sections on renormalisation, section |21 and also on the algorithm for the loop 
integrals in sectional As already mentioned the latter module is a critical part of an automatic 
code at one-loop. This also explains why almost every code comes with its own solution to the 



•^Some issues related to regularisation and the treatment of 75 are discussed in section ^21 
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problem. The implementation of this part has a bearing on the size and hence performance of 
the system, see section |3 and section ^31. This calls for parallelisation and vectorisation of the 
codes as discussed in section |7| 



Phase space integration 



The integration over phase space is carried via BASES^Tj, a Monte-Carlo integration package. 
One can test the convergence of the integration routine and get cross sections and distributions. 
The simulation and event generation is done through the package SPRINGjEl- GRACE includes 

a number of kinematics routines, through a kinematics database , for processes with up to 6 



particles in the final state. The user can select the appropriate kinematics routine from the 
available library depending on the singular behaviour of the process. This singular behaviour 
can be due to a t-channel photon exchange for example or some other peaking behaviour like 
the crossing of a resonance. Some of the issues that need to be addressed in the code for the 
automatic calculations of cross sections as regards integration over phase space will be described 
in section 13.61 



2.1 Specialised codes and building blocks 

As mentioned earlier, some specialised codes exist that only tackle one of the steps above. 
QGRAF[2J is a very powerful Feynman diagram generator. Most of the codes however are matrix 
element generators that work once a specific diagram is supplied to the code. Examples include 
MINCERi24j, MATADjlSl and SHELL2j2H|, all using FORM as a symbolic language. Though being 
devised for up to three-loop diagrams, their domain of application is limited to graphs with 
a restricted hierarchy of masses and momenta. Xloops treats more general graphs, up to 
two-loop, but again on a diagram-by-diagram basis. GEFICOMj^ can be considered as a master 
program that interfaces QGRAF with MINCER(24j and MATAD(2ni- DIANA[nS] is another Feynman 
diagram analyser with graphics facilities that exploits the power of FORM and is based on QGRAF 
for the generation of the Feynman diagrams. However it does not include, for applications to 
loop calculations, all the ingredients we listed in the diagram generators of the complete auto- 
matic system since it lacks the module for loop integration. For a full review of such packages, 
see HO]. 

The new code aiTALcjHOI on the other hand, could be considered as a fully automated code 
for the calculation of cross sections. Built on DIANA, it goes a step further by adding the 
L00PT00LS|41j library. The code is restricted to the evaluation of 2 — > 2 processes at one- loop 
with only external fermions. It does not include hard bremmstrahlung. 

In a different context, MicrOMEGAs is a code written in C for the automatic evaluation of 
the relic density of dark matter in super symmetry. By default the lightest neutralino is assumed 
to be the lightest super symmetric particle, but the user can set any supersymmetric particle 
(not including Higgses) to be the LSP. The code is built on Gale HEP j62j which generates, while 
running, the subprocesses needed for a given set of AiSSM parameters^. 

■^CalcHEPlHS is an outgrowth of CompHEPf^. 
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2.2 Implementation of the model file 

The model file contains all the information about the theoretical model where the calculation 
is to be performed. In particular the whole set of particles of the model and their interactions 
must be transcribed in a machine readable format. To perform the calculation one then needs 
to define the particles and write down the Feynman rules in terms of all possible vertices needed 
to build up the Feynman graphs. For the electroweak theory this means, among other things, 
writing all the set of Feynman rules as listed in Appendix B. At the loop order, the definition 
requires that one specifies counterterms after having set the renormalisation procedure. This 
means for instance that the loop order of a vertex be also specified. The set of one-loop vertices 
that need to be generated at one- loop in the electroweak theory is listed in Appendix IfI 

Figure 2: An example of a model file in GRACE. 

%======================================= 

Order={ELWK, QCD[qcd]}; 

Version={2,2,0}; 

PPhase=2; 

°/o gauge bosons 

I 

Particle=W-plus ["W+"] ; Antiparticle=W-minus ["W-"] ; 
Gnaine={"W", "W-+", "W"-"}; 

PType=Vector ; Charge=l; Color=l; Mass=ainw; Width=agw; 
PCode=2; KFCode=24; Gauge="wb"; 
Pend; 

•/ 

Particle=Z["ZO"] ; Antiparticle=Particle ; 
Gnaine={"Z-0"}; 

PType=Vector ; Charge=0; Color=l; Mass=ainz; Width=agz; 
PCode=4; KFCode=23; Gauge="zb"; 
Pend; 

Particle=Higgs ["H"] ; Antiparticle=Particle ; 
Gname={"H"}; 

PType=Scalar ; Charge=0; Color=l; Mass=ainh; Width=agh; 
PCode=31; KFCode=25; PSelect="higgs" ; 
Pend; 

Vertex={Higgs , W-plus, W-minus}; ELWK=1; FName=chww; Vend; 

Vertex={Higgs, Z, Z }; ELWK=1; FName=chzz; Vend; 

Vertex={chi-minus , W-plus, Z }; ELWK=1; FName=cwzni; Vend; 



Most automatic systems read model files in which the information on the particles and the 
vertices, through the Feynman rules, are coded manually. Figure [2 shows an example of a 
model file for the GRACE system. At first some options are specified. In this example, the name 
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of the coupling constants are defined by Order= .... Each particle is described through a set 
of properties such as the name of the particle, spin, electric charge, representation of internal 
symmetries, whether it is massive or massless and so on. The interactions of the particles are 
defined by a set of vertices. This example shows the definition of some scalar-vector-vector, 
SVV, vertices which consist of the list of interacting particles, the order of the coupling constants 
and the name of the coupling constants used in the generated code. For a model such as the 
SM or any renormalisable model, all types of SVV vertices have a common Lorentz structure. 
As with other types of vertices FFV,VVV etc., see Appendix B, the Lorentz structure is exploited 
in building up the Feynman amplitude, as will be shown below. 

Although it is easy to code a model file by hand for a simple model such as QED, it is 
not always an easy task for a more complicated model. For example, the MSSM consists 
of more than 80 particles resulting, at tree-level alone, in more than 3, 000 vertices. This is 
the reason why it is more desirable to construct model files automatically from a Lagrangian, 
with a minimum of human intervention. In this case the set of Feynman rules and vertices 
is generated automatically. Using dedicated programming languages, software packages have 
been developed for the automatic generation of Feynman rules. One can cite beside LanHEPjnni 
originally designed to work with CompHEP, the codes included in FDC|32] or gss|;64 . In this case, 
a Lagrangian is usually given in a quite compact and symmetric form. Gauge fixing terms and 
ghost terms are added to it. The latter can even be implemented automatically by first defining 
the symmetries of the theory. For instance, in LanHEP, using the BRST transformation j65| 
(see Appendix A) as done in section EUl one can make the system automatically derive the ghost 
Lagrangian and the corresponding Feynman rules. In the following we sketch some of the steps in 
automatically deriving the Feynman rules from an algebraic implementation of the Lagrangian 
and by applying some simple set of rules. We take LanHEP as an example without going into 
the details of the procedure, the interested reader should consult the manual of LanHEP (63]. 
One first needs to define the particles of the model as shown in Fig. |21 for the bosons of the 
electroweak model. 

Figure 3: Particle description in LanHEP 

vector 

A/A: (photon, gauge), 

Z/Z:('Z boson', mass MZ = 91.1875, gauge), 
'W+'/'W-': CW boson', mass MW = MZ*CW, gauge), 
scalar H/H:(Higgs, mass MH = 115). 

Figure |1] shows how parts of the Higgs Lagrangian are entered once the Higgs doublet pp ( . f 
refers to the Goldstones), and the covariant derivative Dpp~mu~a are defined. The command 
Iterm specifies a term in the Lagrangian. Shifts to introduce the wave function renormalisation 
are performed through the command transform, see Fig. |21 It is also possible in LanHEP to 
introduce a command brst for the BRST transformations . In the example of Fig. IHl by 

acting on the gauge-fixing function, this generates the ghost Lagrangian. 
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Figure 4: Entering the Higgs interaction in LanHEP 

let pp = { -i*'W+.f', (vev(2*MW/EE*SW)+H+i*'Z.f ')/Sqrt2 }, 
PP=anti(pp) . 

Iterm -2*lambda* (pp*anti (pp) -v**2/2) **2 
lainbda=(EE*MH/MW/SW)**2/16, v=2*MW*SW/EE . 

let Dpp~mu~a = (deriv~mu+i*gl/2*B0~mu) *pp~a + 

i*g/2*taupm~a~b~c*WW~mu~c*pp~b . 
let DPP~mu~a = (deriv'mu-i*gl/2*B0"mu) *PP~a 

-i*g/2*taupin~a~b~c*{'W-' "mu.WS'mu, 'W+' ~mu}~c*PP~b. 
Iterm DPP*Dpp. 



Figure 5: Introducing wave function counterterms in LanHEP 

transform A->A*(l+dZAA/2)+dZAZ*Z/2, Z->Z*(l+dZZZ/2)+dZZA*A/2, 
'W+'->'W+'*(l+dZW/2) , 'W-'->'W-'*(l+dZW/2) . 
transform H->H* (l+dZH/2) , ' Z . f ' -> ' Z . f ' * (1+dZZf /2) , 

'W+.f '->'W+.f '*(l+dZWf/2) , 'W-.f '->'W-.f '*(l+dZWf /2) . 



Figure 6: Introducing gauge-fixing and ghosts in LanHEP 

let G_Z = deriv*Z+(MW/CW+EE/SW/CW/2*nle*H)*'Z.f ' . 
Iterm- G_Z**2/2. 
Iterm - ' Z . C ' *brst (G_Z) . 



2.3 Feynman Diagram Generation 

The automatic generation of Feynman diagrams necessary for the computation of a process 
within a model proceeds after the user has defined the process through an input file. Figd shows 
an example of an input file specifying a process. In this example it is e'^'e" W~^W~^ at tree- 
level . The input data specifies a theoretical model (here sm.mdl), the order of perturbation by 
counting the power in the coupling constants (here ELWK=3) and by choosing the initial (e"^, e~) 
and final particles {"f^W~^,W~). In the GRACE input file, Kinem="2302" identifies a choice of 
kinematics in the kinematics library to be used for the integrated cross section^. 

A typical algorithm for the generation of the Feynman diagrams will be the following: 

1. Generate the number of vertices. 

The number of vertices is restricted by the order of the coupling constants for the physical 
process and is given by the input file, see Figd Each vertex has a fixed number of 
propagators and external particles to be connected. 

'^For details about the format of the input file, see Ref |84| . 
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Figure 7: An example of the input file in the electroweak SAi for the scattering process 
e~^e~ — > W^W~'y at tree-level. 

0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/ 0/0/ 0/ 0/ 0/ 0/0/ 0/ 0/ 0/ 0/ 0/ 0/ 

/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^ 

Model="sm.mdl" ; 

/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o 

Process ; 
ELWK=3 ; 

Iiiitial={electron, positron}; 
Final =-[photon, W-plus , W-minus}; 
Kinem="2302"; 
Pend; 



2. Connect vertices witli propagators or external particles. 

There are multiple ways to connect vertices. All possible configurations are to be gener- 
ated. 

3. Particle assignment. 

Particles are assigned to propagators confirming that the connected vertex is defined in 
the model. As there will be many ways to assign particles to propagators, all possible 
configurations are to be generated. In generating the Feynman diagrams, conservation 
laws such as electric charge and fermion number conservation will be employed in order 
to avoid fruitless trials. 

The generation of Feynman diagrams borrows heavily from graph theory. In the following let 
us call a vertex or an external particle a node. Similarly, let an edge be a connection between two 
nodes, which may be a propagator or a connection between a vertex and an external particle. 
Thus an edge is expressed by a pair of two nodes (which are connected by the edge). The 
diagram or graph generation process is to construct edges in all possible ways. 

We show an example of an output file generated by the GRACE system in Fig. |H1 for the 
process we defined in the input file of Fig. namely e~^e~ W~^W~'y. The first part of this 
file describes the information about the physical process. The first generated diagram, Graph=l, 
consists of the 5 external particles and the 3 vertices. To this corresponds 7 edges that connect 
these nodes; for example, the final photon (labelled as node 2 and 3rd external particle) is 
connected to the 'yW~^W~ vertex (node 7) by edge 3. Particles are defined as incoming to the 
node. This information is used by another program such as the code for drawing diagrams. This 
reconstructs the structure of diagrams, places nodes on a graphic device and connects them by 
edges as shown in FigEI 

It is to be noted that although vertices of the same kind are not distinguished from each 
other, for instance in our example the vertex W~^W~^ appears twice in Graph 1 (see Fig. inj, 
they will be distinguished in a program, usually through a sequence of numbered labels, in our 
case node 6 and 7 (see Fig. |H|. On the other hand, since a Feynman graph is a topological 
object, it is independent of the way one assigns the sequence of numbers to nodes. Such simple 
algorithms will produce diagrams, which are topologically the same, many times. The problem 
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Figure 8: An example of information about a generated diagram in GRACE. The first part 
below defines the external particles. It is followed by defining the connections between the 
nodes and the vertices in the form "node number={edge numbers }". 

Process=l; External=5; 

0= initial electron; 
1= initial positron; 
2= final photon; 
3= final w-plus ; 
4= final w-minus ; 
Eend; elwk=3;Loop=0; 

Graph=l; Gtype=l; Sfactor=-l; Vertex=3; 
0={ 1 [positron] } ; 
1={ 2 [electron] }; 
2={ 3 [photon]}; 
3={ 4 [w-plus]}; 
4={ 5 [w-minus] } ; 

5 [order={l ,0}] ={ l[electron], 2[positron], 6[photon]}; 
6 [order={l ,0}] ={ 4[w-minus], 6[photon], 7[w-plus]}; 
7 [order={l ,0}] ={ 3[photon], 5[w-plus], 7 [w-minus] } ; 
Vend; Gend; 

Graph=2 ; 



of diagram generation is not so much to generate diagrams but to avoid such duplications. 
Although it is not so difficult to avoid generation of duplicated diagrams when limited to special 
cases such as QED JHj or tree processes, the problem of redundancy gets exasperated when loop 
diagrams are generated. 

A general method of graph generation avoiding duplicated graphs, applicable to any process 
with any order in the coupling constants, has been developed by graph theorists [HZI- With this 
method, one can calculate the symmetry factors of the Feynman diagrams at the same time. 
Unfortunately, such a method is not efficient enough for a large number of external particles 
or loops jnU- With some optimization, this method was first applied to the generation of 
Feynman diagrams in the code QGRAF 21 . Another technique of optimisation and acceleration 
was proposed in ^5^. The GRACE system exploits the latter algorithm and acceleration for 
the generation of Feynman diagrams. Appendix |j] presents in some detail the issue of graph 
duplication and calculation of the symmetry factors and the need for optimising and accelerating 
the generation of Feynman diagrams. Examples beyond one-loop and outside the electroweak 
model are given to illustrate these issues. 
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Figure 9: An example of a diagram drawn by GRACE based on the file in Fig. O Here 
each dot is a node. In this figure we have only labelled those nodes that correspond to the 
external particles according to the listing in Fig. The numbers in parentheses correspond 
to the edges as defined in Fig. also. 



Graph 1 




produced by GRACEFIG 



3 Automatic systems at tree-level 

We first start by giving tlie main lines of how the automatic code builds up and calculates 
Feynman amplitudes from the Feynman graphs at tree-level. It is very much a transcription of 
a calculation by hand. 

3.1 Squared matrix elements 

One method is to calculate squared amplitudes. Here projection operators, Dirac's gamma 
matrices and the sum over the polarizations of vector particles can be handled symbolically. This 
symbolic treatment of the mathematical expression is performed with a symbolic manipulating 
system such as REDUCE or FDRMjl2j in the case of GRACE or a specially developed package in 
the case of CompHEP Symbohc calculation in FeynCalc/FormCalc(2SI uses MATHEMATICA 
and FORM. An automatic system generates a source code for a symbolic manipulating system 
in accordance with the Feynman rules of the model. As the number of particles increases, the 
number of Feynman graphs grows very rapidly. This means that this method is not suitable in 
these cases as it requires computing the square of the number of Feynman diagrams. Moreover 
this method is clearly unsuitable if one requires information on the polarisations. 
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3.2 Spinor technique: a worked out example in GRACE 

Another method employs the spinor technique. Spinors and gamma matrices are dealt with 
in a numerical way. A library of all possible types of vertices and propagators are defined as 
subroutines to be called for a numerical calculation. The CHANEL [^ZI library for GRACE |25j and 
HELAS |7Uj library for MadGraph 3,2, are examples of such subroutines. The automatic system 
calls these libraries according to the structure of the Feynman graphs. Output of the system is 
a complete code for the numerical calculation of differential cross sections without calling any 
other package. Since this method calculates helicity amplitudes directly, it is natural to calculate 
polarized cross sections. The necessary CPU time is proportional to the number of Feynman 
graphs in this method. 

We will here show in some detail how the helicity amplitude method can be implemented 
in an automatic code such as GRACE. The method is purely numerical. The amplitude corre- 
sponding to each Feynman graph is first decomposed into vertex sub- amplitudes. Each of these 
sub-amplitudes is read from a pre-defined model file library so that one only has to call the 
library, in this case CHANEL, where these sub-amplitudes are defined. In order to achieve this, 
the propagators that appear as internal lines are expressed as a product of wave functions. 



Figure 10: A Feynman graph contributing to e^e W^W 7. 

— '^^^^ — fvvvww ^"(^-^^(^^)) 



e iP2,h2) 




W-iq2,e2i\2)) 



7(^-,e3(A3)) 



Consider the scattering amplitude corresponding to the Feynman graph shown in Fig. as an 
example, pi, P2, qi 92 and k are momenta of e"*", e~, W~^, W~ and 7, and hi and /12 are helicities 
of e+ and e~, and ei(gi, Ai), e2{q2, A2) and e3{k, A3) are polarization vectors of W~^, W~ and 7, 
respectively with the corresponding helicity Ai^2.3- 

The scattering amplitude for this contribution is given by 

Tfi = v{pi,hi)cly^ein{qi)SF{-pi + qi)c^^wu{p2,h2) 

xDFf„yiq2 + k) c^^wjil^ + ^' ~^2, -k) e2p(^2) (^3a{k), (3.1) 

where cli;^ ^wW'y express electron-H^ and photon-VK couplings, respectively, and they are 
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given by: 

eMz _ , 1-75 , . 

/2(M| - Ml,) 2 



and 



c 



= e[(p - qYg^P + {q - rf gP'^ + (r - p) Vl- (3.3) 
In CHANEL, propagators are expressed as a bi-linear form of wave functions: 

= E„.«w^°(ft':'.y')g°('.'".p'") (3,,) 

pZ — 

and 

D,Ap) = ^-"-f(^)f (3.5) 

where Wi^j and tUj are c-numbers and weight factors for the decomposition of propagator. [/" 
represents either a spinor n or f depending on the value of index a. Momenta p^*^ are calculated 
from the (off-shell) fermion momentum p. 

The amplitude then writes as a product of vertex sub-amplitudes 



7* 



D{-pi + gi, 0) D{q2 + k, mw) X! '^OL,i X! 



where 



D{p, m) 



yT/(".*) T/(".*.0 T/(0 



1 



9 9 ' 



and 



= ^(m,/ii)c^vyei,(gi)t/"((-pi + gi)»,/i»), (3.6) 
^ r(p»,/iW)c::^6«(<Z2 + fe)^(P2,/i2), 

= c^H/7(^2 + fe,-g2,-A;)e[f)(g2 + A:)e2p(g2)e3<7(A;). (3.7) 



Further details of how the spinors and polarisation vectors are represented, together with 
the weight factors can be found in |^. The fermion- fermion- vector (FFV) of the SM vertex 

parts and V^^^l are calculated with the help of the CHANEL subroutine SMFFV while the 

3-vector (VVV) vertex part V^^,^ by the subroutine SMVVV. These subroutines calculate all 
combinations of helicity states and of spinors u and v. Therefore given a set of momenta and 
helicities, these vertex parts return a number automatically. 

Summation over indices pertaining to the fermion propagators is made through a call to the 
subroutine SMCONF: 
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The last free index /, which corresponds to the propagator of a vector particle is through a call 
to the subroutine SMCONV: 



The program generates a sequence of subroutine calls in the following order: 

1) For a given process, subroutines SMEXTF and SMEXTV are called first, encoding information 
on the external fermions (SMEXTF) and external bosons (SMEXTV) in a form suitable for the 
calculation of the vertices. Since this part is common to all graphs, it is generated once 
before calculating an amplitude of the first graph. 

2) Momenta of propagators are defined by taking linear combinations of external momenta 
based on momentum conservation. Then the denominators of the propagators are calcu- 
lated by the subroutine SNPRPD. 

3) The decomposition of the numerator of the propagators as a bi-linear product of wave 
functions is carried through SMINTF for fermions and SMINTV for vectors. 

4) The different subroutines such as SMFFV or SMVVV are called for calculating vertices sub- 



5) These vertices are connected in correspondence with the propagators of the graph, which 
is realized by appropriate subroutines such as SMCONF and SMCONV. 

6) Before the amplitudes are summed over all graphs a reordering of the particles, so that 
they appear in the same sequence for all the graphs, is carried through a special routine. 

7) After summing over all diagrams, the helicity amplitudes are squared. Summation over 
spin states can also be performed automatically by a call to a dedicated routine. 

The generated code for this example is shown in Fig. 111! A more detailed description of the 
CHANEL library is given in |25j . 

3.3 Dealing with QCD and supersymmetry 

Extension of both the squared matrix elements and spinor techniques to super symmetric mod- 
els requires a proper treatment of Majorana particles. Calculation with Majorana particles 
involves the charge conjugation operator. When a symbolic manipulation system is employed, 
the symbolic treatment of this operator must be implemented. A simple algorithmjTlj has been 
constructed where this operator is confined to the conjugated vertices , so in effect we only deal 
with the usual Dirac propagator while the vertices do not explicitly invoke charge conjugation. 
In the case of the spinor technique, we only need to add new appropriate subroutines to a library, 
such as CHANEL. 

For QCD, the automatic calculation refers to the calculation of the partonic hard scattering 
part of the matrix elements. The colour factor of a graph is separated out. However, Feynman 




amplitudes. 
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Figure 11: The code for the helicity amplitude corresponding to the graph in Fig. \1(A 

He*********************** ************** ****** 

* Graph No. 25 - 1 

* Generated No . 25 
************************************************************************ 

subroutine ag25 
implicit real*8(a-h,o-z) 

include 'incll.h' 

complex* 16 atmp 

* 

* Denominators of propagators 

aprop = 1 . OdO 

call snprpd(pphase , aprop, vn21 , ajnnu(l) **2 ,0 . OdO) 
call snprpd(pphase , aprop, vn23 , amw**2 , amw*agw) 

* Internal momenta 

call smintf (amnu(l) ,pf 21 , vn21 , ex21i ,pt21i , cf 21i) 
call smintv(lepinv, amw,pf 23 , eq23b , ev23b, vn23 , igauwb) 

* Vertices (8) 

call smf fv(lextrn,lintrn,lepinv,ex21,ex21i,amlp(l) ,ajnnu(l) , 
& cwnl(l,l) ,cf21,cf21i,pt21,pt21i,eq23b,lt5,av5) 

call smf f vdintrn, lextrn, lepexv,ex21i ,ex41 ,amnu(l) , amlpd) , 
& cwln(l,l) ,cf21i,cf41,pt21i,pt41,eql7b,lt6,av6) 

call smvvvClepexa, lepinv, lepexv, -1 , -1,-1, caww,pf 9 ,pf 23 ,pf 30, eq9d, 
& eq23b,eq30b,lt7,av7) 

call smconf (lt5 , lt6 ,2,1, ex21i , av5 , av6 , ItS , av8) 

call smconv (lt7 , lt8 , 2 , 2 , ev23b , av7 , av8 , It , av) 

sym = - l.OdO 

cf = + l.OdO 

aprop = cf*sym/aprop 

indexg(l) = 3 

indexg(5) = 4 

if (j cpol (3) .ne . 0) call smcpoKl, It, av) 
call amporddt, av, indexg, agcwrk) 
ancpCj graph) = O.OdO 
do 500 ih = , lag-1 

atmp = agcwrk ( ih) *aprop 

agc(ih,0) = agc(ih,0) + atmp 

ancp(jgraph) = ancp(jgraph) + atmp*conjg(atmp) 
500 continue 
return 
end 
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rules for the four point vertex of the gluon is expressed in a mixed form of colour matrices, 
through the 5*^7(3) structure constants, f°''"^ and Lorentz parts: 



9's 



^ jace jbde jade ycbe^ g ^ _j_ ^ yafee jcde jade jbce 

abe fcde 



_j_ ^ jace jbde jabe j 



There are two ways for the separation of these factors. One is that the Feynman graph is 
generated with the original quartic vertex but the amplitude generator expands this vertex into 
a sum of three terms where the colour index is factored out. The other method, also used in 
CompHEP for example, introduces an auxiliary field which interacts through a vertex with two 
gluons. The Feynman graph generator constructs the quartic gluon vertex as the sum of graphs 
consisting of s, t and u type exchanges of the auxiliary field. The auxiliary field method is also 
best suited in the case of the MSSM where additional coloured particles with 4-point vertices 
exhibiting a variety of colour structures are needed. 



3.4 Checking the result at tree-level 

The SM and the MSSM are gauge theories. One can exploit the gauge freedom to check the 
result of the calculation of a cross section or helicity amplitude. For the calculation based on 
the spinor technique, it is easy to keep the gauge parameters as variables in the generated code 
without increasing the CPU time. For the squared matrix technique that relies on symbolic 
manipulation, checks that are based on varying the usual 't Hooft-Feynman parameter lead to 
very complicated and very large expressions. Checks are much more manageable with generalised 
non- linear gauges that will be discussed in detail later in Section |21 



3.5 Automatic tree calculations without Feynman graphs 

Let us briefly mention that there exists codes for the automatic calculation of amplitudes that 
do not require the standard textbook approach based on Feynman graphs. ALPHAjHSl uses an 
algorithm based on an iterative solution of the equation of motion whereas HELAC [SHI is based 
on solving the Schwinger-Dyson equation. An acceleration method is developed in [SZj and is 
applicable also to the methods using Feynman amplitudes. All these methods exploits a binary 
representation of the momenta involved in the tree-level amplitude calculation. 

Let Pi, « = 1, A'^ be the momenta of the external particles for a certain amplitude. Mo- 
mentum conservation imposes the restriction X^jPi =0. A momentum g of a propagator in a 
tree graph is expressed as g = J2i'^iPi^ where is either or 1. Thus the possible number of 
momenta appearing in the calculation is bounded by 2^~^ — 1 and only a finite limited number 
of vertices appear in the calculation. This helps to construct tree A^-point Green's function 
combining these vertices in a kind of matrix operation. 



3.6 Kinematics and event generation 

Having obtained the amplitudes or the squared matrix elements, the last step is the integration 
over phase space to obtain the cross section. One may also want to generate events for the 
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analysis of experimental data. For multi-particle final states, phase space integration means 
integrating over many variables. Since the amplitudes themselves, the integrand, are most often 
lengthy and have a complicated structure in these variables, Monte Carlo integration packages 
such as VEGAS f72] and BASES ^ are used. Even if the integrand were a smooth function over 
the whole phase space without narrow peaks, naive use of a numerical integration package may 
cause problems as it would be difficult to reach the required accuracy especially for a large multi- 
dimensional space or when one deals with complicated boundary configurations. Narrow peaks 
can occur for example where the momentum of an internal particle is such that it approaches 
its mass-shell. The amplitude could blow up. As these peaks give an important contribution to 
the results, it is necessary to catch the detailed structure of these singularities and to integrate 
around these peaks within good accuracy. In this case the Monte Carlo integration package 
needs to accumulate enough sampling point around these peaks. When singularities run parallel 
to integration axes, VEGAS or BASES optimizes the distribution of sampling points looking at the 
projected values on the axes. 

It is fair to say that it is this stage of the automatic calculation of cross sections in high 
energy physics which requires most intervention by the user. This is due to the fact that there is 
at the moment no general purpose integrator that can automatically catch the different types of 
singularities and suitably adapt itself to properly treat all the singularities present in a process 
with enough accuracy. Several integration packages have been developed [7311^ ESI handling 
many singular situations but none is general enough that it can be trusted to run independently 
and give an accurate result in any situation. The user should select good integration variables 
so that the singularities are arranged along the integration variables. This selection of the 
integration variables depends, of course, on the process. Some systems, such as GRACE , include 
a library of kinematics suitable to handle a number of topologies for the singularities for a given 
number of external legs. 

In the case of QCD, the partonic cross sections and amplitudes need, for example, to be 
convoluted with parton distribution functions. The code can also provide packages for parton 
shower and hadronisation. 

Once a numerical code of a differential cross section is obtained, one can generate unweighted 
events by a program package such as SPRING ^7] using the result of the integration calculated by 
BASES. Unweighted simulated events, generated in accordance with the theoretical predictions are 
very useful tools to analyse the experimental data and take into account the complex structure 
of the detectors. 

4 Automatic systems for one-loop processes 

Extension of an automatic system from tree-level to one-loop processes is, in principle, straight- 
forward because all the ingredients needed for one-loop calculations in a renormalisable theory 
are known. In practice, however, quite a few of the new features not met when handling pro- 
cesses at tree- level pose a real challenge and are fraught with technical difficulties. This explains 
why although there has been a proliferation of automatic systems for tree processes since the 
beginning of the 1990s, only a couple can tackle one-loop processes. Currently two systems are 
available, GRACE-loop and FeynArt/FeynCalc/FormCalc[23 EHl EOl EI] which apply both to 
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the electroweak sector of the SAi and the AiSSAi. Apphcation to QCD at one-loop is at the 
stage of development by some groups but a general purpose code is still not completed yet. In 
this section we describe what kind of problems one encounters in the construction of an auto- 
matic system at one-loop based on the experience we gained while developing GRACE-loop for 
electroweak processes in e~^e~ colliders. 

An automatic system for one-loop should be composed of the following ingredients: 

1. Proper graph generator. This should provide all one- loop diagrams for a given process. 
Most graph generators that work for tree-level handle at least one-loop diagrams. We 
have already described the general algorithm of graph generation in Section 12.31 

2. Amplitude generator. This very much depends on how to calculate one-loop amplitudes. 
Both the helicity amplitude formulation and/or the squared matrix elements technique 
are possible as explained later. 

3. The model file should include a library for counterterms. This should contain all terms 
necessary to compensates ultraviolet divergences. For a given theoretical model, one can 
have different libraries reflecting the fact that one can choose between different renormal- 
isation and regularisation schemes. In GRACE, ultraviolet (UV) divergences are regularised 
through dimensional regularisation and the electroweak sector is renormalised on-shell, 
see section 16.11 Note in passing that there should also be a prescription for handling in- 
frared(IR) divergences. In GRACE this is done by giving the photon an infinitesimal mass, 
see section lR.1.21 These libraries have a structure similar to those defined for the tree-level 
model files. 

4. A library for one-loop integrals. This should include analytic formulas for 2-, 3- and 4- 
point scalar integrals, see section 17.11 5- and 6-point functions can be expressed by a 
sum of 4-point functions, section 17.61 As we will see in section 17.11 a scalar integral is an 
integral over the loop momentum where the integrand is a product of the denominators 
of the various propagators inside the loop, so that no loop momentum appears in the 
numerator of the integrand. When products of loop momenta are involved one speaks of 
higher rank tensors. 

5. Beside the aforementioned reduction of 5- or 6-point diagrams to 4-points, one needs 
another library which decomposes higher rank tensors of a box diagram to relevant scalar 
integrals and surface integrals, namely 3-point functions. The reduction formulas for 
higher rank tensors of vertex diagrams is needed as well. There are a few algorithms dealing 
with the last two points. One needs very efficient algorithms and numerical routines here, 
since this steps can take up a large fraction of the CPU time. This is developed further 
in this article in section [731 

6. Ways to check the results. This is essential and must be performed. We consider it as 
the most important part of the system because the calculation is much more complicated 
than the tree processes. Lacking these tests one cannot be convinced that the obtained 
results are correct. Usually, one checks UV and IR finiteness. Gauge invariance is also 
a very powerful tool, however the familiar 't Hooft gauges are not suitable for multi-leg 
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one-loop processes. GRACE incorporates a series of parameters derived from a generalised 
non- linear gauge. 

The theoretical background and the algorithms necessary for the construction of the libraries 
related to the points 3,4,5 above will be developed in detail in sections l8. 1 . ![ IKTT^ 17. 1(17.61 and 
|721 respectively. The gauge-parameter independence check in the last point above is crucial and 
could be considered as pivotal. Beside being a most powerful check on the calculation it also 
has a bearing on how the libraries for the reduction of the one-loop integrals are constructed. 
The choice of special gauges here is almost mandatory. This is the reason we pay so much 
attention to the implementation of the non-linear gauge in this article. Although we will, in 
the next sessions, get back to the details of all the points listed above and before we point to 
some features related to the technical implementation of the one-loop radiative corrections in 
the electroweak part of the SAi in the automatic code GRACE, it is worth to briefly expose the 
gauge-parameter check. 



4.1 Gauge parameter independence 

The most stringent test for the consistency of the calculation will be provided by the gauge- 
parameter independence of the cross section. In tree-level calculations with an automatic code 
any gauge fixing will do. This is because the inclusion of any modifications brought by the gauge 
fixing Lagrangian and associated ghost Lagrangian is always easy since no loop integration is 
needed. Typically, the standard gauge fixing is used where the gauge parameter hides in the 
propagators of the gauge particles 7, Z, W^, and unphysical scalars of the theory. For example 
in a general i?^ gauge, the Z propagator writes 



(4.8) 



This formulation of gauge fixing, however, is not well suited when loop calculations are in- 
volved. The "longitudinal" part, in the second term of the propagator, constitutes an obstacle 
in carrying out the loop integration, due to the appearance of not only unphysical thresholds, 
^z^'zi functions but also the introduction of higher rank tensors, q^(f ■, and higher order Ap- 
point functions since this part involves "two" denominators. Since as we will see, the reduction 
of tensor loop integrals and higher A^-point functions is very much time consuming it is desir- 
able to find another formalism for gauge fixing which allows the propagators to take as simple 
form as possible, = 1- This can in fact be realized in the non-linear-gauge(NLG) fixing 
formulation f^l'l. With this kind of gauge fixing, the propagators can be taken to correspond 
to = 1 for all gauge particles but additional gauge parameters, through derivatives in the 
gauge-fixing, modify a large number of vertices where the gauge and unphysical (Goldstones 
and ghosts) are involved. We will take a generalised gauge with five gauge parameters in total, 
a, /?, (5, e, K. Although this means some more work in redefining the model file as compared to 
the linear gauge which only modifies a few propagators and also because renormalisation is 
slightly more involved than in the by-now textbook linear R^ = \ gauge, implementation of the 
non-linear gauge is worth the investment if one is interested in a general purpose automatic loop 
calculator. One then has the possibility of performing gauge-parameter independence checks, 
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involving the 5 parameters, on a large scale calculation. We will see in section 15^ that, though 
each diagram out of a set of some hundreds diagrams shows a dependence in a gauge parameter, 
summing over all diagrams gives a gauge independent result. This is a very powerful check on 
the calculation and on the automatic system. 

4.2 Generation of matrix elements 

As is in the case of tree processes, in general the matrix elements can be obtained either through 
the helicity formalism or one can obtain directly the squared matrix elements with the help of 
some symbolic manipulation programs. GRACE-loop has both versions. 

The procedure to get the helicity amplitudes is as follows. Once a diagram is generated a 
REDUCE code is produced which contains an expression for the corresponding amplitude. This is 
just a "reading of the amplitude" in REDUCE. We do not take any trace of 7 matrices unless a 
fermion loop is formed inside the diagram. The role of REDUCE is only to rearrange the generated 
terms. First, the Feynman parameters for the loop integration are introduced, see Sectiondand 
Eg. 17.21 The loop momentum I, see Eg. 17. II is shifted so that one has denominators containing / 
only through P. The amplitude now contains not only the loop, external momenta and Feynman 
parametric variables but also strings of 7 matrices. 

Operations are done in n-dimension. All the Lorentz contractions are taken in n-dimension if the 
pair of indices does not bridge two fermion lines. Then, without taking the trace, every product 
of 7 matrices along the fermion lines is replaced by a corresponding symbol which works as a 
function when the code is converted to a FORTRAN source code. The remaining contractions and 
the calculation of the products of 7 matrices are left to be done numerically in particular by the 
CHANEL library that we already described in the case of tree- level amplitudes. As stated, REDUCE 
is used, in this approach, to get the final form of the amplitude, but only for the rearrangement 
of terms and the functions such as spur (Trace) and index (contraction) are applied only to 
quite a limited part of the manipulation. 

Now the obtained amplitude can be regarded as a polynomial with respect to the loop integration 
parameters. Each monomial is replaced by a symbol which represents a one-loop integral with 
the numerator corresponding to that monomial. Thus the amplitude can be calculated when a 
library of loop integrals is supplied, see Section [7| Since we do substantially nothing with the 
amplitudes, it is quite easy to get them. 

On the other hand, in the second approach traces of all 7 matrices and contractions of indices 
are applied systematically on the product of a tree amplitude with a loop amplitude. Again 
dimensional regularisation is assumed in the case of UV divergent diagrams. The coefficient 
of a monomial of the parameters for the loop integration is merely a polynomial composed of 
external momenta. That monomial is also replaced by a symbol which is a substitute for the 
corresponding loop integral. This approach generates much larger expressions than the helicity 
approach, see Tabled It may happen that REDUCE is not able to complete all the manipulations 
because too huge intermediate expressions are produced. 

The helicity amplitude formalism is certainly desirable because it can provide the spin infor- 
mation. However, the sums on spin states required for external and internal particles produce a 
certain number of arithmetic operations. This would not cause any problem for 2^2 processes 
where the total cross section is obtained by integrating the matrix element with respect to only 
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one variable, the scattering angle. When the final state contains more than 2-body, however, one 
needs to perform a multi-dimensional integration over phase space with, for example, 4 variables 
for 2^3 and 7 for 2^4 processes. If one relies on integration packages that are based on 
Monte-Carlo algorithms, one has to compute the loop matrix element many times, typically 
in the order of 50K — IM times in total. To get the 0{a) corrections the contribution of the 
one-loop diagrams must be combined with the cross section for real photon emission. Usually 
the latter almost compensate the virtual corrections. Thus the number of points also should be 
properly chosen to achieve good accuracy. Hence, at least at present, it is not always realistic 
especially for complicated processes such as 2 ^ 3 and 2 — > 4 processes to use the method based 
on the helicity amplitudes for the computation of the total cross section as this will require a 
lot of CPU time. The helicity amplitude approach is, however, helpful for comparing and check- 
ing, at several points in phase space, the squared matrix elements obtained the symbolic way, 
through REDUCE for example. Another limitation of the symbolic way, however, is that it cannot 
be applied to a process that has no tree-level contribution such as 77 77 and 77 — > 7Z. For 
such cases one has to use the helicity amplitude approach. 

4.3 Regularisation scheme and the issue of 75 

An important step in the proof of renormalisability of the SAi is that the symmetries that 
are present at tree-level are still preserved at the loop level. With a theory like the SA4 that 
involves both vector and axial currents, it has been known that all gauge symmetries may not 
be preserved at the loop level due to the Adler-Bell-Jackiw anomaly j76j. In the case of the 
SM, the charge assignments of all particles within a family is such that the anomaly is not 
present [77^. It is important to stress that the existence of a genuine anomaly has nothing to do 
with regularisation in the sense that there could be no regularisation that can make the anomaly 
disappear while maintaining all other quantum symmetries. Nonetheless an unfortunate choice 
of regularisation can induce an apparent anomaly that would violate some Ward identities that 
could stand in the way of a proper renormalisation programme. The use, even in QED, of a 
naive cut-off is such an example. These apparent anomalies can be removed by the introduction 
of extra counterterms usually not obtained through multiplicative renormalisation. This, of 
course, unduly complicates the renormalisation procedure and can make the implementation 
at the automation level more problematic. Pauli-Villars regularisation 75^ that works so well in 
QED fails in the case of non-Abelian gauge theories because a mass term is introduced in a naive 
way. Dimensional Regularisation PSI I39j is a very powerful and extremely practical scheme for 
regularising divergences in non-Abelian theories^. In DR, the loop integrals as well as the Dirac 
and tensor algebra are calculated in a space-time with arbitrary n dimensions. It is relatively 
easy to implement in an automatic code for one-loop calculations. The divergences are isolated 
as poles in n — 4 through, for example, the variable Cuv as is done in GRACE, with 



^For a step-by-step presentation of the method, see for example [ZH] or [HOj. The chiral anomaly is 
also nicely exposed in these textbooks. 



Cuv = 7£; log 47r, n = 4 - 2e. 
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For the electroweak SAi , the problem with DR is the treatment of 75 which is a 4-dimensional 
object and hence can not be continued to an arbitrary n-dimension. With calculations involving 
75 one therefore needs some additional scheme or prescription. In the original 't Hooft-Veltman 
prescription j38| l8T] . 75 anticommutes with the 7 matrices that carry one of the n = 4 dimensions 
but commutes with the rest. This split is not satisfactory as it breaks Lorentz invariance in 
the full n-dimensional space and does not manifestly respect the conservation of the gauge 
current [H21 IHUl- Many variations on the original scheme, that may be considered as not being 
fully mathematically consistent, have been applied in different calculations, see|83j for a review. 
For applications to the anomaly- free SM, the most efficient and practical scheme, especially 
from the point of view of implementation in a computer code, is the so-called Naive Dimensional 
Regularisation whereby 75 is taken as fully anticommuting with all the n-dimensional 7 matrices. 
In GRACE-loop this is what is implemented. Let us also point out that the cyclicity of the 
trace is not used. For loop calculations in supersymmetry DR is not quite appropriate since 
it breaks supersymmetry by splitting the number of fermionic and bosonic degrees of freedom. 
An alternative is Dimensional E,eduction|84| l85] . A very nice review is given in |86[ l87]. At the 
one- loop level an equivalent prescription is Constrained Differential Renormalisation'88' , CDR. 
Reference j3Ul I31j shows how both DR and CDR can be implemented in the automatic code, 
FormCalc for the calculation of one- loop diagrams. 

4.4 Size of the generated programs 

One of the difficulties with running codes of one-loop processes produced by an automatic system 
has to do with the huge increase in the number of Feynman diagrams, as compared to tree- level 
processes, combined with the equally large size of the generated computer file for the matrix 
elements of each diagram. We could then easily end up with a total size of a program for a given 
process which is so huge that it is not always easy to compute the cross section in a realistic time 
scale. To give a feel for the size of such programs, we show in Tabled examples for some processes 
that have been calculated so far with up to 6- leg final states (2 — > 4 processes). Here the total 
number of diagrams in a general non-linear gauge is shown. Once the gauge parameter check is 
successfully passed, for the cases of 2 ^ 3 and 2 — > 4 processes the amplitudes are generated 
again in the linear gauge, which corresponds to a particular choice of the non-linear gauge as 
will be explained and by switching off the electron Yukawa coupling. This determines the set 
of amplitudes for the "production job" to be supplied to a phase space integrator. We note in 
passing that switching off the electron Yukawa coupling does not mean that all occurrences of 
the electron mass are set to zero. In particular the electron mass is kept in order to regularise 
collinear divergences, see section 15^ or for the renormalisation of the fine structure constant in 
the charge counterterm, see Appendix d 

The size of the generated codes refers only to the amplitude and does not take into account 
common necessary libraries. In the case of 4-fermion production, to save CPU time the masses 
of the light external fermions are neglected, as long as they appear in the numerator of the 
matrix element. At some points in phase space the full set and the production set are compared 
to confirm that the latter is correctly generated. The actual CPU time for e^e~ — > e^e^H, for 
example, amounts to 150K hours equivalent to a single IBM RS/6000(375MHz, Power 3) CPU 
in quadruple precision. 
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Table 1: Size of the generated code for some cross sections measured by the number of 
lines. We also show the total number of tree-level and one-loop diagrams for each process in 
some general gauge. The numbers in parenthesis have the small electron Yukawa couplings 
switched off. "helicity" refers to the method of generating the helicity amplitude. REDUCE 
refers to the direct computation of the full matrix elements squared. 
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5 The Standard Model in a general non-linear gauge 
5.1 The classical Lagrangian 

To help define our conventions and notations we first introduce the classical Lagrangian of the 
SM which is fully gauge invariant under SU{2) x U{1). 

We denote the gauge fields of the theory of the SU{2) x U{\) group as W'i^{a = 0,1,2,3). 
The weak isospin triplet refers to 1,2,3 and the hypercharge singlet to the component. The 
corresponding gauge couplings are g"". The gauge invariant field strength F^^ writes in a compact 
form as 



F^,-d,W,-d,W^+g f W^W, ^ 9 - [g' ^a = 0) (otherwise) " 

(5.1) 

This leads to the pure gauge contribution Cg 

Cg = -\f-,.F'^^^ . (5.2) 
The gauge interaction of the matter fields is completely specified by their isospin and hypercharge 
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(y) quantum numbers, such that the electromagnetic charge is Q = T^ + Y, and the covariant 
derivative, 



3 
a=l 

with = cr"/2, where a" are the usual Pauli matrices. Left-handed fermions L of each gener- 
ation belong to a doublet while right-handed fermions R are in a SU (2) singlet. The fermionic 
gauge Lagrangian is just 

CF = iJ2L-f''D^L + iJ2R7^D^R, (5.4) 
where the sum is assumed over all doublets and singlets of the three generations. 

Mass terms for both the gauge bosons and fermions are generated in a gauge invariant way 
through the Higgs mechanism. To that effect one introduces a scalar doublet with hypercharge 
y = 1/2 that spontaneously breaks the symmetry of the vacuum through a non-zero vacuum 
expectation value v 



* = ( „ ;V43 ) = ((»+ H -L)/v/2 ) • (°i*> = ( ./v-2 ) ■ p-^' 

The scalar Lagrangian writes 



£5 = P^<^)^W)+>Cpot , >Cpot = -n</') = /^V0-A(<^W = -A(</'V-y (5.6) 



The Nambu-Goldstonc bosons X^;X3 hi Cs get absorbed by the Z and to give the latter 
masses {Mz,w), while the photon A remains massless. The physical fields W^,A,Z relate to 
the original W quartet as 



with 



A.. 



cw 



V2 ' 



^ swW^ + cwW^ 

cwWl - swW?, 



(5.7) 



12 ■ 



sw 



(5.8) 
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the electromagnetic coupling e 



e , e 
9 = , 9 = — , (5.9) 



and the masses 

Mw = ^, Mz = -^^. (5.10) 
Cs also defines the mass of the Higgs 

Mfj = 2^?. (5.11) 

//, A and v are not all independent parameters, v is defined to be the minimum of the scalar 
potential. This is equivalent to requiring no tadpole in Cs- In other words we require the 
coefficient, T, 

r = t>(//2 - At;2), (5.12) 

of the term linear in H in Eq. 15.61 to be zero, T = 0. We will impose this requirement to all 
orders. 

Fermion masses require the introduction of a corresponding Yukawa coupling, fu {fo) for an 
up-type fermion / (for a down-type fermion) 

Cm = - ^ fuLu4>Ru - X! fDLD(t>RD + {h.c), cj) = ia'^cj)* , mu,D = "^^'^^ ■ (5-13) 

up down ^ 

Instead of the original set of independent parameters {9, 9' , X, n'^ , fu,D}, it is much more 
advantageous to revert to an equivalent set of physical parameters that is directly related to 
physical observables, namely {e,Mw,Mz,MH,mij,m£)}. In this respect note that the weak 
mixing angle is just a book-keeping quantity that will be defined, at all orders of perturbation, 
in terms of the masses of the vector bosons: 

CW = jj^. (5.14) 

If one allows v to be an independent parameter we will trade it for the tadpole, T, which we 
will add to the list of independent parameters that specify the theory. 



5.2 Quantisation: Gauge-fixing and Ghost Lagrangian 

As known because of the gauge freedom in the classical Lagrangian Cc, Cc = Cq+Cf+Cs+Cm, 
a Lorentz invariant quantisation requires a gauge fixing. We generalise the usual 't Hooft linear 
gauge condition to a more general non-linear gauge that involves five extra parameters j54|. 



G+ 

Cgf = -^\{^^^ - ieaA^ - igcwPZ^)W^'+ + ^{y^{v + 6H + ifix3)x'' f 
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= _J_G+G--^(G^)2-^(G^)2 (5.15) 
iw 

Note that it is not essential for the Feynman parameters i'y^ z that appear within the func- 
tions G^'^ to be equal to those that appear as factors of G^G~ {^w) and G\{^z)- However in 
this case Cgf does not cancel, at tree-level, the mixing terms x-^-, X^-^-, see for instance |5nil5n] . 
To avoid this unnecessary complications we will stick to z — ^w,z- 

To construct the ghost Lagrangian we will require that the full effective Lagrangian, or 
rather the full action, be invariant under the BRST transformation (the measure being invariant). 
The required set of transformations needed to construct the ghost Lagrangian together with the 
definition of the ghost fields can be found in Appendix This is a much more appropriate 
procedure than the usual Fadeev-Popov approach especially when dealing with the quantum 
symmetries of the generalised non-linear gauges we are studying. This implies that the full 
quantum Lagrangian 



CQ = Cc + CGF + CGh, (5.16) 
be such that (Jbrs-Cq = and therefore 5brsCgf = -^brs^g/i- 

Moreover we appeal to the auxiliary 5-field formulation of the gauge- fixing Lagrangian Cgf- We 
will see later that this formulation is also very useful to extract some Ward-Takahashi identities. 
Within this approach 

Cgf = iwB+B- + yl^^P + yl^^P + B-G+ + B+G' + B^G^ + B^G^. (5.17) 

From the equations of motion for the i?-fields one recovers the usual Cgf together with the 
condition B^ = - J (^ = {iw,iz.U])- 

Defining the anti-ghost, c*, from the gauge fixing functions, we write 

6bksc' = B\ (5.18) 

Then by identification 

Cgh = - (c+<5brsG+ + c-<5brsG- + c^^brsG^ + c^<5brsG^) + ^Rs-Cch 

= Cfp + SBRsJ^Gh ■ (5.19) 

That is, one recovers the Fadeev-Popov prescription, Cpp, but only up to an overall function, 
^BRS^Ghi which is BRST invariant. The complete Feynman rules we list assume Cgh = which 
is sufficient for one-loop calculations. For higher orders a counterterm not of the Fadeev-Popov 
type, but which is BRST invariant on its own, may be required to renormalise a quartic ghost 
vertex, in this case one can take Cgh = ^^ijk^BRsic^c^ c'^) ■ The full set of Feynman rules derived 
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from Lq is relegated to Appendix El These Feynman rules are derived with an arbitrary set 
C = (d, /3, 5, K, e), = (CvF) '^Z) '^a) although for one-loop applications and for all our tests we 
stick with the 't Hooft-Feynman gauge = = = 1 where the gauge boson propagators 
take a very simple form. Only their "transverse" part g^y contributes, see Eq. 11.11 This also 
greatly simplifies the calculations not only because the expressions get more compact but also 
because of the fact that the longitudinal parts introduce a high degree of (superficial) ultraviolet 
divergences. Although these greatly simplify when adding, all diagrams the cancellations are 
very subtle and may be not efficiently handled when implemented numerically. In practical 
calculations one can also tune a,/?, .. so that one reduces the number of diagrams and simplify 
some of the vertices. For instance for photonic vertices a = 1 is to be preferred since there is no 
W^X^ A vertex and also because the WWj simplifies considerably. One can also choose (3 so 
that W^x^Z vanishes. 

6 Renormalisation and counterterms 
6.1 Renormalisation constants 

The renormalisation procedure follows very closely the on-shell renormalisation scheme, carried 
injQOj in the case of the usual linear gauge. The set of physical input parameters includes all the 
masses of the model together with the value of the electromagnetic coupling as defined in the 
Thomson limit. As explained above we also add the tadpole, T, to this list. Renormalisation of 
these parameters would then lead to finite S-matrix elements. For the mass eigenstates and thus 
a proper identification of the physical particles that appear as external legs in our processes, field 
renormalisation is needed. S-matrix elements obtained from these rescaled Green's functions will 
lead to external legs with unit residue. Therefore one also needs wave function renormalisation of 
the fields. In the linear gauge with all ^ = 1 this also renders Green's functions finite. Especially 
for the unphysical sector of the theory, the precise choice of the fields redefinition is not essential 
if one is only interested in S-matrix elements of physical processes. We will therefore concentrate 
essentially on the renormalisation of the physical parameters and physical fields, although we 
also introduce field renormalisation for the Goldstone bosons. 

All fields and parameters introduced so far in section 3 are considered as bare parameters 
with the exception of the gauge fixing Lagrangian which we choose to write in terms of renor- 
malised fields. Care should then be exercised when we split the tree-level contributions and the 
counterterms. In Appendix El we also present the alternative approach where the gauge-fixing 
term is also written in terms of bare parameters. Differences between the two approaches, of 
course, only affect the unphysical scalar sector. In Appendix O we derive some useful Ward 
identities that constrain the two-point functions in this sector. 

For the renormalised quantity X, the corresponding bare value will be defined by an under- 
lined symbol, X, and its counterterm by 5X . 

For the physical parameters, and the tadpole, we define 

Kir = M^ + 5M^, 
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A/f2 _|_ j: ]\/r2 

— IVl^ -\- OlVl^, 




= rnf + 6mf, 


Ml 


= Ml + 5Ml, 


e 


= Ye = {l + 5Y)e 


T 


= T + 5T. 



We now turn to the fields and the wave function renormalisation constants. 
1. Gauge fields 

= , = 1 + Szli\ 

A, 

1/2 ./^.' „. _X7l/2 




(6.1) 



(6.2) 



2. Fermions 



For simplicity we will assume no fermion mixing and therefore no CV violation. In this 
case the i 
see later. 



1/2 

case the wave functions renormalisation constants ^Zjj^ jrj^ can be taken real as we will 



-L,R ^'^JL.R J L,R ■> ^^Il,R ^^"^fL,R 

3. Scalars 

S = VZsS , VZs = l + 6Z'J^ , S = H,x^,X3. (6-4) 

Because we are only presenting an application to processes at one-loop, there is no need to be 
specific about the renormalisation of the ghost sector. This is sketched in Appendix lO Suffice to 
say that to generate the ghost Lagrangian including counterterms one needs to re-express Cqf 
in terms of bare fields to first generate, through BRST transformations, Cqh with bare fields. 
This is because Cgf is written in terms of renormalised fields and as such does not induce any 
counterterm. However the BRST transformations are defined for bare fields. 

The generated counterterm library for all 3 and 4-point vertices is listed in Appendix^l The 
counterterms are fixed through renormalisation conditions that are set, with the exception of 
the e^e~ A vertex^, from the two-point functions to which we now turn. 



6.2 Two-point functions at one-loop including counterterms 

We work in the on-shell scheme closely following [^0] for the determination of the renormalisation 
constants. The renormalisation conditions on the parameters are essentially derived by properly 
defining the masses of all the physical particles A, Z,W^, f and the electromagnetic constant. 

^In fact, as we will see, due to a Ward identity, the counterterm for the charge can be expressed in 
terms of two-point functions also, see Eq. 16.231 
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They are all set from the propagator and the e'^e~A vertex. Let us first turn to the propagators 
of the fields of the theory. 

The counterterm contribution will be denoted by a caret while the full contribution (coun- 
terterm and one-loop diagrams contribution) is denoted by a tilde, so that for the vector boson 
we may write 

n = n-MI. (6.5) 

Moreover it is necessary to decompose these contributions according to their Lorentz structure. 
For our purpose we will only consider the case of no mixing (and hence no CV violation) in the 
fermionic sector. The decomposition of two-point functions is as follows. 



type 


formula 




vector-vector 




= {9,. - '-f) mq') + ^n,(g^) 


scalar-scalar 


n(g2) 




vector-scalar 




I {q is the momentum of the incoming scalar ) 


fermion-fermion 


S(g2) = 


Kil + + + ^57^75 



Complete one-loop results for all two-point functions in the generalised non-linear gauge are 
collected in Appendix IHl 

The contribution of the counterterms to the two-point functions writes 
1. Vector- Vector 



WW 




= 6M^ + 2{M^ - q^)6Zli' 




fiY- 


= 6M^ + 2M^6Zli^ 


zz 




= 5Ml + 2{Ml-q^)5zYz 




nf^ 


= 5 Ml + 2Ml5Z^Jz 


ZA 


nf^ 
nf^ 


= {Ml - q^)6ZYi - q'SZY^ 

= MlSZ'Jl 


AA 


nr 


= 



2. Scalar-Scalar 



HH 




= 2(,2 - M%)5zf -5Ml + ^-^ 


X3X3 


nx3 


= 2qHz]l^ + ^ 


XX 




= 2qHz]!^ + ^ 
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3. Vector-Scalar 



Wx 




= Mw{6Mw/Mw + 5Z]i^ + SZ]/"^) 


ZX'i 




= Mz{6Mz/Mz + SZ]/^ + ) 


^X3 


n^x3 


= MzdZ'J^ 



4. Fermion-Fermion 

At one- loop, this sector is unaffected by the parameters of the non-linear gauge and thus 
all functions are as in the usual linear gauge case. As mentioned earlier, all wave functions 
constants are real since we do not consider mixing in the fermionic sector. 

Ki = —nif [bZ^j^ + 6Zyj^^ — 6m f, 
K5 = 0, 

K, = {6Z}i' + 6Zy^), 

k,, = =- (5Z]'y - 5Z]'^) . (6.6) 
6.3 Renormalisation Conditions 

Leaving aside the renormalisation of the electromagnetic charge, these two-point functions give 
all other renormalisation constants. Before deriving these let us first turn to the tadpole. 

1. Tadpole 

The counterterm for the tadpole contribution, 6T, is defined such that the tadpole loop 
contribution T^°°p and the counterterm ST combine such that the tadpole T = T^°°p + 
(5T=0. Then 



5T 



rploop 



(6.7) 



The tadpole contribution in the electroweak SM is sometimes necessary. An example is 
the loop two-point functions of the massive vector bosons and the Higgs in order to check 
the BRSTjHSllEE] or the Slavnov- Taylor 191] identities, see for example Appendix IdI 

2. Charged vector 

The conditions specify that the pole-position of the propagator is M^, and that the 



residue of the propagator at the pole is 1. 



0, 



M^(M4, 

This gives the following relations. 

5M^ = -^eU^iM^r), 



d 



. 



5Z, 



1/2 

w 



(6.8) 



(6.9) 
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3. Neutral vector 



The conditions to be imposed on the photon-photon and Z — Z self-energies are the same 
as for the W-W transition. In addition we require that there should be no mixing between 
Z and the photon at the poles = 0, M|. 



d 



(M|) = 0, ^3?en|^(g2^ 



n^^(o) = o, 



= 0, 



0, 



(6.10) 



(6.11) 



n|^(0) = 0, Mf^(M|) = . (6.12) 

Among these 6 conditions, ^^^^(O) = produces nothing, except that it ensures that the 
loop calculation does indeed give llj^^{0) = 0. One then derives. 



5M| = (Ml), SZ'J^ = ^5ie^n|^(g2) 



g2=M| 



,1/2 ^ 



5z'Jt = -^n^^(o) , 



5Zyi = -nr(0)/M|, 5Z)l^ = (M|)/M| . 



tZA 



1/2 



tZA, 



(6.13) 

(6.14) 
(6.15) 



4. Higgs 

The conditions specify that the pole-position of the propagator is , and that the residue 



of the propagator at the pole is 1, 



d 



This gives the following relations. 

3(5T 



1) Zi Ci(^ 



(6.16) 



(6.17) 



5. Fermion 

The conditions for pole-positions and residues arc the same as for the other physical 
particles. Also the vanishing of 75 and 7'^75 terms at the pole is required. These conditions 
read 



d 



mf^eKj{mj) + ^eKi{mj) = 0, — 3fJe UK^{q^) + Ki{q') 

^ek^{m)) = 0, ^ek^^{m)) = 0. 



= 0, 



(6.18) 
(6.19) 
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CV invariance leads to = 0. In this case, one can take both 5Z^j^ and SZ^J^ to be real 
using the invariance under a phase rotation. We obtain the following relations. 



f 

(6.20) 

6. Charge 

While there are many vertices in the theory, if the charge e is properly renormalised, we 
do not need any further renormalisation conditions. The condition can be imposed on any 
vertex. The most natural one is to fix the e^e~ A vertex as is usually done in QED by 
relating it to the Thomson limit. The condition requests that the coupling is — e when q, 
the momentum of the photon, is 0, while the with momentum p± are one shell, 

{e^e^ A one loop term + e^e~ A counter term)! „ o = . (6-21) 
The counterterm is defined in Appendix IFl From this, we obtain 5Y . In fact we will see 

ity\-ro 1 /2 

that due to a Ward identity, see for example |^, 5Y writes as a combination of 5Zjl^ 
and SZ^J"^ which is valid in all gauges. 

7. The unphysical sector 

Because we are interested in applications to physical processes, the renormalisation of 
this sector is not adamant. Nonetheless one may choose to work, as far as possible, with 
finite Green's functions involving the Goldstones and the longitudinal modes of the vector 
bosons. With a linear gauge-fixing condition in the 't Hooft-Feynman gauge, and in the 
approach we are taking where the gauge-fixing Lagrangian is written in terms of renor- 
malised fields from the outset, all divergences in this sector are taken care of by properly 
choosing Z^^^±. Therefore, following j90|. we define the wave-function renormalisation for 



(6.22) 

C(7^^— part 



where ^^{l^)\cjjy-pa,rt divergent part of the Goldstone two-point functions. We 

extend the same definition in the case of the non-linear gauge, see for example Eq. lE.lfll In 
our approach, where the gauge fixing term is expressed in terms of renormalised quantities 
from the onset, this is not sufficient to make all the unphysical scalar two-point functions 
and mixing finite in the non-linear gauge. In fact in our approach and with the non- 
linear gauge, even nj^*, which does not involve 6Z^^-^ , is not finite. However as shown 
in Appendix |DJ there is a strong constraint on the two-point functions of the unphysical 
scalars. For our purpose of using these kinds of gauge fixing to check the gauge-parameter 
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independence of the results, this also is a non trivial test on the finiteness and gauge 
independence of the results. In Appendix IeI we show explicitly how one may choose to 
have finite two-point functions in the Goldstone sector at the expense of renormalising 
the gauge parameters. This method could be followed but it introduces a few extra 
renormalisation constants, which may slow the code for the cross section evaluation. In 
any case, although we can calculate the cross sections in any gauge, the gauge-parameter 
independence check is systematically applied on some random points in phase space. When 
this is passed we generally calculate the cross section in the linear gauge, with all 't Hooft- 
Feynman parameters being equal to one. In this case, linear gauge condition with ^ = 1, 
both approaches are equivalent. 



6.4 Some remarks on the explicit form of the renormahsation 
constants 

The renormalisation procedure outlined above together with the exact and complete compu- 
tation of all two-point functions permits to derive in a straightforward manner the explicit 
expressions for all parameters and wave function renormalisation constants. The complete ex- 
pressions for all two-point functions are defined in Appendix^ From the conditions imposed 
in section 16., 31 one immediately extracts all the necessary counterterms. Since the general ex- 
pressions for these are lengthy and can be read off from Appendix^ we do not list all of them 
here, but just comment on some important general features. 



6.4.1 Mass shifts and charge renormahsation 

We first verify that all counterterms to the input parameter of the physical particles, namely 
the masses of all particles, SM\y^z,hj, are gauge-parameter independent. This also applies to 
the tadpole counterterm. This constitutes a strong check on our results. 

The same holds for the charge renormalisation constant, 5Y . Although this may be derived 

1/2 1/2 

from the knowledge of SZj[^ and SZ^j^ through a Ward identity, it is easy to compute it directly. 
This is done explicitly in Appendix ^ We find the gauge-parameter independent result 

6Y = -bzil + '-^bz'il . (6.23) 

While both bZ^H and bZ^^l are gauge-parameter dependent, see Eas. l6.25l6.26l below. the above 
combination is universal. 

Sy = ^ -logM^.) - i ^lY.Q){Cvy -\o^rr?^ . (6.24) 

C\jY is defined in Eq. 14.91 



6.4.2 Wave function renormahsation constants 

Since SZ^^ and SZ^^^ are crucial for charge renormalisation and since their expressions are 
rather simple we give them explicitly. 
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We have 



5Z 



1/2 
AA 



a 
A-K 

a 
47r 



)1 2 
[Cuv - log Ml) + - _ - ^ q2 [cuv - log m) 

2(1 - a) + ^) - log M^.) + 1 - ^ ^ Q2 ^^^^ _ ^ 



where a summation on ah fermions of charge Q/ is performed. 



5Z 



yi = -f^(l-a) (Cf/y-logM, 



W 



(6.25) 



(6.26) 



This shows that although both 6ZjIj^ and ^Z^j^ are gauge-parameter dependent, the combination 
that appears in the charge renormalisation is not. Moreover, observe that the choice a = 1 gives 
a vanishing Z-A transition at one-loop. This is due to the residual U{1) gauge symmetry which 
remains after gauge fixing the charged sector, with this particular choice of the gauge parameter. 

The remaining wave functions constants are not very illuminating and involve lengthy expres- 
sions that we can extract from Appendix^ Here we just list the gauge-parameter dependence 
of those of the physical particles which can be expressed in a rather compact form as 



6Z'J^ = 5Z'J^ + ^^PiCuv-^eFo{W,W,Z)), 
bz'il = SZy^ + ^^P{Cuv-^eFo{W,W,Z)), 

6Zll^ = SZlir^ -^^{sl.a{Cuv-^eFo{A,W,W)) + clrp{Cuv-^eFo{Z,W,W))}, 

5Zll' = SZl{^ + -^^r6{Cuv-^eFo{W,W,H)) + -^{Cuv-^eFo{Z,Z,H))] . 

(6.27) 

where the quantities with correspond to the linear gauge result with all Feynman parame- 
ters set to 1. The function Fq is defined in Appendix O As known [HHl I94| FJH] the requirement 
of having the residues of the renormalised propagators of all physical particles to be unity leads 
to a (very sharp) threshold singularity in the wave function of the Higgs at the thresholds corre- 
sponding to Mh = 2Mw^2Mz- This singularity is all contained in the explicit derivative term 

1 /2 

in SZjj and is therefore gauge-parameter independent. Solutions to smooth this behaviour |96j. 
like the inclusion of the finite width of the W and Z, do exist but we have not implemented 
them yet in the present version of GRACE. Therefore when scanning over Mh it is sufficient to 
avoid these regions within IGeV around the thresholds. 
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6.5 Issues of renormalisation for unstable particles 

6.5.1 Wave function renormalisation for unstable particles: absorptive part 
and gauge dependence 

As can be explicitly seen in the previous paragraph wave function renormalisation constants, 
contrary to the counterterms for the physical parameters such as the masses and couplings, are 
gauge-parameter dependent. This is a reflection that fields are not physical observables. In fact, 
since at the end we are only interested in S-matrix observables, we could have defined an ap- 
proach where we could have done without the introduction of wave function renormalisation but 
at the expense of not having finite Green's functions^. We discussed this aspect when we defined 
the wave function renormalisation of the unphysical Goldstone particles in section 16.31 and we 
argued that the wave function renormalisation of the Goldstones, which only appear as internal 
particles, cancels out exactly. The argument applies to all particles that only appear as internal 
particles. Indeed, the wave function renormalisation constant would appear both in the cor- 
rection to the propagator and to the two vertices to which this propagator attaches. It is easy 
to see that the effect from the propagator cancels that from the vertices. The wave function 
renormalisation therefore only applies to the external particles. However there is a problem 
when one is dealing with the wave function renormalisation of unstable particles. One aspect of 
this problem was just pointed out in the previous paragraph, section 1231 related to the singu- 
larity brought about by the wave function renormalisation of the Higgs near the threshold for 
WW and ZZ production. The other problem has to do with the fact that the wave function 
renormalisation has been defined to be real so that the Lagrangian be Hermitian. But obvi- 
ously one is, for unstable particles, applying Hermitian quantum field theory to non-Hermitian 
problems 597'. 

In the standard on-shell approach, because of wave function renormalisation, loop and coun- 
terterm insertions on external legs are simply and conveniently not taken into account since they 
are thought to cancel each other. The fact is, if one insists on only using the real part of the wave 
function renormalisation then one does not completely cancel the self-energy and counterterm 
insertions on the external legs of an unstable particle. There remains in particular a contribution 
from the absorptive part of the self-energy. These absorptive parts occur for unstable massive 
particles when thresholds are crossed, they could also correspond to unphysical thresholds that 
occur in gauges where the gauge parameter is ^ / 1. 

Let us show an example where a gauge dependent threshold appears explicitly even for 
physical thresholds in the non-linear gauge. We take again a Higgs heavy enough to decay into 
a pair of W and consider the one-loop amplitude. In particular we concentrate here only on the 
eventual 6 dependence. This is contained in the two selected graphs in Fig. El consisting of the 
fish-type one-loop diagram and the HWW counterterm. 

The fish-type diagram develops a 5 dependent absorptive part which, at the amplitude 
level, does not cancel against the counterterm contribution if one insists on real wave function 
renormalisation for the Higgs as given by Eq. 16.271 The absorptive 6 dependent contribution 
from the fish would cancel on the other hand if 5Zh is defined to contain both the real and 

®In fact, as stated in section l^T!^ not all two-point functions and mixing involving the unphysical scalars 
are finite in the general non-linear gauge in our approach starting with a renormalised gauge-fixing term. 
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Figure 12: Selected 6 dependent diagrams for loop(left) and the vertex countert- 
erm(right). 




produced by GRACEFIG 



imaginary part. The inclusion of the imaginary part of the wave function renormalisation is in 
fact just a convenient short-cut. More correctly had we taken into account the Higgs self-energy 
insertion (together with the counterterm) on the external leg of the Higgs ^, see Fig. 1131 the 
5 dependence would drop completely. In the non-linear gauge we found a few examples for 
other non-linear gauge-fixing parameters and for other particles besides the Higgs, like the top 
in t bW+. 

In principle, if one is only interested in corrections to the total cross section, at the one-loop 
level it is irrelevant whether the imaginary part of the wave function renormalisation is included 
or not. This is because the effect of the wave function renormalisation is an overall multiplicative 
factor to the amplitude. For the total cross section one only needs ^e{Ttree x T^i./oop); where 
Ttree and Ti^ioop refer to the tree and one-loop amplitude respectively. Since the wave function 
renormalisation contribution, Zp for a particle p, to the one-loop amplitude is Zp Ttree, only ?R.eZp 
would be picked up. But in applications where CP violation is an issue one might need the 
correct gauge invariant one-loop amplitude. This effectively requires that the imaginary part of 
the wave function renormalisation be implemented if one is to avoid external leg insertions. This 
confirms that it is essential to require gauge independence at the level of the amplitude rather 
than in the cross section. The gauge dependence argument makes the inclusion of the absorptive 
part of the self-energy a necessity, although the gauge dependent part must cancel against other 
contributions, in all generality this should not be the case for the gauge independent part 
which should then be observable. For yet another aspect of the wave function renormalisation, 

^Some care should be exercised for external leg insertions, such as factors of 1/2 and apparent 0/0 
divisions. For detailed worked out examples see [HH] , 
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Figure 13: Self-energy insertion on the external leg of the Higgs (right) with its countert- 
erm (left) for H ^W+W'. 




see [1001 • 

6.5.2 Width implementation for resonant diagrams and gauge invariance 

Since it is difficult to reconcile the concept of an asymptotic state with an unstable particle, 
it has been argued [lUlj not to consider S-matrix elements for external unstable particles. This 
provides an easy way out to the problem we have just discussed concerning the wave function 
renormalisation, especially that such an unstable particle decays into stable particles. Nonethe- 
less the treatment of unstable particles even as internal particles poses problems. Even when 
calculating tree-level processes one has to regularise the propagator of an unstable particle if 
one is close to the resonance region. This regularisation brings in elements which are outside 
the order at which the perturbative calculation is being carried out. This problem is exac- 
erbated when performing loop calculations with unstable particles. Take a tree- level matrix 
element where part of the contribution is due to the exchange of a massive particle which, in 
the following, we will take to be a gauge boson. 



^' = 7^2+ns) (6.28) 

s is some invariant mass and M is the mass of the particle. R{s) represents the "resonant" 
contribution and the remainder T some non resonant contribution. Although there are a few 
instances where R{s) and T{s) are separately gauge independent, in the most general case this 
need not be the case. M° in Eq.ESHl is gauge independent. The residue of the pole R{M'^) and 
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the non pole remnant T{s), see below, are however independently gauge invariant jlU21llU3lllU4l 

[Tn5] 

f{s) 



Therefore the naive way of implementing a constant width T, 

although numerically regulating the resonant behaviour, breaks, in general, gauge invariance 
and can lead to disastrous results, see Ref. |lU6j for such an example. Often the use of a 
running width is made T r(s)- The latter is justified on the basis that it emerges from Dyson 
summation. However this summation, that moreover mixes the orders in perturbation theory, 
is only made on the self-energy two-point function. It therefore breaks gauge invariance since 
only one part of the total contribution to the amplitude is corrected. It has been found pi?7] . 
for tree-level process, that in fact in most cases the running width does so much worse than 
the constant width. One proposal to remedy this situation while insisting on using the running 
width was to generalise the Dyson summation to include accompanying corrections to vertices. 
This scheme, the fermion pole scheme flOHl ll()7( ll()9( [1111] . is in fact only part of the full one-loop 
calculation which is quite unpractical especially from the point of view of an automatic code. 
Moreover the scheme only takes into account fermion loops for the gauge boson propagators and 
is only meant as an effective means to correct tree- level predictions. There is also an effective 
Lagrangian approach to implement this scheme with some considerations to the unstable top 
quark and Higgs (llll I112j . 

Another scheme that is easily adapted to an automatic implementation and that can be 
carried beyond tree-level is the factorisation scheme. Starting with the gauge independent 
matrix element of Eq. 16.281 one endows the total contribution with the overall factor 
(s - M2)/(s - + iMT), such that 

A^o ^ '~^\,^ .M°. (6.31) 

s - M2 + iMT ^ ' 

Although gauge invariant, this is unsatisfactory as it puts all non-resonant contributions to zero 
close to the resonance. The pole scheme {1(121 ll()3( ll()4| ll()5j based on an expansion around the 
pole offers some insight. At tree-level only the pole term is regulated by the introduction of a 
width, while the non-resonant gauge invariant remnants are not put to zero at the resonance. 
This amounts to consider 

- + i ^(^) - + Tis)\ (6 32) 

In practice however, especially when one deals with many gauge bosons and for multi-leg pro- 
cesses, this procedure becomes untractable taking into account that there might be a clash 
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in reconstructing the invariant s from otlier momenta especially if kinematical cuts are to be 
applied [in5j . A variant of the pole scheme [ll.Sj in an effective Lagrangian approach has also 
been advocated but its effectiveness has not been fully demonstrated for processes of practical 
interest. 

The majority of the schemes we have reviewed so far have mainly been applied to tree-level 
processes although some implementations like the fermion loop scheme are attempts at including 
parts of the higher order corrections as are some of the effective Lagrangian approaches. It is 
fair to say that there is as yet no fully satisfactory solution for one-loop process. Nonetheless, 
especially for neutral current processes that require the introduction of a width, an automatic 
implementation can be performed. 

Treatment of the Z width in GRACE-loop: e+e^ — > e^e^ H as an example 



Figure 14: Contributing diagrams at tree-level in terms of the s-channel type, left panel, 
obtained from e~^e~ — > ZH , and the t-channel type from ZZ fusion. 




Let us take as an example the case of e~^e~ — > e~^e~H at one-loop |42j. At tree-level, in the 
unitary gauge, the e^e~ — > e~^e~H process is built up from an s-channel diagram originating 
from e^e~ — > ZH and a t-channel diagram which is a fusion type, see Fig. 1141 Each type 
constitutes, on its own, a gauge independent process. In fact the former (neglecting lepton 
masses) can be defined as e'^e~ — > ^'^fj,~H. This case therefore falls in the category where 
R{s) (s-channel here) and T{s) (t-channel here) of Eg. IH.281 are separately gauge invariant. In 
principle it is only the Z coupling to the outgoing lepton, in this s-channel contribution, which 
can be resonating and thus requires a finite width. Nonetheless in our code we dress both Z in 
the s-channel type diagrams with a constant Z width. We apply no width to the Z taking part 
in the ZZ-fusion diagrams. 

To help understand our implementation at one-loop it is instructive to display a selection of 
some of the contributing diagrams at one- loop. 

The introduction of a width is required only for the Z coupling to the final electron pair. This 
s-channel contribution is much smaller than the t-channel contributions for which we do not 
endow the Z propagators with a width. Building up on the implementation of the width at tree- 
level, we include a constant width to all Z propagators not circulating in a loop for the s-channel 
type diagrams. For example we add a width to all Z in graphs 349,762,4311 of Fig. 1151 For 
those one-loop diagrams with a self-energy correction to any Z propagator, represented by graph 
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Figure 15: A small selection of different classes of loop diagrams contributing to 
e+e" —>■ e^e~H . We keep the same graph numbering as that produced by the sys- 
tem. Graph 4311 belongs to the corrections from self-energies, here both the virtual and 
counterterm contributions are generated and counted as one diagram. Graph 349 shows 
a vertex correction. Both graphs are considered as s-channel resonant Higgs-strahlung 
contributions. Graph 762 represents a box correction, it is a non resonant contribution, 
which can not be deduced from e+e" — >• ZH , but applies also to the correction to the 
s-channel e'^e~ fj^^/i^H. Graph 1481 is also a box correction counted as a correction 
to the ZZ fusion. Graph 1575, Graph 1741 and Graph 1757 are fusion type corrections 
involving 77, Z^ and ZZ fusion. Graph 2607 shows a pentagon correction which also 
counts as an s-channel since it is induced for e~ — > ix^ix~H. Graph 3157 on the other 
hand is a pentagon correction that only applies to e'^e~ —>■ e'^e~H . 
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4311 in Fig. 1151 we follow a procedure along the lines described in |114j . We will show how this 
is done with a single Z exchange coupling to a fermion pair of invariant mass Sjj. 
First, it is important to keep in mind that our tree-level calculation of the s-channel is done 
by supplying the width in the Z propagator. Therefore it somehow also includes parts of the 
higher order corrections to the Z self-energy which should be subtracted when performing a 
higher order calculation. The case at hand is as simple as consisting, at tree-level, of one single 
diagram. The simplest way to exhibit this subtraction is to rewrite, the zero-th order amplitude, 
A^^*'^, before inclusion of a width, in terms of what we call the tree-level (regularised) amplitude. 



Sfj - Ml Sfj - Mi + iTzMz \ 'n-Ml J " ^ ' 

The contribution will be combined with the one-loop correction while ATz will be counted 
as being beyond one-loop. 

At the one-loop level, before the summation d la Dyson and the inclusion of any "hard" 
width, the amplitude is gauge-invariant and can be decomposed as 

^a, _ np(.,,) ^ ^ nz 

- ^ - (^-^ - - - ''i^^ ■<^'^^' 

The different contributions in are the following. The first term proportional to the tree- 

level contribution is due to the renormalised transverse part of the Z self-energy correction 
n^^, including counterterms. Such a transition is shown in Graph 4311 of Fig. 1151 The term 
n^^(s ,-) 

proportional to Az'y— — comes from the renormalised transverse part of the Z-'j self-energy, 

11^^ , with the photon attaching to the final fermion (this type is absent for neutrinos in e~^e~ 
Vgp^H ). The Rz terms combine one- loop corrections which nevertheless still exhibit a Z- 
exchange that couples to the final fermions and hence these types of diagrams can be resonant, 
an example is Graph 349 of Fig. We can write Rz = Zzu + where Zzh corresponds to 
the part containing the correction to e^e~ Z*H, while contains the corrections to the 
final Zjrj vertex. Zzh{s = M^) corresponds to e^e~ ZH and is gauge invariant at the 
pole. The term C contains all the rest which are apparently non-resonant^'', an example here is 
Graph 762 of Fig. El Both M^^'^ and M^^^ are gauge invariant. 

Our procedure, in effects, amounts to first regularising the overall propagator in Ea. l6.34l bv 
the implementation of a constant Z width and then combining the renormalised Z self-energy 



^"^Strictly speaking we, here, deal only with the pure weak corrections. In the infrared limit some of 
the QED diagrams can be resonant and require a Z width even in a loop. This is discussed in Ref. |42j . 
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part in Eq. 16.341 with the part of Eq. 16.331 Since our on-sheh renormahsation procedure is 
such that ReI\.j'^{M^) = and since = — Jmnf^(Mf ), see |115j . our prescription is to write 



Sjj -Ml + iVzMz s^j - M| + iVzMz 



(liP(.^j) - n|^(M|)) / lipisfj) 



+ {s^j-Ml)C. 



(6.35) 



The above prescription is nothing else but the factorisation procedure avoiding double counting. 
It is gauge invariant but puts the non-resonant terms to zero on resonance. In practice in the 
automatic code, we supply a constant Z width to all Z not circulating in a loop and by treating 
the one-loop ZZ self-energy contribution as in Eg. 16.351 Up to terms of order 0{Tzct) this is 
equivalent to Eq. 16.351 In particular the contribution of the C term does not vanish on reso- 
nance, since its overall factor is unity rather than the factor (sj-j — M^)/(s^j — -|- iTzMz) 
that would be present in the original factorisation prescription. 



The complex mass scheme 
This scheme stems from the very simple observation that the parameter in Eq. 16.281 could 
be taken from the outset as having an imaginary part. This imaginary part should be included 
consistently even when it enters through couplings and mixings as is the case for the electroweak 
mixing defined in Eq. 15.141 through the ratio of the mass of the (unstable) W and Z. All the 
algebraic relations will therefore be maintained in this analytic continuation, in particular gauge 
invariance is not broken. Identifying the imaginary part with the implementation of a width, at 
tree-level this scheme has been used in |116j . This simple analytical continuation is rather easily 
implemented in an automatic code for the calculation of Feynman diagrams. The roots of this 
idea 1 11 7] pre-date its first practical usage in a tree-level calculation and, in fact, emerged from 
considerations at the loop level. The suggestion is that when splitting the real bare parameters, 
as is done in Eq. 16.11 one may well take the renormalised and the counterterm parameters to 
be both complex and similarly for the fields and the wave function renormahsation as suggested 
already by our discussion in section l6.5.11 The only subtle problem now is that the renormalised 
Lagrangian is not Hermitian. This poses then the problem of perturbative unitarity and how one 
defines the Cutksoky cutting rules jll5] . Barring this issue, a full implementation of this scheme 
has been worked out at one-loop and first applied to e'^e~ — > 4/OP. Its implementation in an 
automatic code at one-loop is rather straightforward since it amounts to define the counterterms 
for the parameters of Eq. 16.11 and wave function renormahsation of Eq. 16.21 by including the 
imaginary parts of the two-point functions, whereas the usual on-shell scheme is based on taking 
only the real parts at the appropriate renormahsation scale. The use of complex masses as we 
have stressed needs to be carried out consistently by analytical continuation. This means that 
even the corresponding masses that enter the loop diagrams need to be complex. This calls for 
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the extension of the loop function, which we treat in the next chapter (Chapter [Tj), to include 
complex mass arguments. This automatically regularises some of the infrared resonant loop 
integrals, a result that had been also arrived at by the direct inclusion of the width in the 
loop so that factorisation of the infrared factor be maintained jll4j . For more details about this 
scheme applied at one- loop see Ref. As it can be inferred by this presentation this scheme 
lends itself to an easy implementation in a code for the automatic calculations of Feynman 
diagrams. 

7 Evaluation of the loop integrals 

Figure 16: General structure of the N-point loop integral. I is the loop momentum, Mi 
are the masses of the particles circulating in the loop. Pi are the external momenta. $2 is 
a combination of external momenta, see Eq. [y.jjj 

Pi 




The evaluation of the loop integrals is one of the most important ingredients of a loop 
calculation. This is also one of the most time consuming especially as the number of external 
legs increases. A generic loop integral involving external particles is depicted in Fig. 1161 The 
tensor integral of rank M corresponding to a A^-point graph that we encounter in the general 
non-linear gauge but with Feynman parameters ^ = 1 are such that M < N. The object in 
question writes in DR as 

M 
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where 



i 

Di = {l + Sif -Mf, Si = Y,Pj, 50 = 0. (7.2) 

i=i 

Mi are the internal masses, pi the incoming momenta and / the loop momentum. 

The A^-point scalar integrals correspond to M = 0. All higher rank tensors for a A^-point 
function, M > 1, can be deduced recursively from the knowledge of the A^-point (and lower) 
scalar integrals. In GRACE-loop all tensor reductions of two, three and four-point functions 
are performed by solving a system of equations obtained by taking derivatives with respect to 
the Feynman parameters. All higher orders parametric integrals corresponding to the tensor 
integrals can then be recursively derived from the scalar integral, as will be described below. 
It is important to stress that this reduction is different to what is usually done through the 
Passarino-Veltman |119j or the Brown- Feynman [120j reductions. It is also different from the 
approach of Bern, Dixon and Kosower |121j who exploit differentiation of the scalar integral with 
respect to a set of kinematical variables. 

Although the present review mainly describes the methods of one-loop calculations for up to 
2^2 processes, where only < 4 loop integrals are needed, we will describe briefly the very 
recent development in the calculation of one-loop processes with 5 external legs and in one 
instance 6 legs and in particular how the A^ = 5,6-point functions are treated in GRACE-loop. 
For A^ > 5 all integrals can be reduced to A^ = 4-point functions. 

Since the computation of the scalar integrals, especially for A^ < 4 is central let us first 
describe their implementation in GRACE-loop. Let us note that, in the intermediate stage of 
the symbolic calculation dealing with loop integrals (in n-dimension) , we extract the regulator 
constant Cuv defined in Eq. 14.91 We treat Cjjv as a parameter in the subsequent (numerical) 
stages. We regularise any infrared divergence by giving the photon a fictitious mass, A. By 
default we set this at A = 10~^^GeV. 

7.1 Scalar integrals for < 4 

The two-point integrals are implemented using simple analytical formulae and evaluated numer- 
ically. This allows to achieve a quite high precision. The scalar 3-point function and all but the 
infrared divergent 4-point scalar functions are evaluated through a call to the FF package jl22j . 
Although the FF package has been extensively used and checked by many authors, we have also 
tested its accuracy and implementation in GRACE-loop by comparing its results against our own 
numerical approach to loop integrals |123j . 

For the infrared four-point function, see Fig. I17| we supply our own optimised routines. 
A purely numerical approach would lead to instabilities and would prevent a complete and 
satisfactory cancellation of infrared divergences between these loop functions and the infrared 
factors from the real soft bremsstrahlung part. Luckily some rather simple analytical results 
have been derived in this case jll4|ll24| . These can be further simplified when the box involves 
quite separate mass scales as often occurs in e~^e~ (smallness of the electron mass). In the 
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Figure 17: The left panel shows the general configuration of the infrared four-point func- 
tions of the scattering of incoming particles with masses mi,m2 to particles with masses 
m3,m4. s,t are the usual Mandelstam variables. The right panel shows some examples 

that need to he treated carefully. In the first example, the Z can he close to the resonance. 

e- F 



mi 




s -} 



1 

M 



m2 




example e^e~ FF shown in Fig. \17\ close to the Z resonance, s ~ Aff, one needs to take 
into account the width of the Z circulating in the loop. This implementation ensures that, even 
close to the resonance, the infrared divergent part exactly cancels against the soft bremsstrahlung 
correction. The calculation of some of these photonic boxes is detailed in jll4j and |124j . 

7.2 Reduction of the tensor integrals for A^^ < 4 

The tensor integral of rank M corresponding to a A^-point graph is defined in Eq. 17.11 The use 
of Feynman's parameterisation combines all propagators such that 



1 



DoDi 



D 



T{N) [dx 



1 



(l>iXi + D2X2 + ---Do{l-i: 



Xi 



N 



50 



JV-2 

r rl l-l-Xi rl- '^^ Xi 

[dx] = / dxi / (ix2 • • • / dxTv-i- (7-3) 

J Jo Jo Jo 

Because the integrals are regulated, we first deal with the loop momenta, before handling the 
integration over the parametric variables and write 

T^^V. ..p = m) J [dx] rp. . . p, with 

M M 



4^11.. p = [A ^''^^ M<N (7.4) 



Integration over the loop momenta / is done trivially. We may write a compact formula that 
applies up to the boxes. Introducing 

N-l N-l 

A = ^ QijXiXj + ^ LiXi + Ao, Qij = Si.Sj, U = -sf + {Mf - Mg), 

i,j=l i=l 
N-l 

Ao = Ml P = -J2 ^i^i (7-5) 

i=l 

we can write 

j-{N) ^ fWT{N-n/2) with f^^) = t}f^^^-{N-n/2) 
^ ' ' (27r)"r(iV) 



j-iN) = r^N) (ri^N -n/2)P^P,-\g^,mN -l-n/2] 



Tlu) = ^^""^ (T{N-n/2)P^P,Pp-^{g^,Pp + g^pP,+g,pP^)T{N-l-n/2) 
'^iupl = ^^""^ {v{N -n/2)P^P,PpP„ + ^{g^,gp^+g^pg,^+g^^g,p)T{N -2-n/2) 

PpPu)T{N -l-n/2) 



(7.6) 



It now rests to integrate over the Feynman parameters contained in the momenta P = P{{xi}). 
We show how this is done for the box (N = 4) and triangle (A^ = 3) integrals. As pointed out 
earlier the case = 2 is straightforward and is implemented analytically, some examples are 
given in Appendix O The problem now turns into finding solutions for the parametric integrals 

4"^ = Jidxf-K(N^ (7.7) 
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4^) = l[dx]xtlogA a = 0,1. (7.8) 



solely in terms of the scalar integral I^^^ = J[dx]A (for which we use the FF package ^21)- 

The appearance of the integrals J-^'' stems from expanding Eq. 17.61 around n = 4 — 2e and 
originates from the e independent terms in e(l/e + 0(6*^'^)). In fact one only needs J*-^^ = J^^ 

and J^.^p All these integrals are derived recursively. The integral in Eg. 17.71 will be referred 
to as the parametric integral of rank-M for the A^-point function. It will then be expressed in 
terms of lower rank tensors and lower N integrals. 

7.3 Reduction of the higher rank parametric box integrals 

Let us first show how the tensor box integrals are implemented. Note that one needs 15 different 
integrals for the rank-4 box, 10 for the rank-3, 6 for the rank-2 and 3 for the rank-1. The trick 
is to use the fact that 



dx 




d 

with di = —, 1 <a + + j = M < N -1 (7.9) 



is a surface term that can be derived from the parametric integrals of the triangle. On the other 
hand expanding the partial derivative generates parametric integrals of order M + 1, beside 
parametric integrals of rank A4 — 1 and rank M. To wit. 



(/yi Of /yi 'y \ /ytf^ rift ^ ryt'^f 1 i \ 

^ j \ j jn ) 

{axl~''x^ixl^5ki + I3x'^~^ x^xl^dii + -ix'^^^^'x^x^ 5mi) (7.10) 

The term on the left-hand side can be trivially integrated and can be expressed in terms of the 
triangle integral of rank M < 3. The terms proportional to Li on the right-hand side are box 
integrals of rank M whereas the term proportional to Aq corresponds to boxes of rank M — 1. 
We combine all these terms into Ci-jkh where the first index i shows that a derivative has been 
applied on the index i . The terms proportional to Qij are boxes of rank M + 1 that we want 
to derive. In particular, to generate the highest rank integrals for the box M = 4, we apply 
Eq. I7.in1 with a = /? = 7 = 1, (M = 3). This amounts to solving a system of equations for the 
integrals 4^^;, 



^ijkl 



jnkm ' jnkl J 



[dx_ 



XiXjXi^Xi 

A2 ■ 



(7.11) 
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One can thus solve for a system of equations for the higher rank parametric integrals of 
order M + 1 in terms of (previously derived) integrals of rank M — 1 and Af for the box and 
rank M for the triangle. One drawback of this approach is that one ends up with a larger set 
of equations than needed to solve the system, especially as the rank of the parametric integral 
increases. This however also shows that one can in principle carry consistency checks. To pick 
up a system of linearly independent equations we first construct 

Ci;ikk = 3 ^ QjrJjJkk^ (7-12) 
i jn 

in order to form the set^^ 



Ci-kkk — I Ci-kkk ~ ^ki Ci;ikk ] — Qij^jkkk' ■^'^) 

For the highest rank, M = 4, this provides 3 independent sets (one for each value of k) each 
consisting of three independent integrals /j^^;. (for j = 1,2,3). Therefore one only deals with 
3 simple 3x3 matrices which help solve 9 out of 15 integrals. We may refer to this set as 
the diagonal integrals. The remaining integrals are provided by the set of the 6 independent 
equations Ci-jkk where k are all different from each other, 

Ci;jkk = -2j2Qinll!l,, i^j^k (7.14) 

n 

It is obvious that the same trick applies to solving /j^], and provides 9 out of the 10 independent 
integrals. The remaining integral in this case is provided by any Ci-jk where all indices are 
different. We also apply Eo. l7.13| to the set ij^"* and obtain C^k = Yl,j Qij^j^ ■ 

This method shows that to solve for the 15 independent integrals of rank-4 one does not deal 
with a 15 X 15 matrix. Rather the previous formulation shows that this splits into simplified 
3 03 03 06 systems of equations. The 6x6 matrix is also easy to deal with since each 
row consists of only 2 non-zero elements. For the rank-3 integrals the system of 10 equations 
decomposes into 3 03 03 01, while for the rank-2, the system of 6 equations decomposes into 
30201. 

^^This could have also been arrived at more directly had we used 

j j 

This enables to re-express the second terms on the right-hand side of Eq. 17.101 involving the Kronecker 
symbols for the case q; = /3 = 7=1, as lower rank terms and triangle integrals. Namely we can write 

j 
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7.4 log(A) terms for the box and triangle 

To extract J^^^ for the box and most of the results for the reduction of the higher rank integrals 
for the triangle, we start by giving a general representation for the logarithm. Take A^r where 
N is to just remind us that it comes from a A^-point function. We can write 



-j^ N-l 



xflogAjv = ^^-^-y ^{a,(xfx,logA;v)-xfx,-5,(logA^)} (7.15) 



Specific formulae needed for the boxes and triangles are 

2 



XilogA = 1 \d, {xa, log A) - X, + ^-(-^^•^^+ I , AT = 3, a = 1 (7.16) 
logA = i^^j9,(x,logA)-l + :^^^^^|,iV = 3, a = (7.17) 
logA = \ {g, (X, log A) - ^ + ^^L,x,+2/^I^,) I ^ ^ ^ ^ ^ (7^8) 

Eq. EH shows that j(^) can be expressed in terms of the "lower" integrals J^.-^ and I^_q 2 3' 

(3) (3) 
In turn, all J^,^^ -^^ are expressed in terms of two-point functions and the integrals i 2- 

7.5 Reduction of the higher rank parametric integrals for the 
triangle 

To generate the triangle ^^=1 2 3'^^ analogy with Eq. 17.101 we use 

9i(x^xflogA). (7.19) 
For example, for M = 3 exploiting Eq. 17.171 we get 



5,(xfcXzlogA) = ^^|l, + 2^Q 



+ l^hiY. (5i(xzx,log A) - X, (l - b^^i^y^ +{k^l) 



(7.20) 



All terms with Lj, Aq or partial derivatives are lower order terms (either in or M). Grouping 
all these as Ci-ki leads to the master equation 
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Ci-ki -"^Qijljki- (7-21) 

j 

Following the same strategy as with the box, we choose the set Ci-kk which furnishes 2 "orthog- 
onal" systems of 2 equations each. Therefore instead of handling a 4 x 4 matrix we only deal 

(3) 

with simple 2x2 matrices. Similar results are obtained for /j^ . The three needed integrals 
are arrived at by first solving for a reduced system of only two independent integrals and then 
deriving the third from a single equation. 

Finally for M = 1, one solves a system with a 2 x 2 matrix. Note that the solution of all 
these equations involves the determinant of the same 2x2 matrix, namely Qij . 

Note also that the system of equations as described here leads to analytic solutions in terms 
of the scalar integrals. In GRACE-loop we implement these analytical solutions. 



7.6 Reduction of 5- and 6-point integrals 

Five point functions are calculated as linear combinations of four point functions |125l[T^ . Vari- 
ants and new techniques have also been worked out very recentlvjlH l I127| [l28j . We will describe 
the methods that have been implemented in GRACE-loop. The reduction takes advantage of the 
fact that for > 4 not all the external momenta are linearly independent. For N = 5, the 
set of vectors {sj} with i = 1, • • • 4, see Eq. 17.21 forms an independent basis of 4-vectors, which 
allows to expand any 4-momentum, particularly the loop momentum I as 

(7.22) 

where the 4x4 matrix Qij is defined as in Eg. 17.51 with i,j = 1, • • • 4. 
From Eq. l|7.22j) we express P as 

4 

which helps define an identity between the denominators of the propagators in the 5-point 
function. Using the same notation as in Eas. 17.21 - 17. 51 we rewrite 

1 4 

Dq + M^ = -J2 {Q'%iDi -Do + Li){l ■ Sj), (7.24) 
to arrive at the identity 

4 4 

4M2 - ^ {Q-%iU)iDj -Do + Lj) = -4Do + 2 ^ (g-^).,(A - Do){l ■ sj). (7.25) 

i,j=l i,j=l 
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This demonstrates that a 5-point function with a numerator of the form N{1) = • • • is 
expressed as a sum of five box integrals, 



dH N{1) 



{2tt)H P5 



(7.26) 



where 



^^5 = I?on^- (7-27) 



Putting A^(/) = 1 one gets the reduction formula of the 5-point scalar integral to a sum of 
four-point integrals. This method has been applied to e^e~ Ve^eH 2 and e^e~ tiH j45j 
where the highest rank tensor of the pentagon is M = 2. We will refer to this technique, for 
short, as the scalar-derived reduction. 

Although this technique can be directly applied to a 5-point integral of any rank M < N = 5, 
we should note that the presence of the term I ■ sj in Eq. 17.261 raises the rank of the integral by 
one unit. This causes a superficial UV divergence for M > 3. For M = 4 the reduction requires 
the evaluation of M = 5 box diagrams, a case that is not covered by our reduction formulae 
of the tensor boxes, see section 17.31 Furthermore when this formula is used to get the matrix 
elements in a symbolic way, the resultant FORTRAN code usually becomes very lengthy. We have 
developed another algorithm for the reduction of higher rank tensors which we first applied to 
e+e" ZHH We apply the identity Eq. (f7^ to the numerator N{1). We have 

4 

N{1) = I'^H''^ ■■■1'"= = J2{Q~%j{l- Si)s'^Hf'' ■■■I'"', 

= 2 ^ -Do + U)sfl>'^ ■ ■■I'^K (7.28) 



Then 



d^NiJ) _ dH ( A - Z)o)/^2 . . . 



(27r)4i P5 2 .^/^ J (2^)4i 



On the right-hand-side the rank of the numerator is lowered by one unit, though there still 
remains a sum of 5-point functions. This reduction can be repeated until one is left with a 
scalar 5-point function and box integrals. We will, for short, refer to this method as the vector- 
derived reduction. An advantage of this method is that the final expression in FORTRAN code is 
about ten times shorter than that obtained by the previous technique. 
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Let us also very briefly describe how the 6-point functions are implemented in GRACE-loopj52j. 
In this case, we first express S5 in terms of the set of four linearly independent vectors Si, i = 
1, • • • , 4 and construct the product 

4 

I -85= E('3~')^i(^5-^^)(^-^j)' (7-30) 

hence 

4 4 
-L5+J2 (<3"')u(s5 • Si)Lj = D5-Do- J2 iQ'%is5 ■ s^){Dj - Do). (7.31) 

i,j=l i,j=l 



Combining with 



Ve = Dol[ Di, (7.32) 



we find 



dH N{1) 
(27r)4i Ve 



dH (D, - Do)N{l) ^ 1 ^ /• (Dj - Do)Nil) 



(7.33) 



This is the standard reduction of a general 6-point function to a sum of 5-point integrals. 
Moreover for further reduction of the 5-point integrals we use a combination of both the scalar- 
derived Eq. 17.221 and vector-derived Eq. 17.281 reductions for the loop tensor A'^(/) of rank-M to 
arrive at a reduction which in a compact form writes as 

M ^ M-2 M=0 M=0 

iV(0 = R^DiDj Di + ^ (7.34) 

i,j=l i 

The first term of the right-hand side is tensor of rank M — 2 which corresponds to box 
diagrams. The next two-terms correspond to scalar 5-point and 6-point function for which we 
use the algorithm of the scalar-derived reduction and the one based on Eq. 17.331 for the 6-point 
function. Once the reduction of the 5-point and 6-point tensor function have been brought down 
to box integrals with a lower tensor rank, we use the algorithm developed for < 4 as described 
in detail in section [7T1 and UT^ Infrared resonant 5-point and 6-point functions that require the 
introduction of a width in the loop integrals require more care, see|52j. 



In actual computations the matrix element for 5- and 6-point diagrams is most time con- 
suming, since as we have just seen the reduction algorithms are quite involved and go through 
various steps. For example the percentage of CPU time needed to calculate A^-point diagrams 
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is summarized in table [2 for a 2 ^ 3 and 2^4 process. One can see for instance, in the case 
of e^e~ — > VePeHH that almost as many as 2000 three-point loop integrals require 8% of the 
CPU time whereas the 74 6-point integrals require 67% of the CPU time. 



Table 2: Percentage of CPU time among various N -point loops. The number of diagrams 
in each class of N -point diagram is shown in parenthesis. "Others" stands for two-point 
functions and counterterms. 



process 


6-point 


5-point 


4-point 


3-point 


others 


e"'"e~ — >• e^e~H 




33% 


11% 


47% 


9% 






(20) 


(44) 


(348) 


(98) 


e'^e~ — > Vgpf.HH 


67% 


13% 


10% 


8% 


2% 




(74) 


(218) 


(734) 


(1804) 


(586) 



As we need a lot of computer power, it is essential to develop software engineering techniques 
in order to reduce the execution time. In GRACE we have developed a parallelised version that 
exploits message passing libraries, such as PVM or MPI jl29j . Concerning the loop calculation, it 
is efficient to distribute calculations of each Feynman diagram among many CPUs because it is 
hard to create a single executable file from too long source codes. 

Another technique must be vectorisation which can make CPU much shorter. Thanks to the 
recent rapid development of microprocessors and their easy availability even for the common 
public one should think of building or adapting codes to run as vectorised codes that can be 
quite effective. Once one has enough memory/cache, the vectorisation of the amplitudes could 
lead to a much more powerful tool. It must be said however, that vectorisation of existing codes 
is not always straightforward or even possible. For example, GRACE relies on the package FF |122j 
for the evaluation of the one-loop scalar integrals (boxes, triangles). Because of its structure 
that employs too many if ... then ... else ... end if statements, the FF package is 
not fully vectorised in the current GRACE system. 

7.7 New techniques for the loop integrals 

The reduction formalism that we have outlined both for the reduction of the A'^-point functions 
with A^ > 4 to the lower A^-point functions and the tensorial reduction, even for A^ < 4, to the 
scalar integrals involves implicitly the inverse of the determinant of the matrix Qij, the Gram 
determinant. This is most apparent in our formulation of section 177)1 see for example Eg. 17.281 
which is expressed in terms of Q~^. For kinematical configurations where DetQ is very small or 
vanishes this can lead to severe numerical instabilities. It must be said that these exceptional 
configurations are, in a Monte Carlo sampling, hardly met. Moreover, the numerical instability 
around these singularities can be cured if one reverts to quadruple, or higher, precision. It 
has been shown [T^ how this solution can be optimised with a dedicated FORTRAN library for 
multiprecision operations such that it does not require much CPU time while keeping the benefit 
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of exploiting the same standard reduction formalism. 

For < 4 and for DetQ = 0, there also exist very efficient algorithms |13H I132[ I133j that 
are amenable to an automatic computer implementation. Expansions about vanishing Gram 
determinants are also possible |1321ll341ll28j . Recently there has been a lot of activity to improve 
this aspect of the loop calculation by finding new, improved and efficient algorithms for the loop 
integrals especially to avoid this problem. One approach is, instead of reducing the system to 
the master set of the scalar integrals with < 4, to use other bases for the master integrals 
which can include a tensorial integral for example, therefore avoiding the appearance of Gram 
determinants before a numerical calculation is performed. These approaches jl27( I128( I4( H I133| 
I135j combine both an algebraic reduction with an efficient numerical implementation. Other 
approaches are, to a very large extent, essentially based on a numerical computation of all the 
loop integrals |1361 ITWl 11381 [1391 ll4Ul I141j . One example is based on the contour deformation 
of the multi-dimensional parameter integrals |14flj and may be implemented in GRACE-loop. On 
the other hand, for loops with internal massless particles as would be the case for applications 
to QCD, some powerful algebraic methods of the loop integrals are being derived [1421 1143j . Let 
us also mention that most of the methods extract the ultraviolet and infrared divergences, so 
that the set of basis integrals is amenable to an efficient implementation. 

It should be kept in mind that the majority of the new techniques have not been implemented as 
fully working codes yet nor has their robustness been tested in practical calculations, especially 
as concerns multi-leg processes. Apart from|4()( I128j which, in fact, is an extension of the 
standard reduction that has been applied successfully to e~^e~ — > 4/|^, it remains to be seen 
how the other new techniques perform when handling the complete one-loop calculation of a 
physical process, of interest for the LHC or the LC. 

8 Tests on the loop calculation 

The results of the calculations are checked by performing three kinds of tests. This concerns 
the ultraviolet and infrared finiteness as well as the gauge-parameter independence. These tests 
are performed at the level of the differential cross section before any phase space integration is 
performed for several points in phase space. These tests points are chosen at random. Usually 
for these tests one keeps all diagrams involving the couplings of the Goldstones to the light 
fermions, such as xs^^e". For these tests to be passed one works in quadruple precision. After 
these tests have been passed one can switch off these very small couplings, involving the scalars 
and the light fermions, when calculating the total (integrated) cross section and hence speeding 
up the computation time. Results of these tests on a selection of the 26 processes for the 2 — > 2 
reactions displayed in Table. 01 are made available at this web lo cation |1 44 j . This list involves 
both purely vector bosons scattering, heavy as well as massless fermions scattering into gauge 
bosons as well as a few processes involving the Higgs. Therefore, as we will see, all the ingredients 
that enter the calculation of radiative corrections in the SA4 are covered by this list. 
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8.1 Ultraviolet and infrared finiteness checks 

8.1.1 Ultraviolet finiteness 

We first check the ultraviolet finiteness of the results. This test applies to the whole set of the 
virtual one-loop diagrams. The ultraviolet finiteness test gives a result that is stable over 30 
digits when one varies the dimensional regularisation parameter Cuy defined in Eq. 14.91 This 
parameter is kept in the code as a free parameter. This parameter could then be set to 
in further computations once the finiteness test, or Cuv independence test, is passed. When 
conducting this test we regularise any infrared divergence by giving the photon a fictitious mass 
that we fix at A = 10~^^GeV. The finiteness test is carried out for a random series of the gauge 
fixing parameters that include the linear gauge as a special case. 

8.1.2 Infrared finiteness and calculation of the soft-bremsstrahlung factor 

When discussing the calculation and implementation of the loop integrals, some diagrams involv- 
ing a photon exchange require special treatment. These diagrams lead to an infrared divergence 
caused by the fact that the photon is massless so that its energy could vanish. These infrared di- 
vergences in the case of the photon, either in QED or in the electroweak theory, can be regulated 
by giving the photon a small mass A. As known the dependence in this fictitious mass cancels 
against the one contained in the soft bremsstrahlung [l45j and do not hinder the renormalisation 
procedure. For a textbook introduction see for instance|Sn] or jl46j . In the non-Abelian case 
where there is no smooth mass limit, an example being QCD, this regularisation of the infrared 
divergence by giving the gauge boson a mass would fail jl47j . In this case one can revert to 
dimensional regularisation jl48j . In GRACE, for the electroweak radiative corrections we use the 
simple trick of the fictitious photon mass to regulate the infrared divergences. 

The second test that we perform relates to the infrared finiteness by checking that when the 
virtual loop correction and bremsstrahlung contributions are added there is no dependence on 
the fictitious photon mass A. We indeed find results that are stable over 23 digits, or better, 
when varying A. 

The soft bremsstrahlung part consists of the tree-level process with an additional photon 
of very small energy, < kc, and requires the introduction of the photon mass regulator, A. 
The hard photon radiation with > kc is regular and will be discussed in section 18.31 The 
soft photon contribution is implemented in the system following an analytical result based on 
factorisation and which can be generalised to any process. The bremsstrahlung differential cross 
section factorises as 



kc is assumed sufficiently small so that the tree-level dao does not change rapidly when the 
soft photon is emitted. In some cases, for instance around a resonance, special care must be 
exercised, see for examDle |114j . The factor (5soft is completely determined from the classical 
(convection) current of a charged particle and does not involve the spin connection. Therefore 
this factor is universal and only depends on the charge Qi and momentum pi of the particles of 
the tree-level process. 
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5so{t = -e^f TT^^^E^^e^^?A7r^T^ = E^U' E^ = V^T^. (8.2) 
J\k\<k,2Ey{2TT)^ ^ ' {k.pi){k.pj) ^ 

where Si = ±1 depending on whether the particle is incoming (+1) or outgoing (—1). Very 
general expressions for Kij have been derived most elegantly in jl49j . Let us here just recall a 
few special cases and refer the reader to |114j for more details. For instance, for the diagonal 
term Ka from a charged particle with \Q\ = 1 of momentum p = (S, "p^),!?^ = rr? and P = I'/fl, 
one gets the very simple result 



Another quite useful result is the contribution, Rpair, from a pair of particle-antiparticle 
of mass m and charge ±1 in their centre-of-mass system with total energy y/s. The radiator 
function writes, with /3 = a/1 — 4m-^/s 



vr IV s/3 V1-/3/ / VA; 2/3 V1-/3 



2s/3 V V1 + /3/ V1-/3 

and 



Li2(z) = - r<iti^%^, (8.5) 



is the Spence function. This factor would represent the initial state bremsstrahlung part in 
e~^e~ processes and is usually written (for s S> m^) as 



The same factor in Eg. 18. 41 can be used as the bremsstrahlung contribution for 77 — > W~^W~ 
(m Mw)- 



8.2 Gauge-parameter independence checks 

For this check we set the value of the ultraviolet parameter Cuy to some fixed value. To tame 
the infrared divergence contained in the virtual corrections we give the photon a fictitious mass 
A = 10-^^GeV. Moreover we also set all widths to zero so that no extra gauge breaking due to 
the introduction of a width is generated. We thus choose a non-singular point in phase space, 
away from any resonance, for this check on the differential cross section. 
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Table 3: Accuracy measured by the number of digits for the gauge-parameter checks on 
the 26 processes for all five gauge parameters. The numbers that appear in the last five 
columns represent the number of digits which are stable when varying the corresponding 
gauge parameter. An empty entry means that the process does not depend on the gauge 
parameter. Only one parameter is varied at a time here. We also show the number of 
diagrams both at tree-level and at the one-loop level. The number of diagrams depends on 
the choice of the gauge parameter, for examples in some gauges some vertices are absent. 
The number of diagrams that we list corresponds to the gauge which leads to the maximum 
number of diagrams. 



3hs(Loop X Tree) 


a 


P 




e 




46 X 2 


— 


30 


— 


— 


— 


112 


X 


3 


31 


31 


31 


— 


32 


150 


X 


4 


32 


31 


31 


31 


31 


288 


X 


4 


32 


30 


30 


31 


31 


334 


X 


4 


27 


27 


30 


31 


— 


336 


X 


3 


33 


29 


31 


31 


— 


341 


X 


3 


30 


30 


31 


31 


30 


162 


X 


3 


27 


27 


28 


— 


28 


213 


X 


4 


31 


29 


30 


— 


30 


196 


X 


3 


30 


29 


31 


31 


31 


239 


X 


4 


22 


25 


29 




29 


284 


X 


4 


31 


22 


31 




32 


285 


X 


4 


29 


28 


21 


26 


30 


267 


X 


2 


24 


34 


30 






338 


X 


3 


30 


29 


31 


31 




354 


X 


4 


30 


26 


31 


31 




355 


X 


4 


30 


28 


29 


29 


31 


619 


X 


5 


22 


24 


32 




31 


657 


X 


3 




24 


31 


31 




680 


X 


5 


28 


28 


31 




31 


840 


X 


6 


26 


24 


29 


30 


29 


925 


X 


7 


27 


26 


30 


31 




823 


X 


5 


29 


25 


29 


26 


31 


830 


X 


6 




23 


24 


20 


31 


827 


X 


6 


29 


23 


22 


23 


30 


805 


X 


4 






29 


27 





processes 



VeVe 

- ti 

— > e~^e~ 

> w+w- 

> w-z° 

■ W-H'^ 
. 1^+7 

W+Z° 

77 tt 
Z'>Z^ ^ tt 
W+W- ti 
Z^H^ ^ ti 
W+W 

w+w- 
^ z^w- 



e+e~ 
e+e~ - 
e+e~ 
e+e~ 

IJ'T^- 

th- 
tb^ 



77 ^ 
ZO7- 

z^w- 
w+w- 

Z^H°- 
Z^H^ 
H^W- 



^ w+w- 
w+w- 

^Z°H° 

H^W- 
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For each process we verify that it does not depend on any of the five non-hnear gauge 
parameter of the set ( = {a,P,S,k,e). Let us remind the reader that we always work with 

= = = ^- The use of five parameters is not redundant as often these different 
parameters check complementary sets of diagrams. For example the parameter /3 is involved 
in all diagrams containing the gauge WWZ coupling and their Goldstone counterpart, whereas 
a checks WWj and 5 is implicitly present in WWH. For each parameter of the set, the first 
check is made while freezing all other four parameters to 0. 

In a second check we give, in turn, each of the remaining 4 parameters a non-zero value (we 
usually take the values (2, 3, 4, 5) for this set) so that we also check vertices and diagrams that 
involve cross terms (like a x S). In principle checking for 2 or 3 values of the gauge parameter 
should be convincing enough. We in fact go one step further and perform a comprehensive gauge- 
parameter independence test. To achieve this we generate for each non-linear gauge parameter 
Ci of the set (, the values of the loop correction to the total differential cross section as well 
as the individual contribution of each one-loop diagram g, dag for a sequence of values for C,i, 
while freezing the other parameters to a fixed value, not necessarily zero. The one-loop diagram 
contribution from each loop graph g to the fully differential cross section, is defined as 

dcTg = dag(C) = 3fie (rj'^'P • r*^^^ . (8.7) 

q-tree jg ^j^g trec-levcl amplitude summed over all tree-diagrams. Therefore the tree-level 
amplitude does not depend on any gauge parameter. Note that in many processes, some indi- 
vidual tree diagrams do depend on a gauge parameter, however after summing over all tree-level 
diagrams, the gauge-parameter independence at tree-level for any process is exact within ma- 
chine precision. Tg°°P is the one-loop amplitude contribution of a one-loop diagram g. It is not 
difficult to see, from the structure of the Feynman rules of the non-linear gauge, that for each 
2^2 process the differential cross section is a polynomial of (at most) fourth degree in the 
gauge parameter. Therefore the contribution dag of diagram g to the one-loop differential cross 
section may be written as 

da,(C) = daf + CdaW + C^daf + C^daJ^) + C^d^'^) . (8.8) 

We have therefore chosen the sequence of the five values ( = 0, ±1, ±2. For each 
contribution dag, it is a straightforward matter, given the values of dag for the five input 
C = 0, ±1, ±2, to reconstruct dcr^^'^'^'^'^^. For each set of parameters we automatically pick 
up all those diagrams that involve a dependence on the gauge parameter. The number 
of diagrams in this set depends on the parameter chosen. Different parameters involve 
different (often) complementary sets. In some cases a very large number of diagrams is 
involved. An example is ZW'^ — > ZW'^ with P 0,k ~ 1 where the set involves 601 
one-loop diagrams out of a total of 840. We then numerically verify that the (physical) 
differential cross section is independent of C 

d<7 = $:d^. = Ed^f ' (8-9) 

9 9 
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Table 4: Non-linear gauge parameter checks on a (all other parameters set to zero), for 
the differential cross section W^W~ W^W~ . For details see text. 



Graph number 






d.f) 


d4^^ 


d4"' 


2 

5 

10 
13 
33 
35 




n 


.2335514E+03 
-.1721d1dE+01 
-.4324751E+01 
-.1721616E+01 


.7789374E+03 

1 1 1 /j /AT71 1 /~1 1 

-.1171640E+01 
.8649502E+01 
-.1171640E+01 
.1048909E+01 
.1048909E+01 


.5615925E+03 
.289325dE+01 
-.4324751E+01 
.2893256E+01 
-.1048909E+01 
-.1048909E+01 


321 
322 
323 
324 
325 
326 
327 


-.3606596E+02 
.2999432E+02 


11 

.1243056E+03 
-.1197862E+03 

n 


.6056929E+03 
-.7411780E-02 
-.2864131E+01 
-.1886059E+03 
.7411780E-02 
.7411780E-02 
-.2864131E+01 


-.1534141E+04 
.2758337E+00 
-.2648726E+02 
.6165932E+03 
-.1453286E+00 
-.1453286E+00 
-.2648726E+02 


-.4615316E+04 
-.2984030E+01 
.2935139E+02 
-.3381954E+03 
.1379168E+00 
.1379168E+00 
.2935139E+02 


493 
494 
495 
496 
498 
499 




U 

-.1798684E+03 
.8331849E+02 
.8331849E+02 


.4305188E+03 
-.2608640E+03 
-.2608640E+03 
.1666370E+03 


-.2277636E+03 
.2717725E+03 
.2717725E+03 
-.3332740E+03 
-.2274438E-01 
-.2274438E-01 


-.8935366E+03 
-.9422699E+02 
-.9422699E+02 
.1666370E+03 
.2274438E-01 
.2274438E-01 


741 
743 
749 
755 

(Do 

764 


.3286920E-31 


u 

-.6573841E-31 
.6445007E+00 

.6615728E+01 

U 


-.2380925E+01 
.4734479E+00 

-.2526752E+02 


.2380925E+01 
.3853975E+01 
.1060865E+01 
.2853713E+00 

yl O y1 '7 O OT? _i_ n 1 
-.4z4 / DZiJlli-\-{ji. 

-.7116714E+01 


.OOOOOOOE+00 
-.8927045E+01 
-.2457305E+01 
-.2853713E+00 

.7139393E+02 


923 
924 




n 




-.1479291E+01 
-.8424135E+00 


-.1127685E+02 
.4331788E+03 


Max{\da^^\) 
sunii 


36.066 

.63168E-28 
.24538E-29 


179.87 

.60757E-29 
.11771E-29 


605.69 

.44209E-28 
.28841E-28 


1534.1 

.69380E-28 
.11464E-27 


4615.3 

.20116 
1.0000 




5 = 

a = 1 
5 = -1 

5 = 2 
5 = -2 

5 = 5 


Results for dag 
928.43820021286338928513117418831577 (input) 
928.43820021286338928513117432490231 (input) 
928.43820021286338928513117410983989 (input) 
928.43820021286338928513117455347002 (input) 
928.43820021286338928513117411023117 (input) 
928.43820021286338928513117695043335 (derived) 
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and therefore that 



sum, = _ =0 , z = 1,2,3,4. (8.10) 



Max(|d4 

As summarised in Table El we find a precision of at least 21 digits on all sum, for all 
the checks we have done. We usually get a much better precision when the number of 
diagrams involved in the check is smaller. 

To appreciate how this level of accuracy is arrived at after summing on all diagrams, we 
show here, see Tabled in some detail the result for W^W~ —>■ W~^W~ for the check on the 
a gauge-parameter independence (all other parameters set to zero). This is extracted from 
the web-page where we have made these checks public [.144j . This process involves some 925 
one-loop diagrams (and 7 at tree-level). Even for this particular example it is not possible 
to list all the entries of the table (that is the numerical contributions for all the diagrams) 
since they would not fit into a single page (the check on a involves some 336 diagrams), 
thus the skip (u) on some of the data. For each graph, labelled by its graph number in 
Table EJ we give all da«. We see that although individual contributions can be of the 
order of 10^, when summed up they give a total of the order of 10"^''. We also show, at the 
bottom of the table, J2g for the input values a = 0, ±1, ±2 and compare these results 
to the result obtained by setting a = 5 in Eg. 18. 81 after da^^'^''^'^'^^ have been reconstructed. 
In this example concerning a we have set the values of the other gauge parameters, f3, S, k, e 
to zero. We have also made a similar test on a allowing all other parameters non-zero. 
The same tests done on a are in turn made for all other parameters. These tests are made 
on 26 processes. More information on the check concerning W^W~ W^W~ and all 
those listed in Table El are to be found at |144j . 

One more note concerning the checks on the non-linear gauge parameter compared to a 
check one would do through the Feynman gauge parameter iw,A,z in the usual linear gauge. 
Having more parameters that clearly affect different sectors differently helps in detecting 
any possible bug. Within the hnear gauge, the usual gauge-parameter dependence is not a 
polynomial, it also involves log^ and other functions of ^. It is therefore almost impossible 
to fit the exact ^ dependence of each graph. Moreover as pointed out earlier one needs to 
prepare new libraries for handling (very) high rank tensor integrals that are not necessary 
in the Feynman gauges. 



8.3 Inclusion of hard bremsstrahlung, kc stability 

A complete 0{a) correction necessitates the inclusion of the contribution from hard pho- 
ton bremsstrahlung. Although this is a tree-level process, in most cases the total cross 
section can not be derived analytically. There is no difficulty in computing the matrix 
elements. In GRACE this is done automatically keeping all particle masses. The integra- 
tion over phase space can get tricky in many cases. In fact in most cases of interest like 
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for e'^e^ processes some care must be exercised. The reason is that though there is no 
infrared problem one still needs to very carefully control the kc dependence and also the 
collinear mass singularity. This kc dependence when combined with the one in the soft 
bremsstrahlung part (based on a analytical implementation), see section 18.1.2^ should 
cancel leaving no dependence on the cut-off kc- This would constitute another test on 
the 0{a) calculation of GRACE. The collinear mass singularity is most acute when the 
mass of the charged particle is very small compared to the typical energy scale of the 
problem, as in e~^e~ at high energies. All these problems are due to the integration over 
the propagators encountered in Eq. 18.21 For example, take the emission from the positron 
with momentum, p. This propagator is defined from 

k.p = E^{E -Pcose^) = E^P 

This becomes extremely peaked in the forward direction, cos^^ = 1. For instance, while 
the term in square bracket is of order one for cos^^ = —1 it is of order ~ 10~^^ for 
linear collider energies of SOOGeV. In GRACE , integration is done with BASES [T7| which is 
an adaptive Monte-Carlo program. For these particular cases one adapts the integration 
variables so that one fully picks up the singularities brought about by the hard photon 
collinear mass singularities. This step is therefore not as automatic as the previous ones 
in the calculation of the radiative corrections since one needs to judiciously choose the 
integration variables. For more details see jll4j . 

Stability of the result as concerns the cut-off kc is tested by varying the value of the 
cut-off kc- We take kc^ X but, usually, much smaller than the maximum energy that the 
photon can have, for example a few percent of the centre-of-mass energy in e^e~ processes. 
This is, typically, of the order of the photon energy that can be observed by a standard 
detector. In some multi-leg processes and for high energy it may be necessary to go down 
to values as low as O.lGeV or even 10~^GeV. One finds agreement within the precision 
of the Monte-Carlo which is at least better than 4 digits. From the view point of the 
computation a numerical cancellation occurs among the contribution of on the one hand, 
the virtual loop diagrams and soft photon emission, and on the other hand the hard 
photon emission. These individual contributions can be 1 ~ 10 times larger than the 
tree cross section while the full 0{a) correction, including virtual, bremmstrahlung and 
hard photon emission, is of order 1%. The individual contributions should therefore be 
calculated extremely accurately. A loss of accuracy can also be caused if there is some 
severe cancellation among loop diagrams, as due to unitarity for example. In such cases 
one reverts to quadruple precision. 

In e^e~ processes where corrections from initial state radiation can be large it is 
possible to sum up the effect of multiple emission of photons, either through a structure 
function approach (see for instance |114p or more sophisticated approaches that even 
takes into account the Pt of the photon like that of the QEDPS approach |15()j . 



m„ 



P{E + P) 



+ (1 — cos 9^ 



(8.11) 



66 



9 Checks on selected cross sections 



The previous sections have shown that the system passes highly non trivial checks for 
the calculation of the one-loop radiative corrections to SA4 processes. All those tests 
are internal tests within the system. To further establish the reliability of the system we 
have also performed comparisons with a number of one-loop electroweak calculations that 
have appeared in the literature. For all the comparisons we tune our input parameters 
to those given by the authors. Therefore one should remember that some of the results 
in the following tables are outdated due to the use of by now obsolete input parameters. 
All results refer to integrated cross sections with, in some cases, cuts on the scattering 
angle so as to avoid singularities in the forward direction. Apart from e"^e^ —>■ tt where 
a complete fully tuned comparison was conducted with high precision and includes the 
effect of hard photon radiation at 0{a), we compare the results of the virtual electroweak 
and soft photon bremsstrahlung (V + S), taking the same cut-off, kc, on the soft photon as 
specified in those references. We note in passing that since the GRACE system is adapted 
to multi-particle production, we can, contrary to some calculations, treat both the loop 
corrections and the bremsstrahlung correction within the same system. Let us also note 
that for all the processes we will consider below, we have taken the widths of all particles 
to zero since we never hit a pole. 

9.1 e+e- tt 

Table 5: Comparison of the total cross section e"'"e~ — > tt between GRACE-loop and fl51^. 
The corrections refer to the full one-loop electroweak corrections including hard photon 
radiation. 



e+e" tt 


GRACE-loop 


m 


^ = SOOGeV 






tree-level(in pb) 


0.5122751 


0.5122744 


0{a) (in pb) 


0.526371 


0.526337 


5 (in %) 


2.75163 


2.74513 


y/s = ITeV 






tree-level(in pb) 


0.1559187 


0.1559185 


0{a) (in pb) 


0.171931 


0.171916 


S (in %) 


10.2696 


10.2602 



This process is an extension of the 2-fermion production program that has been suc- 
cessfully carried at LEP/SLC. The radiative corrections to this process first appeared in |152j 
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and then in |153j . A new computation has appeared very recently |151j . A dedicated 
tuned comparison between GRACE- loop and the program topf it flSlj has recently been 
conducted at some depth including hard photon radiation and with the active partici- 
pation of the authors of |151j . Details of the comparison are to be found in |154j . Here 
we only show the excellent quality of the agreement for the total cross section including 
hard photons and we refer the reader to |154j for other comparisons concerning differ- 
ential cross sections and forward-backward asymmetries. Let us point out however that 
the comparisons at the level of the differential cross sections agree within 8 digits be- 
fore inclusion of the hard photon correction and to 7 digits when the latter are included. 
For the totally integrated cross section including hard photons this quality of agreement 
is somehow degraded but stays nonetheless excellent even at high energies. As Table El 
shows, the agreement is still better than O.lpermil. 

The authors of |151j have also conducted a tuned comparison with another independent 
calculation based on |15Hj . Practically similar conclusions to the ones presented here are 
reached, see |155j . 

9.2 e+e- ^ W+W 

Table 6: Comparison of the total cross section e~^e~ — > W'^W' between GRACE-loop and 
I15b'f . The calculation includes full one-loop electroweak corrections, but no hard photon 
radiation. 



e+e- W+W 


GRACE-loop 




v/i = 190GeV 
tree-level(in pb) 
5 (in %) 


17.8623 
-9.4923 


17.863 
-9.489 


^ = SOOGeV 
tree- level (in pb) 
5 (in %) 


6.5989 
-12.743 


6.599 
-12.74 


^/s = ITeV 
tree- level (in pb) 
6 (in %) 


2.4649 
-15.379 


2.465 
-15.375 



This process is the most important electroweak process at LEP2 and constitutes one of 
the most important reactions for the linear collider. A few independent calculations 
exist and the most recent ones agree better than the permil. To check the results given 
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by GRACE, we have set our parameters to those appearing in Table 2 of the review |156j . 
The results refer to the total cross section but without the inclusion of the hard photon 
bremsstrahlung. As we see the agreement for energies ranging from LEP2 to ITeV are 
about at least O.lpermil. 

9.3 e+e- ZH 

Table 7: Comparison of percentage correction to the total cross section e~^e~ ZH 
between GRACE-loop and fU^ . 





e~^e - 


ZH 




GRACE-loop 


P! 


V 


SOOGeV 


Mh = 


lOOGeV 


4.152S9 


4.1524 


V~s = 


SOOGeV 


Mh = 


SOOGeV 


6.90166 


6.9017 


V~s = 


lOOOGeV 


Mh = 


lOOGeV 


-2.16561 


-2.1656 


V~s = 


lOOOGeV 


Mh = 


SOOGeV 


-2.49949 


-2.4995 


V~s = 


lOOOGeV 


Mh = 


SOOGeV 


26.10942 


26.1094 


V~s = 


2000GeV 


Mh = 


lOOGeV 


-11.541S1 


-11.5414 


V~s = 


2000GeV 


Mh = 


SOOGeV 


-12.S2256 


-12.S226 


V~s = 


2000GeV 


Mh = 


SOOGeV 


11.246S0 


11.246S 



This process is an important discovery channel for an intermediate mass Higgs at a 
moderate energy linear collider and could permit to study the properties of the Higgs. 
Three independent one-loop calculations exist [HSl IHH ESj which all agree beyond the pre- 
cision of any future linear collider. A comparison was conducted against the calcula- 
tion in[nn] where one of the authors has provided us with more precise numbers than 
those appearing in Table 1 of |,93j^^. Table [7| shows that the results given by our system 

^^We thank A. Denner for providing us with the correct Mw masses used in this table. Beside the input 
given in Mw is crucial for a precise comparison. The following Mw masses have been used: Mw = 
80.231815GeV(Afe = lOOGeV), Mw = 80.159313GeV(MH = SOOGeV), Mw = 80.081409GeV(MH = 
800GeV). 
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GRACE-loop and those of [HSl agree on all digits. This means that the radiatively cor- 
rected cross sections at different energies for a Higgs mass ranging from the light to the 
heavy agree within at least 6 digits. The corrections refer to the full one-loop electroweak 
corrections but without hard photon radiation. 

9.4 77 tt 

Table 8: Comparison of the total cross section 77 tt between GRACE-loop and \15^ 
include full one-loop electroweak corrections at one-loop, but no hard (final) photon radi- 
ation. 



77 — s> tt 


GRACE-loop 




.158^ 




y/s = 350GeV 
tree-level (in pb) 
5 (in %) 


0.332477 
-6.889 


0.33248 
-6.88 


= 500GeV 
tree-level(in pb) 
6 (in %) 


0.904371 
-4.824 


0.90439 
-4.82 


v^=lTeV 
tree-level (in pb) 
S (in %) 


0.434459 
-5.633 


0.43447 
-5.63 



The comparison has been made with Table 1 of jl58| without any convolution over 
any photon spectra. Mh = 150GeV so as to avoid the Higgs resonance. As we see the 
agreement is very good, it is just limited by the precision of the numbers provided in 

9.5 77 w+w- 

The first complete calculation of the electroweak radiative corrections to 77 W^W~ 
has been performed in [ih^- Jikia has performed a full 0{a) calculation, including hard 
photon radiation |160j . Comparison has been made on the one hand with Table 1 of 
|159j without convolution over any photon spectra as well as with Table 2 of )160j . In 
both comparisons we considered the total integrated cross section but with no inclusion 
of the hard photon radiation which in any case is not treated in |159j . Because the 



70 



Table 9: Comparison for 77 — >■ W^W~ between GRACE-loop and \159^ and 
GRACE-loop and '.160^. No hard (final) photon radiation is included. When not stated 
the cross sections and corrections refer to the total cross section with no angular cut. 



77 ^ W+W~, Mw = 80.36GeV Mh = SOOGeV 


GRACE-loop 


[inni 


^ = 500GeV 






tree-level(in pb) 


77.497 


77.50 


6 (in %) 


-10.06 


-10.1 


= ITeV 






tree-level(in pb) 


79.995 


79.99 


6 (in %) 


-18.73 


-18.7 


y/s = 2TeV 






tree-level(in pb) 


80.531 


80.53 


6 (in %) 


-27.25 


-27.2 


v/i = 2TeV 60° <9 < 120° 






tree-level(in pb) 


0.39356 


0.3936 


5 (in %) 


-75.6827 


-75.6 


= 80.333GeV Mh = 250GeV 


GRACE-loop 


IT591 


= SOOGeV 






tree-level(in pb) 


77.552 


77.55 


S (in %) 


-3.376 


-3.38 


v/i = ITeV 






tree-level(in pb) 


80.049 


80.05 


S (in %) 


-7.087 


-7.08 



71 



fermionic contribution is extremely small compared to the bosonic contribution in the 
radiative correction to the total cross section, we also looked at the correction with an 
angular cut on the outgoing W as check on the fermionic correction. As can be seen from 
Table El the agreement is just limited by the precision of the numbers provided in |159j 
and [16QJ. Note that in ^160] . the correction is split between the bosonic corrections and 
the fermionic corrections. When considering the total cross section the latter are much 
too small and are below the precision with which the bosonic corrections are displayed 
in |160j . Therefore given the precision of the data in |16Uj the corrections are essentially 
given by the bosonic part for the total cross section. For the entry with the angular cut, 
the fermionic corrections are not negligible. 

9.6 67 WUe 

Table 10: Comparison of the total cross section cj — >■ between GRACE-loop and 



67 — > WVe 


GRACE-loop 


nnn 


^ = SOOGeV 






tree-level(in pb) 


36.5873 


36.587 


6 (in %) 


-12.2803 


-12.281 


y/s = 2TeV 






tree- level (in pb) 


43.9368 


43.937 


6 (in %) 


-19.0917 


-19.092 



The comparison shown in Table. ITUl is made on the total cross section (0° < < 180°) 
based on Table 5.1 of |161j . No convolution on the photon spectra is applied nor is the 
hard photon bremsstrahlung included. The agreement is rather excellent. 

9.7 67 eZ 

The comparison shown in Table. ITTl is made with Table 5.3 of |162j . No convolution on the 
photon spectra is applied nor is the hard photon bremsstrahlung included. The agreement 
is excellent. 
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Table 11: Comparison of the total cross section ej — *■ eZ between GRACE-loop and 116^ . 



cy — >■ eZ 


vjrrtiiL'l-j ±OOp 




v^ = 500GeV , 






tree-level(in pb) 


U. / Uolo 


0. / Ool 


[in /o) 


— /o.Doy 


— zo.oy 


a/S = oUUvoreV , i S S -L'y 






tree-level(in pb) 


i. /oyb 


1 wn 

i. lU 


[IB. /o) 


— ZZ.OiO 


99 


, /q — 9TpV 90° < ft < ^ (sn° 






tree-level(in pb) 


0.046201 


0.04620 


6 (in %) 


-39.529 


-39.53 


= 2TeV , 1° <9 < 179° 






tree-level(in pb) 


0.1170 


0.117 


6 (in %) 


-30.845 


-30.84 



9.8 W+W- W-W+ 

This is one of the most difficuh 2^2 processes in the SAi ever to be calculated. As 
discussed previously the number of diagrams at one- loop is of the order 1000. Moreover 
very subtle gauge cancellations take place especially as the energy of the participating 
Ws increases. The most complete calculation has been performed in |163j and the code 
is freely available at www.hep-processes.de. However, hard photon radiation is not 
included. Following jl63j we have compared our results with those of the code by requiring 
a cut on the forward-backward direction such that the integration over the scattering angle 
is over 10° < 9 < 170°. Moreover we have considered two cuts on the photon energy (for 
the bremsstrahlung part), = .05^/s and kc = .5^/s^^. The Higgs in this comparison is 
light, see |163j for a justification on this issue. Having at our disposal the code, a tuned 
comparison could be performed. We can see from Table IT^ that at centre-of-mass energy 
of the W pair of 2TeV one reaches agreement over 6 digits, 4-5 digits for -^i = 5TeV 
but "only" 3 digits agreement for y/s = lOTeV. Note that even in this case this means 
that the radiative corrections are known to about O.lpermil. This very high energy for 
the WW would probably never be reached. Moreover as the authors of |163j warn, for 
this kind of energy an integration in quadruple precision is probably already mandatory. 

^■^These cut-off photon energies are much higher than those recommended in section [8.81 However we 
stick to these values to comply with those chosen in Ref . |163 | . 
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Table 12: Comparison of the total (unpolarised) cross section W^W —>■ W between 
GRACE-loop and Uf)^ . Mh = 100 GeV. For the cuts see the text. 



w+ 


^ w^- 




GRACE-loop 


.163^ 






A;, = .05v/i 










= 2TeV 
tree- level (in pb) 
6 (in %) 


77.17067 
-21.0135 


77.17067 
-21.0135 


^/s = 5TeV 
tree-level(in pb) 
5 (in %) 


14.2443 
-57.1567 


14.2443 
-57.1556 


= lOTeV 
tree-level(in pb) 
6 (in %) 


3.644573 
-93.9942 


3.644573 
-94.0272 














x/i = 2TeV 
tree-level(in pb) 
5 (in %) 


77.17067 
-17.23988 


77.17067 
-17.23989 


y/s = 5TeV 
tree-level(in pb) 
6 (in %) 


14.24434 
-49.9736 


14.24434 
-49.9724 


v/i = lOTeV 
tree-level(in pb) 
5 (in %) 


3.644574 
-83.9247 


3.644573 
-83.9577 
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Therefore it is fair to conclude that one has for this reaction an excellent agreement. Note 
that GRACE-loop automatic calculation is the first confirmation of the result of |163j . 



75 



10 Conclusions 



Precision measurements in high-energy physics need to be matched by very accurate 
theoretical predictions. This requires that one performs calculations that go beyond 
the tree-level approximation. Moreover with the increase in energy, as will be available 
at the upcoming colliders, some important multi-particle final states beyond the usual 
2 — > 2 processes become important physical observables. Higgs production at the linear 
collider is such an example. Calculation of multi-loop and multi-leg processes involves 
the computations of thousands of diagrams, with the property that the topologies that 
appear as the order of perturbation theory increases become extremely arduous. An 
example is the treatment of the A'"-point functions beyond the box or two-loop diagrams. 
A calculation by hand is obviously no longer possible. Although a few of the steps involved 
in these computations are now done with the help of computers, it has become almost a 
necessity to perform the whole chain in the calculation of these complex cross sections by 
a computer with a minimum of human intervention so as to avoid any risk of error. We 
have in this review taken GRACE-loop as a prototype of such a fully automated system 
and described in detail its workings and performance in evaluating one-loop processes in 
the clectroweak theory. The general strategy of constructing a code like this can of course 
be applied to other systems and we have discussed some of them. We have also reviewed 
in some detail the most important modules and components that an automated system 
must have. Since there are a few fully automated systems at tree-level, the emphasis 
in this report has been on the extension to and the implementation of the one-loop 
corrections. Central to this implementation is the library for the reduction of the tensorial 
loop integrals and also the reduction of the higher A^-point functions (A^ = 5, 6) to lower 
A^-point scalar integrals. We have presented an algorithm which is now fully functional 
in GRACE-loop . 

Another crucial aspect of a fully automated calculation is the possibility of checking the 
results of its output. In an automated tree-level system this is almost trivial and is usually 
provided by switching between the Feynman gauge and the unitary gauge. As we argued 
the unitary gauge or the usual general linear gauge are not suitable at all for an 
implementation in a multi-purpose one-loop automatic code. Other important checks are 
the ultraviolet and infrared finiteness tests, but the gauge-parameter independence check is 
most powerful. This is the reason a part of this review has been devoted to renormalisation 
in the SAi within a non-linear gauge. The non-linear gauge that we exploit introduces 5 
gauge parameters. This gauge fixing modifies a large number of vertices in the bosonic 
sector but can be chosen so as to leave the propagators of all gauge bosons as simple as in 
the standard hnear 't Hooft-Feynman gauge. Technically this means that the structure of 
any one-loop A^-point function is not more involved than what it is with fcrmionic loops 
and therefore that many libraries for these functions need not be extended. To show that 
one-loop automated systems have now become completely functional and trustworthy as 
concerns the treatment of any 2 — > 2 processes we have presented conclusive tests on the 
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finiteness, both infrared and ultraviolet, and the gauge-parameter independence of the 
results pertaining to some 26 processes. These checks are verified with a precision that 
attains at least 20 digits. We have also used the new system to carry further comparison 
on radiative corrections to a few processes that have appeared in the literature. This 
selection includes heavy fermion production, vector boson and Higgs boson production in 
both e^e~ , 77 and 67 machines as well as the very challenging WW scattering process. 
For the latter we provide the first check to the complete calculation that has appeared in 
the literature. In all cases we find excellent agreement. 

This shows that the automatic system of calculating radiative corrections numerically 
has now all the ingredients to tackle 2 — >■ 3 processes confidently. In the last two years, ma- 
jor progress has been made in this area thanks to the automated system GRACE-loop and 
the package FeynArts-FeynCalc-FormCalc. Most important processes for Higgs pro- 
duction at the linear collider, e"'"e~ — ^ vj^eii e"'"e~ e^e~H e+e^ 
ZHH iSlIll!, e+e- ^ ttif HEl EI , 77 ^ tiH^ as well as e+e" ^ e+e-7 W have 
been computed and for most of them checked against each other, thanks to the automatic 
systems. GRACE-loop is now in a position to compute one-loop 2—^4 processes. We 
have in fact already presented some results pertaining to e~^e~ UeUgHH jJU] and some 
preliminary ones for the four-fermion final state e'^e'^ —>■ fi~T7^ud at one-loop Based 
on the FeynArts-PeynCalc-FormCalc package the one-loop correction for a certain class 
of e~^e^ — > 4/ has also been achieved^^ recently. 

To improve the efficiency of the system for applications to one-loop corrections for 
processes with more than 2 particles in the final state, one should seek a derivation based 
on helicity amplitudes. At tree-level this has been nicely implemented in GRACE and 
applied to processes up to 6 particles in the final state. Moreover a derivation based on 
helicity amplitudes allows the implementation of full spin-correlation, for processes when 
the final particle is unstable. At the one-loop level, another advantage is that it would 
allow the calculation by the system of processes that are not generated at tree-level, such 
as 77 — > 77, 77 —>■ ZZ, Z —>■ 37, H '-f'-f. The working version of GRACE-loop can not 
handle such processes since the one-loop correction are calculated as products of tree-level 
and one-loop matrix elements. A version of GRACE-loop which is being developed is based 
on helicity amplitudes. Preliminary results are encouraging. This new version will also 
be used as an additional test on the results given by the traditional method of squaring 
matrix elements. 
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Appendices 



A Specific form of the BRST transformations 

The action of the BRST transformations is derived as a generahsation of the usual gauge 
transformations. The ghost fields corresponding to the four gauge bosons write in terms 
of the ghosts of the SU{2) x U{1) fields as: 

A S S 

C = SwC + CwC 



= cwc^ - swc^ ■ (A.l) 



One obtains 



5BRsWi = d,c^ T ie [({A, + ^Z,) - (c^ + ^c^) 
^BRsA, = d,c^ - le {W^c- - W;c^) . (A.2) 



Likewise by considering the gauge transformation on the Higgs doublet one gets 



&RSX3 = ~f H^^" + + ^^^^'^^^'^ + ^) ' 

To find the transformation for the ghost fields, notice that the BRST transformation 
is nilpotent. For instance from ((^brs)^W^^ = one gets (5brsc*. Indeed more generally one 
has, for any group, 

5brs>1; = D^d = d^c' + g[A^, c]' ^ 5brsc^ = -g\[c, c]' (A.4) 

Care should be taken that 5brs being a fermion operator the graded Leibnitz rule applies: 
^BKs{XY) = (5BRS-^)^=t-^(^BRS^) where the minus sign applies if X has an odd number 
of ghosts or antighosts, note also that (c*)^ = 0. 
In our case this implies 

Sbrsc^ = 5brsc' = -^9 ^ijk c'c'' (A.5) 
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and thus 




+i g cw c c 



(A.6) 



The transformation for the anti-ghost field is defined through the auxihary B field of 
the gauge functions, 



B Feynman Rules 

The basic Feynman rules follow the so-called Kyoto convention!^. A particle at the 
endpoint enters the vertex. For instance, if a line is denoted as W~^, then the line shows 
either the incoming or the outgoing W~ . The momentum assigned to a particle 
is defined as inward except for the case of a ghost particle for which the momentum is 
defined along the flow of its ghost number, as will be shown in the figures. 



(A.7) 
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B.l 
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B.3 Vector- Vector-Scalar 
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B.5 Scalar-Scalar-Scalar 



Pl P2 P3 



Pi 



H H H -e- 



-M 



H 



.• 4 P3 



H X3 X3 -e- — (Mi + 2£M|-ez) 



B.6 Vector- Vector- Vector- Vector 



Pl (/^) P2 (i^) P3 (P) P4 (C^) 



_2^M-^P- + (1 - ayCw)i9^'9'"' + 9^9"') 



P2,V 



w+ w- 




Sw 



w+ w- 



2 



-2g^^g"^ + (1 - ^Vew')(/7^''i7'^'' + 9^" 9 



w+ w- w~ w~ 



83 



B.7 Vector- Vector-Scalar-Scalar 
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B.8 Scalar-Scalar-Scalar-Scalar 
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B.9 Fermion-Fermion- Vector 

Fermion mixing is not shown here. Colour for quarks is also not explicit and should be 
taken into account when appropriate. 
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B.IO Fermion-Fermion-Scalar 
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Ghost- Ghost- Scalar 
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B.13 Ghost- Ghost- Vector- Vector 
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B.14 Ghost- Ghost- Scalar- Scalar 
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C Counterterms in the ghost sector 



Since we deal specifically with processes at one-loop, there is no need to dwell on the 
renormalisation of the ghost sector. We only briefly sketch the procedure without giving 
explicit formulae for the various counterterms and the renormalisation constants. 

There is some freedom for the introduction of renormalisation constants for the ghost 
fields. We use the following convention. 



\ / ^ZZ ^ZA 







) 


_ / Zzz 




\ Zaz 








= c± 




= -c' 








= -c^ 




(C.l) 



To derive the full counterterm Lagrangian for the ghost we appeal to the auxiliary 
S-field formulation of the gauge-fixing Lagrangian £gf, see Eq. 15.171 in section 1^21 As 
stressed previously Cqf is written in terms of renormalised fields and as such does not 
induce any counterterm. However the BRST transformation are defined for bare fields. 
Therefore in order to generate the ghost Lagrangian including counterterms one needs to 
re-express Cqf in terms of bare fields to first generate Cch with bare fields. From there 
one can then derive the counterterm ghost Lagrangian. One exploits the freedom in the 
renormalisation of the 5-fields so that the combination of the B fields and gauge fields 
shows no explicit dependence on wave function renormalisation. We therefore define 

(%)-(^- ^^::T(f^) (^^^' 

The relation is just the inverse of that for gauge fields. 
Then Eq. 15.171 is 



Igf = B+d>^W- + B+iwMwX- + {h.c.) 
+B^d>^Z^ + B^izMzX^ + 
-|- non-linear gauge terms(a, /3, 5, e, k) 
+ S-linear terms. 



(C.3) 



One then makes the identifications 

B-^iwMwX- = B+i^M^X , (C.4) 
B'izMzXz = B\^MzX, + B%^MzX^ (C.5) 
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where we defined the renormahsation of the gauge parameters as 



i^. = iz{Mz/Mz)^fZzA^fz: 



— 1 



X 



(C.6) 



:-l 



is not an independent parameter but the short-hand notation given by Eq. IC.61 
Non-hnear gauge terms can be transformed in a similar way by the renormahsation of 
(a, /?, 5, e, k). Note that the terms bihnear in the B fields would get extra factors. However 
this does not affect the renormahsation program since ^brs-B = 0. This helps define the 
bare G functions. 

We obtain bare G terms by the above equations. 



Except for G^, these are the same as those in 15.151 assuming that quantities are bare 
ones. With these, one defines the bare ghost Lagrangian, that contains extra terms than 
those obtained at tree-level due to the induced mixing ^za- One then readily obtains the 
renormalised ghost Lagrangian. 



D Auxiliary Fields and Generalised Ward-Takahashi 
Identities in the unphysical scalar sector 



Although the renormahsation of this sector is not essential if one wants to arrive at finite 
S-matrix elements, the various two-point functions of the Goldstones and the longitudinal 
gauge bosons as well as their mixing are related. 

To easily derive these generalised Ward-Takahashi identities that constrain the dif- 
ferent propagators it is very useful to introduce Cqf via the auxiliary fields as done in 
Eq. 15.171 The Ward-Takahashi identities are particularly easy to derive if one works with 
the 5-fields and considers the BRST transformations on some specific Green's functions 
(vacuum expectation values of time ordered products). For the two-point function of any 
two fields (pA and 0^, we use the short-hand notation: 



For example, take the generic Green's function (c* B^) which in fact is zero (it has a 
non vanishing ghost number). Subjecting it to a BRST transformation one gets: 



QT — Q^^\^^ _|_ ^^M_Y/X^ + non — linear gauge terms 
qZ _ ()iJ'Z_^ _|_ ^^MzXg -|- non — linear gauge terms 
qA _ _|_ ^^^M_^X2, + non — linear gauge terms 



(C.7) 



{^A 0b) = (O|(T0^(x)0b(i/))|O) 



(D.8) 
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or (G' G^) = (D.9) 



where in the last part we have used the equation of motion for the i?*'s. The above relation 
leads directly to a constraint on the two-point functions of the gauge vector boson, the 
gauge-Goldstone mixing and the Goldstone two-point functions. One novelty compared 
to the usual linear gauge is that these identities also involve correlation functions with 
composite operators. Indeed if we specialise to the ZZ functions, one has with = 1 



{Gzix)Gz{y)) = 

= d^,d^^{Z,{x)ZM) + 2Mzd^{Z,{x)x3{y)) + Ml{xs{x)x3{y)) 
+ [d^^{Z,{x){H{y)xM)) + Mz{x^{x){H{y)xM))] 

Cw 

+ {^f{H{x)xz{x)H{y)xz{y)) = ^- (D.IO) 

It is important to realise that these are the full Green's function and therefore the 
external legs are not amputated. Therefore it is crucial to note that the last two terms (in 
the last two lines) do not have the double pole structure. 

This translates into the following identity, in momentum space, 

g2(nf^ - 2Mzn^^«) + Mlli^^^''] = Azz ■ (D.ll) 



Note that in the approach where the gauge-fixing Lagrangian is expressed in terms of 
renormalised quantities, Eq. ID. Ill also holds for the renormalised two-point function. In 
the linear gauge Azz = for any (tadpole contributions must be included here). An 
explicit calculation gives 



^zz = —^^{q^ - Ml){l{q' -?,Ml){Cuv - Fo{ZH)) 



+2 



q\FoiZH) - 2F,iZH)) - Mf^iCuv - F.iZH)) . (D.12) 



The functions Fq^i are defined in Eq. IG.8I We see clearly that Azz does not vanish 
in the non-linear gauge, but at the Z pole. In fact the contribution Azz can be derived 
directly from the last two terms of Eq. ID. 101 In a diagrammatic form, at one-loop, 
Eq. ID. 101 can be described as 
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(D.13) 

Calculating the new graphs explicitly confirms the identity and is a check on our 
calculation of nf^,n^^^ and n-^^^-^^. Note that for the genuine two-point functions (con- 
tributing to Azz) the loops include all possible particles including matter fields. For the 
latter one gets the same contribution as in the linear gauge. 

For the charged sector the identities go along the same line. 
For the photon, the identities give 

nf^ = at any q'^. (D.14) 

For the AZ transition we get (from (GaGz) = 0) that 

nf^ - MzTi^^^ = . (D.15) 

This identity holds at any g^, but only at one-loop thanks to the fact that the vertex 
AHx3 does not exist at tree-level. This can be easily checked by looking up the explicit 
formulae in sections lH.21 and IH.3I and section IH. 8. 11 



E Renormalising the gauge-fixing functions 

So far we have chosen to take the gauge-fixing term as being written in terms of renor- 
malised quantities. Although this is quite practical and avoids the introduction of more 
counterterms for the (unphysical) parameters entering the gauge- fixing Lagrangian, it does 
not lead to finite Green's functions, in the general case of the non-linear gauge, even when 
all Feynman parameters are set to one, ^a,z,w = 1- However all S-matrix elements are 
finite and gauge-parameter independent. Therefore as we argued, this approach of taking 
the gauge-fixing Lagrangian as renormalised from the outset is easy to implement and at 
the same time acts as a good test on our system since many divergences in a few Green's 
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functions cancel at the level of the S-matrix. If one wants to have finite Green's functions 
one also needs to consider counterterms to the gauge-fixing Lagrangian, Eq. 15.151 The 
purpose of this Appendix is to show how all two-point functions can be made finite if one 
also introduces counterterms for the gauge parameters, beside the renormalisation of the 
physical parameters (masses, electric charge) and the tadpole as well as the wave-function 
renormalisation for all fields as defined in the main text, see section lUTTl For the two-point 
functions, the difference between this approach (taking a bare gauge fixing Lagrangian) 
and that of the paper (taking the gauge fixing as renormalised) , only concerns the un- 
physical scalars (and longitudinal part of the vector bosons). Instead of Eq. 15. 151 we write 
for the charged sector the gauge fixing terms in bare quantities. 



(E.l) 



For the neutral sector one has to allow for A-Z and A-xs mixing. We take 
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CGF,iz,A) = -7^{d.Z + e^{Mz + ^eH)x/ -^{d.A + d^^MzXs? 

+ Sr]Az d.A d.Z . (E.2) 



It will also be useful to introduce and such that ^(y = y^vK iwi^'z = \jiz ^z- At 
tree-level our implementation requires that ^ ^ = ^'^ ^ = 1, 6^'^ = 6rjAz = 0. 6C,'^ and 
61] Az should be considered of order 0{a) and are introduced to avoid that Cgf,(z,a) does 
not induce any non diagonal photon transition. Apart from the S^'a and Sr)Az terms, these 
bare gauge fixing conditions are, of course, formally the same as those we introduced in 
Eq. 15.151 In fact since the gauge fixing for the photon is still linear and that we are working 
with C,A = 1; we could still take the gauge-fixing for the photon to be renormalised. As 
we will see later, the two approaches lead to the same form of the two-point functions 
for fi^^'^^ and fl'^^^. Since we only seek to show how finite two-point functions can be 
arrived at, it is sufficient to only consider the addition of the counterterms to the Feynman 
gauge parameters ^ and C,'- The bare and renormalised Feynman gauge parameters are 
related as 



l^=^^ + 5^^C = C + 5^l^ i = W,Z,A. (E.3) 

Cqp i^z,A) be derived through the auxihary B-ficld formulation and the gauge functions and G'^. 
One writes £gf,(z,a) = (W2)(S^)' + (W2)(S^)' + SfiAzB'^B^ + + B'^G'^ with G^ = d.Z + 

£,'z{Ma + (g/2cty)£ S)x^ and G^ — d.A + 5£,'axMzX3- For the purpose of generating the counterterms 
for the two-point functions at the one-loop order, one makes the identification, 5rjAz — 5f)Az / {(,aS,z) and 
5i'A^5Ux-{^'zl^z)5f,Az. 
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Since we will only work in the 't Hooft-Feynman gauge, it is sufficient to only consider 
the case with all = = 1 apart from = 61] az = 0. 

The generated counterterms for the two-point functions, allowing for the renormali- 
sation of the gauge fixing Lagrangian, and for the approach we take in GRACE-loop of 
considering the gauge fixing Lagrangian in Eq. 15.151 renormalised, are shown below. 

1. Vector- Vector 





Cgf with renormalised quantities 


Lgp with bare quantities 


WW 


n^r = 6M^ + 2{M^ - e)bz]i;' 
= ml, + 2Mi,bz]l^ 


unchanged 

= 2(M2, - q^)bZf + ml, + qHiw 


zz 


flZZ = + 2(M| - q^)5Z'-Jz 
flfz = 5Ml + '^M^dZ^Jz 


unchanged 

nf ^ = 2(M| - q^)bzTz + ml + qHiz 


ZA 


nf^ = (M| - q^)5ZYi - qHZ]l'^ 

= Mldz^l 


unchanged 

MlbZ^il - q\5Z'Jl + 5zTz - ^VAz) 


AA 


tl^^ = -2qHzTA 

nf^ = 


unchanged 

nf ^ = -2qHz]ll + qHU 


Scctlcir-^ 


jcalar 

Lgf with renormalised quantities 


Cgf with bare quantities 


HH 




unchanged 


X3X3 


fl>^3 = iqHZ]!^ + f 


= 2(g2 _ M|)<5ZV2 _ Ml5(z 
-SMI + f 


XX 


fl>^ = 2qHZ]l^ + f 





3. Vector-Scalar 





Cgf with renormalised quantities 


Cgf with bare quantities 


Wx 


tl^x = Mw{5Gw + 5Z]i^ + (5^1/2) 




ZX3 


ti^^^ = Mz{6Gz + SZ'J^ + bZ]l^) 


tlZxz = M^[^Siz-fiz) 


AX3 


fl^X3 = MzbZ^il 


n^X3 = -Mzb^A 
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4. Fermion-Fermion 

This remains the same in both approaches and is given by Eq. Ki.fil 

Note that in both approaches one has n^(0) = n^(0) for all vectors as should be. 
Note also that for the 2-point functions involving photons, n^^'^"^ and fl"^^^, 5^a, S^'a 
and 6riAz can be chosen so that the counterterms in both approaches are the same. In 
particular the Ward identities Ilf^ = and Ilf ^ - M^II^^^ = o, see Eqs. ITTT4irfTT5l are 
maintained after renormalisation. Therefore we can take 



ZA 

1/2 , <:ryl/2 



6vAz = bZ^A^^nz- (E.4) 
Let us turn to the charged sector (the Z transitions go along the same line). Defining 



A^^ = q'{flY - 2MH.n^^) + M^fl^^ , (E.5) 

we find that 

in our approach (as expected) , 

(g2 - M^) [q''6^w + M^dlw + + ^M^SZ]/^ - 2q^bZf) 

with Lqf in terms of bare fields. (E.6) 

This again means that there are constraints on ^vi/ and i-^ they are not inde- 
pendent once the other wave functions have been set. Exactly the same applies for the 
ZZ transition. Also this means that on-shell renormalisation for the unphysical sector is 
possible. That is, that the scalars have poles at the same location as the physical vector 
bosons. For example for the W, this condition {IIl{M^^) = Il^{M^r) = 0) gives that 



ACT 



M^S^w + SM^ = +n>^(M^) + ^. (E.7) 

This shows that the renormalised W-x transition becomes finite in the non-linear 
gauge: 



flWx 



2M, 



nr(M^) 



+ n^(M^ 



w 




finite 



(E. 
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This can be exphcitly shown by using the full expressions for the two-point functions at 
one-loop given in Eqs. IH.14IH.26l and IH.321 
Using Eq. IE.7[ the renormalised XX writes 



Iix = 2{q^ -M'^)(5Z]I^ + -^{h+~5-{2 + l/^^))Cuv + finite) . (E.9) 



In our approach we define 5^^^ so that all Cuv terms proportional to vanish: 



5Z]/^ - -^^^^ = t£|; + '/^^^ - ^ - 5) . (E.IO) 



This is exactly the same result had we required H^ in Eq. IE. 91 to be finite. This result 
would have been arrived at directly had we required that the residue at the pole of the 
XX propagator be 1. The Cuv Part of 5Zy^ would be the same, differences would appear 

A 

in finite terms that have no incidence on S'-matrix. 

Having constrained 5^w and 5Z}^'^, 5^w is fixed. Taking for example only the Cuv part 
of H^, from Eq. IE.7[ one has that 



-Ml - Ml - ^^^^ "i^^^^'^+^^'^ + T + xr^^- 



(E.ll) 



We indeed find, by explicit calculations, this to be verified. Similar results hold for 
the other combinations of two-point functions. 



F A library of counterterms for the vertices 

Here, we list the full counterterms to the vertices after applying the field redefinitions. 
Those for the ghost vertices are not shown since they are not required at one-loop. Those 
for the two-point functions (propagators) and the tadpole have been discussed separately 
for a proper definition of the renormalisation conditions. 

(■ ■ ■) will refer to the tree-level expression of the vertex defined in SecElbut with a = 
P = 6 = e = k = 0. As a result of Z — 7 mixing new vertices, denoted as (new), appear. 

To help write our results in a compact form, we introduce, as is done, in jHO] the 
following "correction" factors 



^CjYij 



5C 



H 



drrij 
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5Gw — 
SGz = 
SH = 



5M^ 
2M| 

2(M| - M^) 
SGi = 5Gw-5H 
5G2 = 6Gz-SH 

6Gs = 6Gz — SGw 
26M^ - 5Ml 



F.l Vector- Vector- Vector 



Pi ifj') 


P2 (J^) 


P3 (P) 




w- 




A 


{6Y + 25zKr^ + 5Z]ll){WWA) + 5Z^Jl{WWZ) 


w- 


w+ 


Z 


{5Y + 5Gi + 25Z]i^ + 5Z^Jz) {WWZ) + 5Z]lz{WWA) 



F.2 Vector- Vector- Scalar 



Pi {p) 


P2 


P3 






A 




{SY + 5Gw + 5^^' + 6Zll^ + 5Z]ll){WAx) + 5ZYl{WZx) 




Z 




{5Y + SH + 6ZU' + SZy^ + SZy^){WZx) + 5Z'J^{WAx) 


w- 


w+ 


H 


{6Y + 6G2 + SGw + 2bZ]P + bZ'^''){WWH) 


z 


z 


H 


{bY + bG2 + + bGz + 2bZ^il + 6Zll^){ZZH) 


z 


A 


H 


5Z^Jl{ZZH) (new) 



F.3 Scalar- Scalar- Vector 



Pi 


P2 






H 






{5Y + 5G2 + 5^^/' + SZl''' + 5Z^')(i/xW) 


X3 


x^ 




(5y + 5G2 + 5zV2 ^ 5^1/2 + 5z]i;'){xzxW) 


x~ 


x" 


A 


{5Y + 25Zy2 + 5ZY1){XXA) + <^4^i(xx^> 


x~ 


x+ 


Z 


(jy + + 25Zy2 + 5ZYz){xxZ) + ^^^^'(XX^} 


H 


X3 


Z 


(5r + 5G2 + + + 5ZV2 ^ szy^){Hx3Z) 


H 


X3 


A 


^Z^li^XsZ) (new) 



98 



F.4 Scalar- Scalar- Scalar 



Pl P2 P3 



H H H 


{SY + SG2 - SGw + SGh + 3SZli') - ST -^—^ 


(HHH) 


H X- X+ 


{5Y + 5G2 - SGw + SGh + SZ^' + 2bZ'l'') - ST J 


{HXX) 


H :y.3 :y3 


[SY + SG, SGyv + SGn+SZ]l^ + 2SZl'^) ST ^ ^ . 


(Hx.X-^) 



F.5 Vector- Vector- Vector- Vector 



Pl ifj') 


P2 (i^) 


P3 (P) 


Pa (o-) 




w+ 


w- 


A 


A 


{2SY + 2SZ\l^ + 2SZ)ll){WWAA) 

+2SZ^Jl{WWAZ) 


w+ 


w- 


A 


Z 


{2SY + SGi + 2SZ]l^ + SZYI + 
+SZ)lz{WWAA) + SZyi{WWZZ) 


w+ 


w- 


Z 


Z 


{2SY + 2SGi + 2SZK,^ + 2SZ^Jz){^^ ZZ) 
+2SZ]lz{WWAZ) 


w+ 


w- 


w- 


w+ 


{2SY + 25^2 + ASZll,^){WWWW) 
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F.6 Vector-Vector- Scalar- Scalar 



Pi {lA 




P3 


P4 




A 




H 




(25y + 5G2 + ^^^^i + szli^ + + 5zV2)(Aiyi/x) 


A 




X3 


x^ 


(25y + 5G2 + bzYi + + + 54/')(Aw^x3x> 
^szTa^^Wx^x) 


Z 




H 


x^ 




Z 




X3 


x^ 


(25r + + 54z + + <^^xf + <^4'")(^W^X3X) 

-r^Z'il{AWxzX) 


A 


A 




X" 


(25r + 25Zyi + 25Zy2)(AAxx) + 254/i(ZAxx) 


Z 


A 


x+ 


x~ 


{2bY + 5G4 + 54z + + 25Zy2)(ZAxx) 

+54/i(ZZxx) + (^^XX) 


Z 


Z 


x+ 


X" 


(2(5F + 2(56-4 + 2(5^^2' + 25Zy2)(zZxx) 
+2<5Zyj(ZAxx) 




w- 


H 




{2bY + 25^2 + 2(5 Zj;^^ + 2bZ]!^^)lyVWHH) 




W~ 


X3 


X3 


(2(5F + 256-2 + 2(5Zi^^ + 25Zy2)(i^i^;^3;^3) 


T T T 1 

w+ 


W~ 


X" 


1 

x+ 


(25y + 2(5^2 + 2(5Z^ + 2(5 2) (pyw^xx) 


z 


Z 


H 


H 


(2(5y + 2(56-2 + 2(5^3 + 2(5^^^ + 26Zll^){ZZHH) 


z 


Z 


Xs 


X3 


(2(5y + 2(562 + 2(563 + 2(5^^^ + 2(5Zy2)(zZx3X3) 


z 


A 


H 




dZ]!l{ZZHH) (new) 


z 


A 


X3 


X3 


^Z^za{ZZx3X2.) (new) 
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F.7 Scalar- Scalar- Scalar- Scalar 



Pi P2 P3 Pa 



H H H H 


{25Y + 25G2 

ST ^ 

swMwMl\ 


- 25Gw + SGh + ^^SZll^) 
{HHHH) 


X3 X3 X3 X3 


{25Y + 26G2 
ST ^ 


- 2SGw + 5Gh + 
(XsXsXaXs) 


X^ x^ 


{26Y + 25^2 

ST ^ 


- 2SGw + 5Gh + 

(xxxx) 


H H X3 X3 


{25Y + 25G2 
ST ^ 


- 25Gw + + 25Zl/' + 25ZV2) 
(MX3X3) 


H H X- 


{25Y + 25^2 

ST 


- 2(5Gh/ + 5Gh + 25Zl/' + 25Z]I^) 
{HHxx) 


X^ X" X3 X3 


{2SY + 2SG2 
ST ^ 


- 26Gw + 5Gh + 25ZV2 + 25Z^^{^) 
(XXX3X3) 



F.8 Fermion-Fermion- Vector 

We define L,R^{1t 75) /2. 



Pi 


P2 


P3 {lA 




7 


f 


A 


+ 5z]ll + 2SZyjl)eQfYR 


7 


f 


z 


{6Y + 5^2 + 5^3 + 54z + 26Zyj^)-^2hYL 
SG2 + 5Gs + SZ^il + 2SZyi)-^{-2sl,Qa^L) 

ZS\yCw 

SG2 + 5Gs + SZ^l + 2SZ)l^)-^{-2sl,Qa''R) 

ZS\yC\y 

+SZ]lzeQn^{L + R) 


U/D 


D/U 




- {5Y + 5G2 + SzIIj)^^, + SZll,)^^, + 8Z]l')^^YL 
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F.9 Fermion-Fermion- Scalar 

L,i? =(lT75)/2 

Pi P2 P3 

J f H 

{5Y + 5G2 + 5Graf - SGw + 6Z]'^ + 5Z]';^ + 5Z]l^) f--^^^) L 

\ Zsw J-Viw / 

+{8Y + 5G^ + 5G^f - SGw + ^Z'ff + 5Z]'^ + 5Z]l^) (-^^) R 



U/D U/D X3 

{6Y + 6G, + 5G^f - 6Gw + ^Z^)^^^ + 5Z%%^, + ) ((-/+)^^) (-^) 
+{5Y + 5G, + 5G^f - 5Gw + ^^J^J^)^ + 5Z\lj)^^^ + ) ((-/+)^^) R 



u 


D 






-i) 












D 


U 


X- (5F + 6G2 + 6GmD - 




-i) 






+{6Y + 6G2 + 6Gmu 







G Properties of two-point functions 

As mentioned earlier, loop integrals are calculated using dimensional regularisation. In 
the following / will be the loop momentum. Since one-point and two-point functions 
(tadpoles and self-energies) are essential in the derivation of the counterterms, we list 
here the properties of these functions. 

For the one-point function, which corresponds for example to the diagram shown in 
FiglTHl(a). we have the well known result: 

where Guv is defined in Eq. 14.91 

A typical two-point function refers to a diagram as shown in FiglTHl 



(G.l) 
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Figure 18: Diagrams for the one-point (a) and two-point functions (b). 

e-q 




This leads to the calculation of 

dH N{1) 



r d"i 
J if27r 



dx 



N 



(G.2) 



i{2T:Y Jo [(1 - x){& - Ml) + x{{i - qf - Ml)Y ' 
where N{1) depends, in general, on the momenta Z, q and the masses. Defining D2 as 



D2 = {l-x)M\ + xMl-x{l-x)s, {s = q^) 
we usually need to compute 

1 



r 1 
J i(27r)" (£2 - J 



{Cuv- log D2) 



I 



(27r)" (£2 - L>2)2 I67r2 

—2D2 (Cuv + ^ -log D2j 
{Cuv + 1 -log 02)9^" 



i(27r)" (£2 - D2)2 167r2 

i^'r _ 1 ^2 
i(27r)» (£2 - D2)^ ~ 167r2T 



Then the integral over the parameter x gives 



i 



'dx L>2 = ^(Ml + M|)-^s 

Z D 



(G.3) 

(G.4) 

(G.5) 
(G.6) 



(G.7) 



FJA, B) = (' dx x"logD2 = t 2;" log [(1 - x)Ml + xM^ - x(l - x)s] (G.8) 
Jo Jo J 
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We do not show the exphcit form of F„ in terms of elementary functions (the result is 
well known). We only encounter n = 0, 1, 2. 

The notation F{A, B) = Fi{A, B) - F2{A, B) = F{B, A) is sometimes used. 



F{A,B)^ CdxD^logD^ 

(G.9) 

= Ml (Fo(^, B) - F,(A, B)) + MlF,(A, B) - sF{A, B) 



We have several relations for the F„ functions as shown below. All F„ can be reduced 
into Fq. 

Exchange of A and B 

Fo(S,^)= Fo{A,B) 

F,{B,A)= Fo{A,B)~F^{A,B) (G.IO) 

F2{B,A)^ Fo{A,B)-2F^{A,B) + F2{A,B) 

Reduction into Fq, A ^ B 



F^{A, B)^l[l + Fo{A, S) + 1 {Ml log M^ - Ml log - + Ml) 

(G.ll) 



+^ [mI logMl + \{Ml - Ml)) - 

Reduction into Fq, A = B 



(G.12) 



FM,A) = \fo{A,A) (G.13) 
F,{A,A) = -{i-^\ Fq{A,A) + ^ log Ml - - (G.14) 



G functions (Derivative of F) 



Gn{A, B) = ^F^{A, B) = r dx "^"'^^^ (G.15) 
as Jo JJo 
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F and G at special energy 



F4A,B;C) = F4A,B 
Fn{A,B;0)=Fn{A,B)l^, 
Gn{A,B;C)=Gn{A,B)l^j,,, 
Gn{A,B;0)=Gn{A,B)\^^, 



(G.16) 



Fn for s = 



FoiA,B,0) 



log Ml 

M| log M|- Ml log Ml 
Ml-Ml 



(A = 5) 
1 (A^B) 



(G.17) 



H Results for the one-loop corrections to the propa- 
gators 

We here give the details on the calculation of the loop corrections to the various propaga- 
tors and mixings. For the vector bosons we present both the transverse and longitudinal 
part as defined in Section 16.31 We also show the various contributions by classifying 
them according to the diagrams of Fig. ^1 Therefore for each propagator we show a 
table containing the two types of diagrams. The fermion contributions are summed over 
all fermion species in the case of the neutral sector and over all doublets in the charged 
sector. In both cases summing over colour for quarks is implied. Before presenting the 
results for the two-point function, we first start by presenting the one-loop contribution 
to the tadpole. Although we will require this contribution to vanish against the tadpole 
counterterm we give its full expression for completeness. Also, the latter is needed for the 
Ward identities. 



H.l The tadpole 

The tadpole contribution only receives contributions from diagrams of the type 

shown in Fig. ITHl(a) where A = W, Z, x, Xs; H, c, c^, /. The result is as follows: 
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( {Cuv - logM^ + 1)(3M^ + ^M^) 



+M| ({Cuv - log Ml + 1)(^M| + ^M|) - M|) 



2M, 



(H.l) 



+ -M^(Cc/y-logM| + l) 

-J2^mjiCuv-logm} + l) 

f 

It is important to note that all dependence on the non-linear gauge parameters (namely 
e and S) vanishes among all diagrams and is therefore the same as in the usual linear gauge. 
This can be considered as a check on the calculation, since the tadpole T can be considered 
basic parameter of the theory. 

From this expression we immediately derive the tadpole counterterm: 



^rp rploop 



(H.2) 



H.2 A- A 

For all the two-point functions we will list, as done in the table below, the contributing 
diagrams where (a) corresponds to the type shown in Fig. ITHl-fa) and (b) to Fig. [THJ-(b) 



(a) 



(A B) = {W, W), {W, x), ix, W), ix, X), (c+, c+), (c-, C-), (/, /) 
A = W,x, c+, c- 



a 



7Cuv - 5Fo{W, W) - 12F{W, W) + - - 4{1 - a){Cuv - Fo{W, W)) 



- 8j2Q}[^Cuv-F{fJ) 



[R.3) 



Note that independently of the gauge parameter 11^"^ (0) = 0. This is just a remnant 
of the QED gauge invariance which is explicit at one-loop. This also gives nf'^(g^) = 
which we explicitly verify. More generally we will also check explicitly that in both the 
linear and non-linear gauges 11^(0) = 11^(0) for all vector- vector transitions. This is 
another check on the calculation and encodes the property that there is no spurious pole 
in the propagators essential for the Goldstone mechanism. 
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H.3 Z-A 



(b) 



(A B) = (W, W), (W, x), (X, W), (x, X), (c+, c+), (c-, C-), (/, /) 



(a) 



A = W,x, c+, c- 



n; 



ZA 



a cw 

47r Siy 



uv 



7 + 



4 3 



F{W, W) + 



2Cyy I 3 



5 + 



-2(1 -P) {Cuv - FoiW,W)) 

-J- E IQ/I (l - ^^IQfl^w) ilCuv - F{f, f) 



9„2 



Fo{W,W) 



-2(1 - _ Ml) {Cuv - Fo{W, W)) 



n 



ZA 



a Cw ^ 
2tt sw 



a)Ml{Cuv-F^{W,W)) . 



(H.4) 



(H.5) 



Note that we do get as a check that nf"^(0) = nf^(O). Moreover for a = 1 this condition 
is even stronger since we get nf"^(0) = Ilf"^ = 0. This is due to the fact that for this 
particular choice of the parameter, the gauge- fixing in the charged sector which contributes 
here (note that fermions do not contribute to Hf"^), there is an additional f/(l)QED gauge 
invariance. This choice is therefore very useful. As we will see this is also responsible for 
the vanishing of the induced A — xs transition, see section IH.8.11 It is also important to 
remark that at the Z-pole the a dependence vanishes. This is also responsible for the fact 
that the counterterms needed for the mass definitions do not depend on the gauge fixing. 

H.4 Z -Z 



(b) 



(A B) = (ly, w), (ly, x), (x, w), (x, x), {h, x), {h, z), (c+, c+), (c-, c"), (/, /) 



a 



zz _ 

47rs^c^ 



rjnZZ ^ rjnZZ ^ _ ^) (^^2 _ ^2\^ ^rj^ZZ 



(H.6) 



ZZ 



uv 



6 



- 2M^ - Ml 



2 2 A T 
3^ 12 



8g^ct^.Fo(iy, W) + q\FQ{W, W) - AF{W, W)) 3c^ + 



4c- 



w 
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+ 2M^Fo(W, W) + Z) - -fFo{H, Z) - ^ ^ F,{H, Z) 

+ M|Fo(//,Z) + ^(M|logM| + M|logM|), 



rpZZ 



] E[((l - ^\Qf\s'wf + l) (^^^c/y - F{f, /)) - m) {Cuv - Fo(/, /)) 



(H.7) 



We have, 



nLnLinear(^l)=nLeax(^l)- 



(H.8) 



n 



zz 

L 



a 



1677 S^irC^ 



Cuve^ - 3 + 2F(i/, Z) - (1 - efFo{H, Z) - AF^iH, Z) 



+A{l-i)Fi{H, Z) 

- AMliCuv - FoiH, Z)) - 2M]jFo{H, Z) - 2(M| - Ml)F^{H, Z) 

+ {Ml log Ml + M| log Ml) - 8M^ {Cuv - Fo{W, W)) (l - 2c^(l - ^)) 

+ 2Y,m]{Cuv-Fo{fJ))\. (H.9) 
/ J 

It is easy to see that, for any choice of the gauge parameters, nf^(O) = 11^^(0) which 
is a check on the calculation. Also note that the e dependence is proportional to (f. 



H.5 W -W 



(b) 



(a) 



(A B) = (Z, W), {Z, x), (A W^), {A, x), (if, X), (ii, W), (X3, X), 
(c^c+),(c^c-),(c^c+),(c^c-), (/,/') 



A = A^,W^,i^,X3,X,c+,c- 



^ 



WW 



a 



(H.IO) 



Ct/y(yg' + 2M^-M| 



6 
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rpWW 



w 



q\F{A, W) - Fo{A, W)) - M^F,{A, W) 



Aq\F{Z, W) - Fo{Z, W)) - 4(M^ - M|)Fi(Z, W) 



Ml 



+ {-7M'z + ^)Fo{Z,W) 



1 r 



MlF^{H, W) + M|Fo(Z, W) + (M^ - Ml)F^{H, W) 
+{M^-Ml)F,{Z, W) 

+^ log + ^ log Ml + ^ log Mj, + 2M^ log M| , 



(H.ll) 



-l E |4(^C^v-F(/,/'))?'-K + m'/)Cc;y + 2mjFi(/',/) 
+2m^2Fi (/,/')} , 



^rpWW 



-2 - Fo{A, W)) + C^/3(Cc;y - Fo{Z, W))) . 

Here also we check that , 



(H.12) 
(H.13) 



ttT 

NonLinear V 



H 



WW 
L 



a 



47rs^ 



Mi 



6 



+ 2s^ q\6FiA,W)-Fo{A,W))-2M^F,iA,W) 



+ 



+ 



+ cw 



2q\6F{Z, W) - Fo{Z, W)) - 4(M^ - Ml)F^{Z, W) 



+ (^-7Ml)Fo{Z,W) 



■-w 



+ - 



Fo{H,W) 



-{Ml-Ml)F^{H, W) 



q\?>F{Z,W)- 



Fo{Z,W). 



) + 2M^Fo{H, W) - M^hHH, W) 



M|Fo(Z, W) - {M^ - Ml)F,{Z, W) 



^ A/f2 M^ 

^ log + ^ log Ml + ^ log Ml + 2M^ log M| 



+ 2s^a 



3g2(Fo(A, Py) - 2Fi(A W)) + M^{Cuv - F^{A, W)) 



+ 2c^^ 3g2(Fo(Z,VF)-2Fi(Z,VF)) + M^((:7c;v-Fo(Z,VF)) 
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4- 






+ 


clrPW^CuV- 


-QFq{Z,W) + 2Fi{Z,W)-2) 


+ 


^q'{Fo{H,W)- 


-2Fi{H,W)) + '^q\Fo{Z,W)-2F,{Z,W)) 


+ 




{H,W)) + ^q\Cuv-Fo{Z,W)) 


+ 


\ E {K + m])Cuv - 2m}F,if, f) - 2m'fF,(f, f)}\. 

doublet ) 



We again have that, for any choice of the gauge parameters, n^^(O) = n^^(O) which 
is a check on the calculation. Also note that the S, k dependence is proportional to q^, 
as is any dependence quadratic in ct, j3. All these dependencies will be present in the 
propagators/mixings of the Goldstones. 

Note that for both the WW and ZZ transition the tadpole contribution is not included 
as it will be canceled against that of the tadpole counterterms. Moreover note that such 
contribution do not depend on the gauge parameter. On the other hand the inclusion 
of the tadpole one-loop correction is needed for the Ward identities relating the bosonic 
two-point functions. 



H.6 H-H 



(b) 


(A, B) = {W, W), {W, x), (X, W), (x, X), {Z, X3), (X3, Xs), {Z, Z), {H, H), 
(c+,c+),(c-,c-),(c^c^),(/,/) 


(a) 


A = W,Z, H,x,X3,c-^,c-,c^ 



Here we explicitly add the tadpole contribution. 



with 



Il^{q') + 



3ST 



a 



w 



(nf + nf + (q^ - M^)n£ 



(H.15) 



H 



H 




Fo{H,H) 



-q^ + 3M| + 



My 



+ 3Mi (1-logM^) 
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n? = E^{y(^t/v-Fo(/,/)) + 2mJ(l-logmJ + Fo(/,/))|, 



n 



H 



-Cuv + Fo{W, W)) S + i-Cuv + FoiZ, Z)) 



9„2 



(H.16) 
(H.17) 
(H.18) 



Again at = Mfj the self-energy is independent of the gauge parameter, which means 
that the shift in the Higgs mass will also not depend on the gauge parameters. 



H.7 /-/ 

At one-loop the result is the same as in the linear gauge, but we give here the full result 
that includes mass effects as well as the contribution of the Goldstones. We have neglected 
all fermion mixing. The Kj have been introduced in section 16.21 and correspond to the 
different Lorentz structures of the fermion propagator. Since we are neglecting mixing 
and assuming CV -invariance = holds. We have also found it convenient to express 
each of these Lorentz coefficients in a basis that corresponds to the various contributions 
to the self energy (photon exchange, W-exchange, etc.). 



(b) 



(a) 



(A, B) = (/, A), (/, Z), if, W), (/, H), if, X3), (/', X) 



None 



\Qf\Kf^'^ + 



9^2 



K 

2 ^2 3 



Z(3) 



m 



AJ2 2 /l#2 



a = 1,7,57), 



:h.i9) 



= mf[-ACuv + 2 + 4Fo(/, A)], 
K^ = Cuv-l-2F,if,A), 



;h.2o) 



i^f = mf[-4Cuv + 2 + 4Fo(/, Z)], Kt''> = K 



K = 0' 



Z(2) 



.Z{1) 
1 , 



Z(3) 



K^^'^ = Cuv-l-2F,if,Z), irfW = irf«, irf (3) ^ 



:H.21) 
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= 0, 

K^^Cuv-l-2Fr{f,W), 



(H.22) 



= nif 



m 



-Foif, H) + Fo(/, Z) - 2^{Cuv - Fo{f, W)) 



m, 



= Cuv - F,{f, H) - F,{f, Z) + Ul + '^] {Cuv - 2F,{f', W)), (H.23) 

Z y Til J J 



m' 



rrit 



/2' 



{Cuv-2F,{f, W)). 



H.8 The Goldstone sector 

We do not need to be explicit about the renormalisation of this sector in order to arrive 
at finite S-matrix elements. Nonetheless we list all the vector-Goldstone mixings and 
Goldstone propagators. 

H.8.1 A-X3 



(b) 



(a) 



(A B) = {W, x), (X, W), (c+, c+), (c-, C-), (/, /) 



None 



There is no fermionic contribution. 



27rs 



w 



1-a) {Cuv - FoiW,W)) 



(H.24) 



As expected this does vanish for a = 1 and is a remnant of the ^7(1)qed gauge invariance 
for this value of the parameter. 



H.8.2 Z-Xs 



(a) 



{A, B) = {W, x), (X, W), {H, X3), {H, Z), (c+, c+), (c", c"), (/, /) 
None 
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■ 

+ 



4(-3 + 4c^(l - /3) - (Cc;y - Fo(W^, W)) 



M2 

^[Fo(i/,Z)-2Fi(i/,Z)] 



2M| 



F,{H,Z)-Fo{H,Z) 



% - Fi(i^, Z) + ^{Cuv - Fo{H, Z) 



2M| 



+ e\Cuv-F^{H,Z)) 



m. 



(H.25) 



H.8.3 W-x 



(b) 



(A, = {H, W), {H, x), {Z, x), {Z, W), {A, x), {A, W), 

{c^,c),{c^c+),{c^,c-),{fj') 

None 



n 



Wx+ 



w 



uv 



f M"^ ~ ~ ~ ~ ~\ 

+ Cuv I + + (5 + 252 + s\y{\'&d? - 12a) + 4/5(4 - 3c^) + ISc^/?^ J 

+ 4s^(4Fi(A,W^)-Fo(AW^)) 

+ 4s2^,a (-6Fi(A, ly) + 6Fo(A, M^) - 5aFo(A, W) + (5Fi(A, W^)) 



+ ly) ( s^(-16 - ^) - 2« - 24c^/? + 4c^/32 



+ Fo(Z,iy) 2-4c^ + 4 



16/3 + 24c^/3 - 204;3' 



+ 



(2Fi(iy, //) - Fo(W^, E){\ + 5)) + 2 (2Fo(W^, //) - //)) 



+ 25(Fi(l^,iJ) -Fo(l^,i/)(l + 5) 
+ 8as^(l-a)-8/3(l-c^(l-;3)) 



+ 4 E 



12 



Cuv-(mjFo{f',f)+mfF,{f,f 



(H.26) 
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H.8.4 X3 - X3 



(b) 



(A B) = {W, x), (X, W), {H, Z), {H, xs), (c+, c+), (c", c"), (/, /) 



(a) 



A = W,Z,X3,X,H, c+,c-,c' 



The second line of II^^ shows the fermionic contributions. Since the tadpole contri- 
bution appears with n-*^^, we present the formula for the sum. We note that C^^, the 
coefficient for the divergent part, is proportional to s in the linear gauge. 



V 



a 



■X3 



m. 



+VE77f (^t/v-^o(/,/)) 



(H.27) 



= - 



(H.28) 
(H.29) 



2 \2 



= ( 2Mi - Ml + 2,^ - - 2e(M^ + g^) - 3e-M| ) , 



;h.3o) 



If 



- Ml \ \ogMl + (m| - Ml) logMl - + 2M| - M 



2 \2 



H.31) 



H.8.5 x-X 



(b) 



(A S) = {E, W), {H, x), (X3, W), {Z, x), {Z, W), {A, x), (A W), (c^ c+), (c^ c"), (/, /') 



(a) 



A^A,Z,W,H,x3,X,c-^,c- 



The second line of 11-^ shows the fermionic contribution. Since the tadpole contribution 
appears with 11'^, we present the formula for the sum. We note that C^, the coefficient 
for the divergent part, is proportional to s in the linear gauge. 
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a 



C^Cuv + d^zwFoiZ, W) + djj^FoiH, W) + d^^F^iA, W) + ^ 



. . rrr . 



fH. 



(H.33) 



d^w 



2-8sl, + — ]q^- 16c^M, 
V '^w J 



c 



w 



-32Ps^wM^ - lep^c'wM^ - 2kq^ + k^M^, 



dlw = -M^ + 2Mij - 



- 25{s + Mfj) - 3S'M^, 



23M^ 



(H.34) 
(H.35) 
(H.36) 



4= logM^(-M| + 2M^-M|)+logMi(i^-M| 



+ log M| ( ^ + 8s^M| - Ml] - + 2Ml + 6M^ - ^ - 6M| 



(H.37) 



I Direct determination of the charge counterterm 

A direct derivation of the charge counterterm SY necessitates the calculation of the vertex 
e^e~ A with on-shcll electrons and in the Thomson limit where the photon momentum 
g ^ 0, r^(0). In this limit one can relate the vertex to the electron self-energy as in 
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Figure 19: Electron self energy and e^e A vertex. 



{a) A,Z,H,Xi {b) ve 




depicted Fig^J [p is the electron momentum and / is the integration momentum). Due 
to the identities 



d 



-1 



-1 



-1 



and 



dpn \ ^+ / — m 



d 



-r 



'-)+/ — m ^+ / — m 

\ -2(j9 + £)^ 



dpf, \ip + ey - A/p J {{p + ey-M^y 



(i.i) 



(1.2) 



The major part of r^(0) can be calculated by the corresponding fermion self energy. Using 
the notation in FiglTI^ 



(a) AZ,H,X3 r:f(0) = (-e) ^S(p2 



(b) W,x r:f(0) = (-e) 



(1.3) 



(1.4) 



It is only which is gauge-parameter dependent, in fact it only depends on a and 
vanishes for a = 1. 



l-a)iCuv-logM^)j>'L. 



(1.5) 
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We note that for this particular value of the gauge parameter there is a residual 
[/(1)qed symmetry and therefore it is no wonder that the naive Ward identity is verified 
in this case in Eq. 11.41 

The gauge-parameter independent part of r^(0) is derived from the (on-shell) electron 
self-energy. We obtain 



d 



{2m,K[{ml) + 2mlK'(ml) + K,{ml))Y + K,,{ml)^>^^, 



(1.6) 



Here Eq ]6.20l is used. 

The counterterm for the e~^e~'A vertex, r(f(0), is defined in Sec lF.81 Adding this to 
the loop calculation we get 



:-eY)[SY + 6Z]/^-'-^6Zl/f] + 

Cw 



6Z'£yl + G^, 



(1.7) 



The second term (within square brackets) vanishes identically. Imposing the renormali- 
sation condition (Ea KT^i f ^(0) = 0, 



6Y = -^zTa + '-^SZl^l 



Cw 



and we find the linear gauge result 



(1.9) 



J Graph theory and optimization in the generation 
of Feynman diagrams 

We shall first introduce some technical terms. In the following let us call a vertex or 
an external particle a node. Similarly, let an edge be a connection between two nodes, 
which may be a propagator or a connection between a vertex and an external particle. 
Thus an edge is expressed by a pair of two nodes (which are connected by the edge). The 
graph generation process is to construct edges in all possible ways. Although vertices 
of the same kind are not distinguished from each other, they will be distinguished in a 
program, usually through a sequence of numbered labels. Since a Feynman graph is a 
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topological object, it is independent of the way one assigns the sequence of numbers to 
nodes. Therefore two graphs are topologically the same when there is a permutation in 
the sequential numbers in a graph which produces another graph. Some permutations 
acting on a graph will keep the graph unchanged. These permutations form a symmetry 
group of the graph whose number of elements is needed to calculate the symmetry factor 
required for the calculation of the Feynman amplitudes. 

Let us consider three-loop vacuum graphs in the 0^ model for example. There are 
two one-particle irreducible Feynman graphs as shown in FigEUl A permutation of any 
two nodes in the graph of FigEOK produces the same graph. However in graph Figl2()b. 
the exchange of nodes 1 and 2 gives different numbering of the topologically equivalent 
graph, while the simultaneous exchange of nodes 1^2 and 3 4 results in the same 
graph with the same numbering. When a permutation of nodes in a graph produces a 
different representation of the graph, one must remove either the original graph or the 
new one obtained by the permutation. When the permutation keeps a graph invariant, it 
is found to be an element of the symmetry group of the graph. This symmetry structure 

Figure 20: Three-loop vacuum graphs in (f^ model 




a b 



of a graph originates from the fact that a vertex is symmetric under the exchange of the 
three legs. When the vertex does not have such a symmetry as in the case of QED, it 
is possible to avoid the duplication of graphs at all orders of the coupling constant [TH]. 
Generation of tree-level graphs is also easy since a node is distinguished from others by the 
set of momenta entering the node. The most difficult cases are vacuum graphs of the 0^ 
or 0^ model at the multi-loop level, since all nodes in a graph are equivalent. Fortunately, 
the calculation of decay or scattering processes involves external particles that make the 
distinction simpler. This also makes the code for the generation of graphs simpler and 
faster. 

The first idea to avoid duplication is to keep already generated graphs in the main memory 
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or the hard disks in order to compare them with the newly generated ones P?] . However as 
the number of graphs grows as the factorial of the number of nodes, it is impossible to keep 
all the graphs in the main memory. One is then forced to move all the generated graphs 
onto the hard disk. This very much slows down the code. A more efficient method has 
been developed by graph theorists to avoid generating duplicated graphs [HZI. It checks 
a generated graph with the one obtained through permutation and therefore requires 
memory space just for two graphs. The outline of this method is the following: 

1. Let P be the set of all permutations among nodes of a graph G and <j{G) be a graph 
obtained by applying a permutation a G P to a graph G. The set of all topologically 
equivalent graphs of G is {(y{G) \ a E P}. 

2. We introduce the connection matrix M of a graph G, whose matrix element Mij is 
1 when nodes i and j are connected by an edge and is otherwise. We apply the 
permutation a to this matrix. When the permuted matrix of a graph is identical to 
the matrix of another graph, these two graphs are topologically the same. 

3. By lexicographical comparison of matrix elements among two connection matrices 
we can introduce an ordering relation, among graphs. Among all topologically 
identical graphs, we keep only the "larger" graph in this comparison . This means 
that we discard graph G when there exists a permutation a such that cr(G) >- G. 

4. The resulting algorithm is as follows: 

(a) When a graph G is generated, generate all possible permutations a of the 
nodes. 

(b) The connection matrix of G is compared with one of a{G). 

(c) If cr(G') )~ G, the graph G is discarded. Otherwise another permutation is tried 
in the same way. 

(d) When G is found to satisfy G >z <^{G) for all possible permutations, G is kept. 

(e) At the same time, the number of such permutations a is counted that satisfy 
G = a{G). 

5. The symmetry factor of G is given by l/(iVs x Ne)., where Ne is the number of 
permutations among edges that keep G invariant. 

This method was first applied to the generation of Feynman graphs in the code 
QGRAF[21 . Even with this method, it is necessary to compare two graphs. Unfortu- 
nately, no efficient algorithm which avoids the factorial growth from all permutations and 
that terminates within a number of steps which is a polynomial function of the size of 
the graph is known The size is of course determined by the number of nodes. In 

^^Such efficient algorithms in graph theory are known as polynomial-time algorithms. 
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order to speed up the code, one has to decrease the number of generated graphs and the 
number of permutations to be tested in a consistent way^^. A method for acceleration by 
a systematic classification of graphs proposed in (Ml is used in the GRACE system. 

Programs can be checked by two counting methods First the number of generated 
graphs for the (fr" or cf)'^ model is compared with what is predicted by a graph theoretical 
method. Second the sum of the number of graphs weighted by the symmetry factor is 
compared with the value calculated in zero-dimensional field theory. 



QGRAPH seems to include some acceleration method but this is not described in 



120 



References 

[1] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. 
Shimizu, Nucl. Phys. (Proc. Suppl.) 116 (2003) 353; hep-ph/0211268. 

[2] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. 
Shimizu, Phys. Lett. B559 (2003) 252; hep-ph/0212261. 

[3] A. Denner, S. Dittmaier, M. Roth and M.M. Weber, Phys. Lett. B560 (2003) 196; 
hep-ph/0301189 and Nucl. Phys. B660 (2003)289; hep-ph/0302198. 

[4] For a review of some of these techniques and achievements, see Z. Bern, Nucl. Phys. 
(Proc.Suppl). 117 (2003); hep-ph/02 12406. 

[5] E. Witten, Commun. Math. Phys. , 252 (2004) 189; hep-th/0312171. 

[6] A. Brandhuber, B. Spence and G. Travaghni, Nucl. Phys. B706 (2005) 150; hep- 
th/0407214. 

[7] V.V. Khoze, hep-th/0408233. 

[8] I. Bena, Z. Bern, D.A. Kosower and R. Roiban, Phys. Rev. D71 (2005) 106010; 
hep-th/0410054. 

[9] M. Veltman, ''SCHOONSCHIP, A CDC 6600 programme for symbolic evaluation 
of Feynman diagrams." CERN Preprint, 1965. 
H. Strubbe, Comput. Phys. Commun. 8 (1974) 1. 

[10] See for example, A.C. Hearn, Commun. ACM 14 (1971) 511. 
A.C. Hearn, J. Comp. Phys. 5 (1970) 280. 

[11] A.C. Hearn: REDUCE User's Manual, version 3.7, Rand. Corp. 1999. 

[12] J. A. M. VermasereniA^eu; Features of FORM; math-ph/0010025. 

[13] S. Wolfram, The Mathematica Book, Wolfram Media, 5th edition, 2003. 
http : //www . wolf ram . com/. 

[14] http: //www. maplesoft . com/products/maple/. 

[15] T. Sasaki, J. Comput. Phys. 22 (1976) 189. 

[16] See for example, T. Kinoshita ed., "Quantum electrodynamics", World Scientific, 
Singapore, 1990. 

[17] S. Kawabata, Comp. Phys. Commun. 41 (1986) 127; iMd., 88 (1995) 309. 



121 



[18] T. Kaneko, S. Kawabata and Y. Shimizu, Comput. Phys. Commun. 43 (1987) 279. 

[19] D. Perret-Gallix and W. Wojcik ed., ''New Computing Techniques in Physics Re- 
search", Edition du CNRS, Paris, 1990. 

[20] For a review see, R. Harlander and M. Steinhauser, Prog. Part. Nucl. Phys. 43 
(1999) 167; hep-ph/9812357. 

[21] P. Nogueira, J. Comput. Phys. 105 (1993) 279. 

[22] L. Briicher, J. Pranzkowski, and D. Kreimer, Comput. Phys. Commun. 115 (1998) 
140. 

[23] M. Steinhauser, Comput. Phys. Commun. 134 (2001) 335; hep-ph/0009029. 

[24] S.A. Larin, F.V. Tkachov, and J.A.M. Vermaseren, Rep. No. NIKHEF-H/91-18 
(Amsterdam, 1991). 

[25] T. Ishikawa, T. Kaneko, K. Kato, S. Kawabata, Y. Shimizu and H. Tanaka, KEK 
Report 92-19, 1993, GRACE manual Ver. 1.0. 

[26] A. Pukhov et al, "CompHEP user's manual, v3.3". Preprint INP MSU 98-41/542, 
1998; hep-ph/9908288. 

[27] R. Mertig, M. Bohm and A. Denner, Comput. Phys. Commun. 64 (1991) 345. 

[28] J. Kiiblbeck, M. Bohm, and A. Denner, Comp. Phys. Commun. 60 (1990) 165; 
H. Eck and J. Kiiblbeck, Cuide to FeynArts 1.0 (Wiirzburg, 1992); 
H. Eck, Guide to FeynArts 2.0 (Wiirzburg, 1995); 
T. Hahn, Comp. Phys. Commun. 140 (2001) 418; hep-ph/0012260. 

[29] R. Mertig, M. Bohm and A. Denner, Comput. Phys. Commun. 64 (1991) 345. 

[30] T.Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153; hep- 
ph/9807565. 

[31] T. Hahn, Nucl. Phys. Proc. Suppl. 135 (2004) 333; hep-ph/0406288. ibid hep- 
ph/0506201. 

[32] T. Stelzer and W.F. Long, Phys. Commun. 81 (1994) 357; hep-ph/9401258. 
http : //madgraph . hep . uiuc . edu/ index . html. 

[33] J.-X. Wang, Comput. Phys. Commun. 77 (1993) 263, 

J.-X. Wang and Y. Shimizu, "New Computing Techniques in Physics Research IIF 
K.-H. Becks and D. Perret-Gallix ed., World Scientific, Singapore, 1994, 517; 
http: //www. ihep. ac . cn/lunwen/wjx/public_html/. 



122 



[34] H. Tanaka et al, Nucl. lustrum, and Meth. A389 (1997) 295. 

J. Fujimoto et al, Comp. Phys. Commun. 153 (2003) 106; hep-ph/0208036. 

[35] F. Caravaglios and M. Moretti, Phys. Lett. B358 (1995) 332; hep-ph/9507237. 

F. Caravaglios and M. Moretti, Z. Physik C74 (1997) 291; hep-ph/9604316. 

[36] A. Kanaki, C.G. Papadopoulos, Comput. Phys. Commun. 132 (2000) 306; hep- 
ph/0002082. 

[37] M. Moretti, T. Ohl and J. Renter, hep-ph/0102195. 
[38] G. 't Hooft and M. Veltman, Nucl. Phys. 44 (1972) 189. 

[39] C. G. BoUini, J.J. Giambiagi and A. Gonzales Dominguez, Nuovo Cim. 31 (1964) 
551. 

C. G. BoUini and J.J. Giambiagi, Nuovo Cim. 12A (1972) 20. 

J.F. Ashmore, Nuovo Cim. Lett. 4 (1972) 289. 

G. M. Cicuta and E. Montaldi, Nuovo Cim. Lett. 4 (1972) 329. 

[40] A. Denner and S. Dittmaier, Nucl. Phys. B658 (2003), 175; hep-ph/02 12259. 

[41] T. Hahn, LoopTools, http://www.feynarts.de/looptools/. 

[42] F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. 
Shimizu and Y. Yasui, Phys. Lett. B600 (2004) 65; liep-ph/0407065. 

[43] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, Y. Kurihara, K. 
Kato, and Y. Shimizu, Phys. Lett. B576 (2003) 152; hep-ph/0309010. 

[44] Zhang Ren- You, Ma Wcn-Gan, Chen Hui, Sun Yan-Bin, Hon Hong-Sheng, Phys. 
Lett. B578 (2004) 349; hep-ph/0308203. 

[45] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. 
Shimizu and Y. Yasui, Phys. Lett. B571 (2003) 163; hcp-ph/0307029. 

[46] A. Denner, S. Dittmaier, M. Roth and M.M. Weber, Phys. Lett. B575 (2003) 290; 
hep-ph/0307193 and Nucl. Phys. B680 (2004) 85; hep-ph/0309274. 

[47] You Yu, Ma Wen-Can, Chen Hui, Zhang Ren- You, Sun Yan-Bin and Hou Hong- 
Sheng, Phys. Lett. B571 (2003) 85; hep-ph/0306036. . 

[48] Chen Hui, Ma Wen-Can, Zhang Ren- You, Zhou Pei-Jun, Hou Hong-Sheng, Nucl. 
Phys. B683 (2004) 196; hep-ph/0309106. 

[49] F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. 
Shimizu, S. Yamashita and Y. Yasui, Nucl. Instrum. and Meth. A534 (2004) 334; 
hep-ph/0404098. 



123 



[50] F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. 
Shimizu and Y. Yasui, hep-ph/0510184. 

[51] A. Denner, S. Dittmaier, M. Roth and L.H. Wieders, Phys. Lett. B612 (2005); 
hep-ph/0502063. ibid Nucl. Phys. B724 (2005) 247; hep-ph/0505042. 

[52] F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara and Y. 
Shimizu, Nucl. Phys. Proc.Suppl. 135 (2004) 323; hep-ph/0407079. 

[53] K. Fujikawa, Phys. Rev. BY (1973) 393 

M. Base and N.D. Hari Dass, Ann. Phys. 94 (1975) 349 

M.B. Gavela, G. Girardi, C. Malleville and P. Sorba, Nucl. Phys. B193 (1981) 257 
N.G. Deshpande and M. Nazerimonfared, Nucl. Phys. B213 (1983) 390 
F. Boudjema, Phys. Lett. B187 (1987) 362. 

M. Baillargeon and F. Boudjema, Phys. Lett. B272 (1991) 158, ibid B317 (1993) 
371; hep-ph/9308342. 

[54] F. Boudjema and E. Chopin, Z. Phys. C73 (1996) 85; hep-ph/9507396. 

[55] A. Denner, G. Weiglein and S. Dittmaier, Nucl. Phys. B440 (1995) 95; hep- 
ph/9410338. 

[56] For some of these issues as weU as a comprehensive review of radiative correction 
in the electroweak theory see, D. Y. Bardin and G. Passarino, The Standard Model 
in the making: precision study of the electroweak interactions, Clarendon Press, 
Oxford, U.K. (1999). 

[57] H. Tanaka, Comput. Phys. Commun. 58 (1990) 153. 

[58] J. Fleischer and O.V. Tarasov, Comput. Phys. Commun. 71 (1992) 193. 

[59] L.V. Avdeev, J. Fleischer, M.Yu. Kalmykov, and M.N. Tentyukov, Comput. Phys. 

Commun. 107 (1997) 155; hep-ph/9710222. 

M. Tentyukov and J. Fleischer, Comput. Phys. Commun. 132 (2000) 124; hep- 
ph/9904258. 

[60] A. Lorca and T. Riemann, Comput. Phys. Commun. 174 (2006) 71; hep- 
ph/0412047. 

[61] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 
149 (2002) 103; hep-ph/0112278. ibid, Comput. Phys. Commun. in Press; hep- 
ph/0405253. 

http : //wwwlapp . in2p3 . f r/lapth/micromegas/. 



124 



[62] A. Pukhov, hep-ph/0412191. See also, 

http : //www . if h . de/~pukhov/calchep . html. 

[63] A. Semenov. LanHEP — a package for automatic generation of Feynman rules. 
User's manual; hep-ph/9608488. 

A. Semenov, Nucl. Inst. Meth. and Inst. A393 (1997) p. 293. 
A. Semenov, Comp. Phys. Commun. 115 (1998) 124. 
A. Semenov, hep-ph/0208011. 

[64] T. Kaneko, Nucl. Instr. and Meth. A502 (2003) 555. 

[65] C. Becchi, A. Rouet, R. Stora, Comm. Math. Phys. 42 (1975) 127. 

See also, B. W. Lee, Les Houches, Session XXVIII, 1975, Methods in field theory, 
eds. R. Balian, J. Zinn- Justin, North Holland, Amsterdam, 1976. 

[66] I.V. Tyutin, Gauge invariance in field theory and statiscal physics in operator for- 
malism, Lebedev Institute Preprint, N° 39 (1975). 

[67] See for example, R.C. Read, Lecture Notes in Mathematics 884 (1980) 77. 

[68] See for example, C. M. Hoffman, Lecture Notes in Computer Science 136 Springer- 
Verlag, 1982. 

J. Kobler, U. Schoning and J. Toran, Lecture Notes in Computer Science 577 
Springer- Verlag, 1992 p.401. 

[69] T. Kaneko, Comput. Phys. Commun. 92 (1995) 127; hep-th/9408107. 

[70] H. Murayama, I. Watanabe, and K. Hagiwara, KEK Preport 91-11 (1992). 

[71] A. Denner, H. Eck, O. Hahn and J. Kiiblbeck, Nucl. Phys. B387 (1992) 467. 

[72] See G. P. Lepage, J. Comp. Phys. 27 (1978) 192. 

[73] K. Tobimatsu and S. Kawabata, Tech. Rep. 85, Kogakuin Univ. (1998) 423. 

[74] E. de Doncker, A. Gupta, J. Ball, P. Ealy and A. Genz., "Proc. of the 10th ACM 
International Conference on Supercomputing," (1996) 149; 
http : //www . cs . wmich . edu/~parint/. 

[75] W. Kilian, http://www-ttp.physik.uni-karlsruhe.de/whizard/. 

[76] S. L. Adler, Phys. Rev. 177 (1969) 2426. 

J. S. Bell, R. Jackiw, Nuovo Cim. 60 A (1969) 47. 

[77] C. Bouchiat, J. Iliopoulos, P. Meyer, Phys. Lett. 38 (1972) 519. 



125 



[78] W. Pauli and F. Villars, Rev. Mod. Phys. 21 (1949) 433. 

[79] L.H. Ryder, Quantum Field Theory, Cambridge University Press, June 6, 1996. 

[80] M.E. Peskin and D.V. Schroeder, An Introduction to Quantum Field Theory, Perseus 
Books, Frontiers in Physics, March 1993. 

[81] P. Breitenlohner, D. Maison, Commun. Math. Phys. 52 (1977) 39. 

[82] G. Bonneau, Int. J. Mod. Phys. A5 (1990) 3831. 

[83] F. Jegerlehner, Eur. Phys. J. C18 (2001) 673; hep-th/0005255. 

[84] W. Siegel, Phys. Lett. B84 (1979) 193. 

[85] W. Siegel, Phys. Lett. B94 (1980) 37. 

[86] I. Jack and D.R.T. Jones, Regularisation of supersymmetric theories, in Perspectives 
on supersymmetry, p. 149, Editor G.L. Kane; hep-ph/9707278. 

[87] D.M. Capper, D.R.T. Jones and P. van Nieuwenhuizen, Nucl. Phys. B167 (1980) 
479. 

[88] D.Z. Freedman, K. Johnson and J.I. Latorre, Nucl. Phys. B371 (1992) 353. 
P.E. Haagensen, Mod. Phys. Lett. A7 (1992) 893; hep-th/9111015. 
F. del Aguila, A. Culatti, R. Muiioz Tapia and M. Perez- Victoria, Phys. Lett. B419 

(1998) 263; hep-th/9709067. 

F. del Aguila, A. Culatti, R. Mufioz Tapia and M. Perez- Victoria, Nucl. Phys. B537 

(1999) 561; hep-ph/9806451. 

F. del Aguila, A. Culatti, R. Munoz Tapia and M. Perez- Victoria, Nucl. Phys. B504 
(1997) 532; hep-ph/9702342. 

[89] I.J.R. Aitchison, An informal introduction to gauge field theories, Cambridge Uni- 
versity Press, Cambridge, 1982. 

[90] K. Aoki, Z. Hioki, R. Kawabe, M. Konuma and T. Muta, Suppl. Prog. Theor. Phys. 
73 (1982) 1. 

[91] A. Slavnov, Theor. Math. Phys. 10 (1972) 99. 
J. C. Taylor, Nucl. Phys. 33 (1971) 436. 

G. 't Hooft, Nucl. Phys. 35 (1971) 167. 

[92] M. B5hm, H. Spiesberger and W. HoUik, Fortsch. Phys. 34 (1986) 687. 
[93] A. Denner, J. Kiiblbeck, R. Mertig and M. Bohm, Z. Phys. C56 (1992) 261. 



126 



[94] B.A. Kniehl, Z. Phys. C55 (1992) 605. 

[95] J. Fleischer and F. Jegerlehner, Nud. Phys. B216 (1983) 469. 

[96] T. Bhattacharya and S. Willenbrock, Phys. Rev. D47 (1993) 469. 

K. Melnikov, M. Spira and O. Yakovlev, Z. Phys. C64 (1994) 401; hcp-ph/9405301. 
B. A. Kniehl, CP. Palisoc and A. Sirlin, Nud. Phys. B591 (2000) 296; hep- 
ph/0007002. 

[97] F. Kleefeld, AIP Conf. Proc. 660 (2003) 325; hep-ph/02 11460. 

[98] D. Bailin and A. Love, Introduction to Gauge Field Theory, Adam Hilger, Bristol 
and Boston, lOP Pubhshing Limiting 1986. 

[99] D. Espriu, J. Manzano and P. Talavera, Phys. Rev. D66 (2002) 076002; hep- 
ph/0204085. 

[100] B. A. Kniehl and A. Sirlin, Phys. Lett. B530 (2002) 129; hep-ph/0110296. 
[101] M. J. G. Veltman, Physica 29 (1963) 186. 
[102] R. G. Stuart, Phys. Lett. B262 (1991) 113. 

[103] A. Sirlin, Phys. Rev. Lett. 67 (1991) 2127, Phys. Lett. B267 (1991) 240. 
[104] H. G. J. Veltman, Z. Phys. C62 (1994) 35. 

[105] A. Aeppli, G. J. van Oldenborgh and D. Wyler, Nucl. Phys. B428 (1994) 126; 
hep-ph/9312212. 

[106] Y. Kurihara, D. Perret-Gallix and Y. Shimizu, Phys. Lett. B349 (1995) 367; hep- 
ph/9412215. 

[107] E.N. Argyres et al, Phys. Lett. B358 (1995) 339; hep-ph/9507216. 

[108] U. Baur and D. Zeppenfeld, Phys. Rev. Lett. 75 (1995) 1002; hep-ph/9503344. 

[109] W. Beenakker et al, Nucl. Phys. B500 (1997) 255; hep-ph/9612260. 

[110] G. Passarino, Nucl. Phys. B574 (2000) 451; hep-ph/9911482. 

[Ill] W. Beenakker, F.A. Berends and A.P. Chapovsky, Nud. Phys. B573 (2000) 503; 
hep-ph/0303105. 

[112] W. Beenakker et a/., Nucl. Phys. B667 (2003) 359; hep-ph/0303105. 

[113] M. Bcnckc, A. P. Chapovsky, A. Signer and G. Zandcrighi, Phys. Rev. Lett. 93 
(2004) 011602; hep-ph/0312331 and Nucl. Phys. B686 (2004) 205. hep-ph/0401002. 

127 



[114] J. Fujimoto, M. Igarashi, N. Nakazawa, Y. Shimizu and K. Tobimatsu, Suppl. 
Prog. Theor. Phys. 100 (1990) 1. 

[115] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. 
Shimizu, hep-ph/0308080. 

[116] A. Denner et ai, Nucl. Phys. B560 (1999) 33; hep-ph/9904472. 

[117] R. G. Stuart, p. 41, in Proceedings of the XXVth Rencontre de Moriond, edited by 
Tran Thanh Van, Editions Frontieres, Gif-sur-Yvette, 1990. 

[118] R. E. Cutkosky, J. Math. Phys. 1 (1960) 429. 

[119] G. Passarino and M. Veltman, Nucl. Phys. 160 (1979) 151. 

[120] L. M. Brown and R. P. Feynman, Phys. Rev. 85 (1952) 231. 

[121] Z. Bern, L. Dixon and D.A. Kosower, Phys. Lett. B302 (1993) 299; Erratum-ibid. 
B318 (1993) 649; hep-ph/92 12308. 

[122] G. J. van Oldenborgh, Comput. Phys. Commun. 58 (1991) 1. 

[123] J. Fujimoto, Y. Shimizu, K. Kato and Y. Oyanagi, Prog. Theor. Phys. 87 (1992) 
1233. 

[124] W. Beenakker and A. Denner, Nucl. Phys. B338 (1990) 349. 

[125] D. B. Mehose, // Nuovo Cimento 40A (1965) 181. 

W.L. van Neerven and J.A.M. Vermaseren, Phys. Lett. 137 (1984) 241. 

[126] P. Nogucira and J.C. RomSo, Z. Phys. C60 (1993) 757 . 

[127] T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon and C. Schubert, JHEP 10 (2005) 
015; hep-ph/0504267. 

[128] A. Denner and S. Dittmaier, Nucl. Phys. B734 (2006) 62; hep-ph/0509141. 

[129] F. Yuasa et al.„ Proceedings AIHENP 99,hep-ph/0006268. F. Yuasa, D. Perret- 
GaUix, S. Kawabata and T. Ishikawa, Nucl. lustrum, and Meth. A389 (1997) 77. 

[130] J. Fujimoto, N. Hamaguchi, T. Ishikawa, T. Kaneko, H. Morita, D. Perret-Gallix, 
A. Tokura and Y. Shimizu, talk given by J. Fujimoto at ACAT05, May 2005, DESY, 
Zeuthen, Germany. 

http : //www-zeuthen . desy . de/acat05/talks/Fuj imoto . Junpei . 3/ 
talk_f uj imoto-1905 . ppt. 



128 



[131] R. G. Stuart, Comput. Phys. Commun. 48 (1988) 367. 

R. G. Stuart and A. Gongora, Comput. Phys. Commun. 56 (1990) 337. 

R. G. Stuart, Comput. Phys. Commun. 85 (1995) 267. 

G. Devaraj and R.G. Stuart, Nud. Phys. B519 (1998) 483. 

[132] F. Boudjema, A. Semenov and D. Temes, Phys. Rev. D72 (2005) 055024; hep- 
ph/0507127. 

[133] F. del Aguila and R. Pittau, JHEP 07 (2004) 017; hep-ph/0404120. 

[134] W. Giclc, E. W. N. Glover and G. Zanderighi, Nud. Phys. Proc. Suppl. 135 (2004) 
275; hep-ph/0407016. 

R. K. Ellis, W. T. Giele and G. Zanderighi, liep-pli/0508308. 

[135] W. T. Giele and E. W. N. Glover, JHEP 04 (2004) 029; hep-ph/0402152. 

[136] A. van Hameren, J. VoUinga and S. Weinzierl, Eur. Phys. J. C41 (2005) 361; hep- 
ph/0502165. 

[137] A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, Nud. Phys. B650 (2003) 
162; hep-ph/0209219. 

[138] E. de Doncker, Y. Shimizu, J. Fujimoto and F. Yuasa, Comput. Phys. Commun. 
159 (2004) 145. 

E. de Doncker, Y. Shimizu, J. Fujimoto, F. Yuasa, K. Kaugars, L. Cucos and J. Van 
Voorst, Nud. Instrum. Meth. A534 (2004) 269; hep-ph/0405098. 

[139] T. Binoth, G. Heinrich and N. Kauer, Nud. Phys. B654 (2003) 277; hep- 
ph/0210023. 

[140] Y. Kurihara and T. Kaneko, hep-ph/0503003. 
[141] C. Anastasiou and A. Daleo, hep-ph/0511176. 

[142] G. Duplancic and B. Nizic, Eur. Phys. J. C35 (2004) 105; hep-ph/0303184. 

[143] Y. Kurihara, hep-ph/0504251. 

[144] http: //minami-home . kek.jp/ jf /nig. html. 

[145] F. Bloch and A. Nordsicck, Phys. Rev. 52 (1937) 54. 

D.R. Yennie, S.C. Frautchi and H. Suura, Ann. of Phys. 13 (1961) 379. 

[146] C. Itzykson and J.B. Zuber, Quantum Field Theory, McGraw-Hill Education, 
September 1980. 



129 



A. Calogeracos, N. Dombey and A. D. Kennedy, Phys. Lett. (1981) B104 444 . 
J. Prenkel, R. Meuldermans, I. Mohammad and J. C. Taylor, Phys. Lett. B64 
(1976) 211 . 

R. Gastmans and R. Meuldermans, Nucl. Phys. B63 (1973) 277. 
W.J. Marciano and A. Sirlin, Nucl. Phys. B88 (1975) 86. 

G. 't Hooft and M. Veltman, Nucl. Phys. B153 (1979) 365. 

Y. Kurihara, J. Fujimoto, T. Munehisa and Y. Shimizu, Prog. Theor. Phys. 96 

(1996) 1223; hcp-ph/9603312. 

T. Munehisa, J. Fujimoto, Y. Kurihara and Y. Shimizu, Prog. Theor. Phys. 95 
(1996) 375; hep-ph/9603322. 

J. Fleischer, A. Leike, T. Riemann, and A. Werthenbach, Eur. Phys. J. C31 (2003) 
37, hep-ph/0302259. 

J. Fujimoto and Y. Shimizu, Mod. Phys. Lett. 3A (1988) 581. 

W. Beenakker, S. van der Marck, and W. HoUik, Nucl. Phys. B365 (1991) 24-78. 

J. Fleischer, J. Fujimoto, T. Ishikawa, A. Leike, T. Riemann, Y. Shimizu, and 
A. Werthenbach, in Second Symposium on Computational Particle Physics (GPP, 
Tokyo, 28-30 Nov 2001), KEK Proceedings 2002-11 (2002) (Y. Kurihara, ed.), 
pp. 153-162; hep-ph/0203220. 

J. Fleischer, T. Hahn, W. Holhk, T. Riemann, G. Schappacher, and A. Werthenbach, 
contribution to the second workshop of the extended EGFA/DESY study Physics 
and Detectors for a 90 to 800 GeV Linear Collider, 12-15 April 2002, Saint-Malo, 
France; hep-ph/0202109. 

For a nice review on the radiative corrections to e+e~ W^W~ and related issues, 
see W. Beenakker and A. Denner, Int. J. Mod. Phys. A9 (1994) 4837. 

M. Lemoine and M. Veltman, Nucl. Phys. B164 (1980) 445. 

M. Bohm, A. Denner, T. Sack, W. Beenakker, F.A. Berends and H. Kuijf, 
Nucl. Phys. B304 (1988) 463. 

J. Fleischer, F. Jegerlehner and M. Zralek, Z. Phys. C42 (1989) 409. 

A. Denner, S. Dittmaier and M. Strobel, Phys. Rev. D53 (1996) 44; hep- 
ph/9507372. 

A. Denner and S. Dittmaier and R. Schuster, Nucl. Phys. B452 (1995) 80; hep- 
ph/9503442. 



130 



[160] G. Jikia, Nucl. Phys. B494 (1997) 19; hep-ph/9612380. 

[161] A. Denner and S. Dittmaier, Nucl. Phys. B398 (1993) 239. 

[162] A. Denner and S. Dittmaier, Nucl. Phys. B398 (1993) 265. 

[163] A. Denner and T. Hahn, Nucl. Phys. B525 (1998) 27; hep-ph/9711302. 



131 



